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Some Aerodynamic Problems of Aircraft 
Engines 


WILLIAM REDE HAWTHORNE* 
Massachusetts Institute of Technology 


INTRODUCTION 


In the creation of a new machine, man’s inventiveness often 
outstrips his understanding. It is surprising how often tech 
nological advances have been made from a state of knowledge 
which in retrospect we would be tempted to regard as cripplingly 
inadequate. Someone has remarked that thermodynamics owes 
more to the steam engine than the steam engine to thermody 
namics; and it is certainly true that the engines which powered 
the industrial revolution were designed without the benefit of the 
first and second laws of thermodynamics. We are told in fact 
that it was the existence of the steam engine which inspired Carnot 
to puzzle over the problem of determining the highest efficiency 
one could get from an engine, and led to the formulation of the 
second law of thermodynamics 

A new machine may therefore not only fulfill a need and cause 
social change, but it can also pose questions in men’s minds which 
lead to new understanding and even to the formulation of pro 
found natural laws. This new understanding in turn provides 
fertile ground for new inventions 

Nowhere is this leapfrogging of invention and knowledge better 
illustrated than in the history of aviation. The idea of flight 
has inspired many inventions and much research. Two inven 
tions, the Handley-Page Slot and the swept-back wing, both of 
Which are almost dismayingly simple ideas, illustrate how the 
understanding resulting from research makes new inventions 
possible. On the other hand, the elegant theory of induced drag 
illustrates the process of improving understanding under the 
pressure of development 

Ina lecture in honor of Minta Martin, whose interest in educa 
tion has led to the establishment of memorial funds for the sup 
port of education and research in aeronautics, it seems appro 
priate to emphasize and try to describe the growth of understand 
ing in the field of engineering with which the lecturer has been 
associated. Hence, while my choice of subject is the jet engine, 
my treatment of the subject will concentrate on basic ideas and 
their application rather than on specific details of its design 

I shall concentrate on some of the aerodynamical aspects of the 
engine, discussing a few of the new ideas or new applications of 
Old ideas which have been important in its development. These 


Based on the First Minta Martin Lecture presented at the 
University of Maryland, March 21, 1956. Revised and received 
May 28, 1957 

* Jerome Clarke Hunsaker Professor of Aeronautical Engineer 
ing (Visiting Now, Professor, Engineering Departinent, Uni 
versity of Cambridge, Cambridge, England 


include ideas on the flow in ducts, the flow in turbines and com 
pressors, some aspects of combustion, and flow with density 
change. A large number of aspects of the aerodynamics of gas 
turbines are naturally omitted by this selection, not because they 
are unimportant It has not been possible, for instance, to dis 
cuss heat transfer, supersonic compressors, jet noise, blade 
flutter, or intakes and supersonic diffusers, in spite of their 


recognized importance 


THE INFLUENCE OF COMPRESSIBILITY 


A IMPRESSION of some of the features of turbojet 
engines may be gained from Figs. 1, 2, and 3. The 
diagram in Fig. | shows the Mach Numbers in the flow 
at various points in an early Whittle-type engine. In 
many places the Mach Numbers are of the order of 
unity, showing the effort made by the designer to reduce 
frontal area, size, and weight. The tip speed of the 
double-sided centrifugal compressor at nearly 1500 
ft./sec. was considerably higher than that of normal 
industrial compressors and led to a large Mach Number 
at inlet to the diffusers. The choice of a double-sided 
impeller was made in order to reduce the Mach Number 
at the tip of the impeller inlet which even then was at 
the fairly high value of 0.89 relative to the rotating 
vanes. By this means the air swallowing capacity of 
the compressor per unit of frontal area was made 
satisfactorily high, but the air flow through the engine 
became limited by the Mach Numbers and stresses in 
the turbine blading. The author has shown the interre- 
lation between the design choices and stress and Mach 
Number limitations in turbojet engines with centrifugal 
compressors elsewhere.*? In fact, increases in air mass 
flow and thrust output were mainly achieved in these 
engines by lengthening the turbine blades. 

Fig. 2 shows the pressures, temperatures, and Mach 
Numbers at various points in a later turbojet engine 
with an axial compressor. The general layout of the 
single rotor engine shown in Fig. 2, often with a com 
pressor of about 12 stages and a single stage turbine, is 
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Fic. 1. Typical Whittle-type turbojet engine showing values 


of Mach Number at various points for designed r.p.m. sea level 
test bed conditions. 





» Ces résultats bizarres ont déja été signalés en partie sous une autre 
forme par M. Coriolis (*). On ne les évite pas en faisant varier le nombre 
fini m. On n'y échappe pas non plus en tenant compte du refroidisse- 
ment de lair, du a la dilatation quill éprouve jusqu’’ son passage par 
lorifice; car si l'on suppose, avec Lagrange et M. Poisson, que la pres- 
de la densité J, l'équation 1) 
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(*) Rapport fait le 19 feévrier 1838 sur divers Memoires de M. de Pambour. [es 
developpements verbaux que M. Coriolis a bien voulu donner, alan des auteurs , sur 
un passage de ce Rapport (Compie rendu, 1* semestre , 1838, p. 22g), ont ete la prim- 
cipale ovcasion de ces recherches. 


Fic. 4. Reproduction of St. Venant and Wantzel’s formula for 
compressible flow through an orifice.! P, A stagnation pressure 
and density, P; downstream pressure, y ratio of specific heats, 
m a constant, and V, volume flow rate per unit area based on 
stagnation density. 
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e 3 ¢ § |} 
STAGNATION PRESSURE, psia} 45 | 13.8 | 54.5 | 51.8 | 67 
STAGNATION TEMPERATURE,°F |-40 | 150 | 500 {1600 | 1122 
STAGNATION AXIAL VELOCITY | 548 | 481 | 407 | 490 | 920 


Fic. 3. Typical two-spool turbojet showing values of Mac 
Numbers at various points at design conditions. Table show 
typical values for a 12.1 pressure ratio engine. 
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TOTAL PRESSURE...... 14.7 TOTAL PRESSURE.__-_ 8| 


TOTAL TEMPERATURE..518 TOTAL TEMPERATURE..9 40 


TEMPERATURE RISE...1O60 


TOTAL PRESSURE...78 TOTAL PRESSURE...-..- 33 


TOTAL TEMP....2000 TOTAL TEMPERATURE. ..|675 


AXIAL MACH NO......... O44 AXIAL MACH NO._____. 0.28 INLET VELOCITY.---.... 110 AXIAL MACH NO......---.047 
REL. MACH NO..........089 PRESSURE RATIO___.5.5 REL. MACH NO... 0.85 PRESSURE RATIO__....-.24 
(REL. TO! ROTOR BLADE) (REL. TO BUCKET) 
UNITS 
PRESSURE ......... PSIA 
TEMPERATURE - - --- 
VELOCITY. Ft/Sec. 


Typical turbojet with axial compressor showing 


conditions at various points in engine 
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enerally typical of the engines which have largely re- 
nlaced those with a centrifugal compressor in the last 
decade. The lower frontal area of the axial compressor 
and the ease with which higher pressure ratios can be 
achieved by adding stages makes the higher Mach 
Numbers less essential to the attractiveness of the en- 
gine as a basic type than in the case of the engine based 
on a single centrifugal compressor. There is still, how- 
ever, the competitive drive to increase Mach Numbers 
and thereby reduce frontal areas and the number of 
stages required. 

As performance and design techniques have improved, 
the opportunity to improve fuel consumptions by in- 
creasing the pressure ratio has been seized. To over- 
come difficulties caused by the large departure of the 
velocities from their design values, most noticeable in 
the first stages of compressors operating at part loads, a 
two-spool design, shown diagrammatically in Fig. 3, has 
been evolved. As might be expected in an engine de- 
signed for low fuel consumption, the Mach Numbers 
have been maintained on the whole at moderate values. 
High subsonic Mach Numbers have, of course, been 
used in turbojet engines to reduce size simply because 
the smallest cross-sectional area of a duct passing a 





given mass flow is obtained when 7 = 1. Increasing 
the circumferential velocities enables the work per stage 
to be increased, and, hence, the weight to be reduced. 
The concentration on Mach Numbers in the flow 
through jet engines is not an incidental matter. Early 
workers (Whittle, for instance) were aware of the im- 
portance of compressibility. One-dimensional gas dy- 
namics has in fact been reorganized and developed to 
become the principal analytical tool of the engine aero- 
dynamicist. Earlier still (in 1894) the Swedish engineer 
De Laval patented a converging-diverging nozzle for use 
in turbines, and textbooks on steam turbines written 
at the close of the last century deal with compressible 
flow in nozzles and blades, although not in terms of 
Mach Numbers, a parameter which only came in to 
extensive use in the 1930's. 

We need to go back to the beginning of the Nineteenth 
Century to trace the growth of understanding of com- 
pressible flow in ducts which has been so important in 
the development of the jet engine. Fig. 4 is a reproduc- 
tion of a portion of a paper by St. Venant and Wantzel 
in the Proceedings of the French Academy of Sciences 
in 1839.! It shows an equation which has now become 
familiar to all those who learn gas dynamics. P is the 
stagnation pressure, A is the stagnation density, P; is 
the static pressure downstream of an orifice, y is the 
ratio of specific heats, m is a constant, |’, is the volume 
flow rate per unit area based on the stagnation density. 
As the footnote shows, the research described in this 
paper was instigated by Coriolis. 

It was well known in hydraulics at the time that the 
rate of flow of a liquid through an orifice of sufficient 
size Was proportional to the square root of the difference 
of pressure across the orifice. On deriving a formula for 
the flow of an elastic fluid through an orifice—a formula 


such as the one shown in Fig. 4—it was noted that, as 
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the downstream pressure, P,, was reduced below a cer 
tain value, the flow given by the formula diminished. 
This appeared to suggest the ridiculous conclusion that 
gas would not flow through an orifice into an evacuated 
chamber! 

St. Venant and Wantzel conducted experiments with 
an orifice discharging into an evacuated chamber. 
They found that as long as the downstream pressure P, 
did not exceed between 0.3 and 0.4 of the upstream 
pressure P, which was fixed, the rate of flow through 
the orifice was unaffected by the downstream pressure. 
They remarked in their paper that there is every reason 
to think that the average pressure in the orifice never 
falls below 3/5 of the upstream pressure. We now 
know that at approximately this critical pressure the 
velocity of sound is reached in the orifice and hence 
changes in pressure downstream are not transmitted 
through the orifice, so that the flow into and through the 
orifice is unaffected. The term choking is often applied 
to this condition. 

The formula in Fig. 4 relates the flow area and pres- 
sure ratio for the one-dimensional flow of a perfect gas 
Remarkably enough it was derived before 
either the first (1834) (1850) laws of 
thermodynamics were satisfactorily formulated. Apart 


in a duct. 
or the second 


from the pressure-density relation, which was also used 
by Laplace in his derivation of the formula for the 
velocity of sound (1816), the derivation depended only 
upon the laws of mechanics from which Bernoulli's equa- 
tion is derived. Zeuner® points out that St. Venant and 
Wantzel’s formula and results remained unnoticed for 
about 30 years, perhaps because Poncelet argued that 
such an odd result must be due to the smallness of size 
of the orifices used, the largest being only 0.059 inches 
in diameter. 

Kloster® rediscovered St. Venant’s work in 1867, and 
he and Rankine’ derived the formula for the critical 
pressure ratio. The fact that the velocity in the throat 
is equal to the local velocity of sound appears to have 
been recognized first by Holtzmann in 1861,? although 
Osborne Reynolds* and Hugoniot# clarified the whole 
question in 1885 and 1886. 

The development of analytical methods and the pro- 
duction of numerical tables for treating the one-dimen- 
sional flow of a perfect gas were accomplished by a 
number of workers. The external and internal aerody- 
namicists have a common interest in isentropic one- 
dimensional flow,'’ but flow with friction and heat addi- 
tion is normally of greater importance to the engine 
aerodynamicist. Examples of these latter flows will be 
discussed below since they contain some unusual fea- 
tures which were only noticed as a result of wartime and 
postwar aircraft engine research and development. A 
systematic analysis in a form convenient for use in in- 
ternal aerodynamics was given by Shapiro and Haw- 
thorne*’ in 1947. The value of such an analytical tool 
is, of course, largely dependent on the assumption that 
the flow in ducts may be regarded for most purposes as 
one-dimensional. This assumption appears to be justi- 
fied by results. 
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AN EXAMPLE OF CHOKING IN FRICTIONAL 
COMPRESSIBLE FLOW 


When the velocity of the gas in a duct reaches the 
velocity of sound in such a way that no downstream 
effects can be transmitted upstream, the condition is 
described as “choking.”’ The most familiar case of 
choking is that described earlier with reference to flow 
through a nozzle. <A set of turbine nozzle vanes choke 
due to the attainment of the velocity of sound at the 
throats. In a turbine it is possible to have choking 
occur in three places simultaneously —first, in the throat 
of the nozzles, second, in the throat at outlet from the 
rotating blades, and third, in the annulus downstream 
when the exhaust pressure is so low that the axia/ ve- 
locity in the annulus is the velocity of sound. In this 
condition the turbine is aerodynamically ‘‘locked”’ and 
further reduction in the exhaust pressure at outlet from 
the downstream annulus will not increase the mass flow 
or the work, nor will it alter the velocity, pressure, or 
temperature anywhere in the blading. Although friction 
modifies these effects, they would occur in the absence 
of friction as a result of the contraction and expansion 
of stream tubes in the compressible flow. 

Choking also occurs in the flow in pipes of uniform 
bore owing to the effect of friction. As the subsonic 
flow proceeds along the pipe, the pressure decreases 
owing to wall friction. This causes a decrease in 
density, which in turn increases the velocity, since con- 
tinuity must be satisfied. The temperature also de- 
creases and the Mach Number increases. 

The effects of friction in compressib’e flow were dis- 
cussed by Grashof and Zeuner, who introduced the 
polytropic exponent or small stage efficiency. Grashof,* 
Zeuner,? and Fanno® showed that the velocity in 
adiabatic frictional flow down a pipe of uniform bore al- 
ways approaches the velocity of sound. The analysis*® 
gives the following differential equations for pressure Pp, 
temperature 7’) and Mach Number .\/, change, in terms 
of the coefficient of friction f, length x, and hydraulic 
mean diameter of pipe D, and the ratio of specific heats 


k. 


dp 7 —kMP?[1 + (Rk — 1)4A?] ip 
p 2(1 — M?) D 
dT 7 —k(k — 1)A3 pe 
T 2241— M) D 


and df? = RM2)1 + [((k — 1)/2]477} ' dx 
Mm (1 — M?) Ip 


from which it will be seen that in subsonic flow pressur 
decreases and Mach Number increases in the directioy 
of flow (x increasing). In supersonic flow the pressur 
increases and Mach Number decreases. It is interesting 
to observe that in subsonic flow the effect of friction js 
to decrease the temperature of the stream. The equa- 
tions show that as J/ approaches unity the Mach Nun. 
ber changes rapidly in a short length of pipe. The 
equations restrict the point at which a Mach Number of 
unity can be obtained to the downstream end of a pipe 
of uniform bore where the cross section is subsequently 
sufficiently enlarged, as, for instance, in a tube discharg. 
ing into the atmosphere. It will be noted that the 
choking is not due to an effective reduction in area due 
to increasing displacement thickness of the boundary 
layer on the walls of the pipe, since fully developed pipe 
flow is assumed in the analysis. 

Hence, when the pressure in the space into which a 
pipe of uniform bore discharges is sufliciently reduced 
the velocity of sound is reached at the discharge of the 
pipe. No further reduction in downstream pressure can 
then have an effect on the flow in the pipe. The Mach 
Number at outlet will be unity, and the Mach Numbers 
at upstream sections of the pipe will have characteristic 
choking values depending on the coefficient of friction 
and the distance from the discharge end. This de- 
pendence is shown in Fig. 5. In subsonic flow the Mach 
Number cannot exceed the choking values. 

This phenomenon is now well understood, but it is oi 
interest to recall how its significance was emphasized to 
some of us in the early British turbojet research and de 
velopment. In one experimental turbojet the cross- 
sectional area of the ducting between the turbine outlet 
and the propelling nozzle was maintained constant. 
To ensure this result the sheet metal exhaust cones were 
especially faired and curved shapes. The Mach Num 
ber at outlet from this particular turbine and hence at 
the inlet to the duct was designed at nearly 0.8 (com- 
pare the value of 0.67, Fig. 1). Fig. 5 shows that with 
this inlet Mach Number, a length of only 4 duct diame- 
ters is required to reach a Mach Number of unity. The 
duct was longer than this and the choking Mach Num- 
ber at inlet to the exhaust system was therefore appre 
ciably lower than 0.8. The engine was faced, there- 
fore, with the difficult problem of discharging the gas 
from the turbine into the exhaust pipe at a lower Macli 
Number. This meant either a lower exhaust velocit\ 
or higher temperature (to increase the corresponding 
velocity of sound). One way to solve the problem 0! 
maintaining the same approximate mass flow was to it- 
crease the density, i.e., the pressure of the gases leaving 
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the turbine. This reduced the available pressure drop 
through the turbine and necessitated an intolerable 
increase in temperature to obtain sufficient work from 
the turbine to drive the compressor. 

The high temperature was the first intimation that 
mything was wrong. It was found that maximum 
turbine temperatures were reached well below the de- 
jgned r.p.m. Checks showed that there were no exces- 
sive losses and the use of a larger or no propelling nozzle 
gave insufficient improvement since the choking due to 
iriction Was so large. 

The reason for the high temperature was soon dis- 
wvered* and exhaust ducts were subsequently designed 
to have increasing cross-sectional area and some diffu- 
sion downstream of the turbine (see Fig. 1, which shows 
a decrease in Mach Number. from 0.67 to 0.47 in the 
exhaust duct). This made the fabrication of exhaust 
ducts easier since simple conical sections gave the de- 
sired diffusion without any difficulty. 

This was an example of a situation in which ex- 
perience with an engine directed attention to earlier 
theoretical work and stimulated its revision in more 
modern terms. 

Another example of frictional choking which oc- 





curred in early turbojet engines, and probably still 
occurs in some, is in impulse turbine blade passages. If 





the blades are closely enough pitched they may be re- 
garded as ducts of constant cross-sectional area. They 
do not generally tolerate subsonic inlet Mach Numbers 


greater than 0.7. 


AXIAL COMPRESSORS 


Axial flow fans have almost as long a history as wind- 
mills, but the axial compressor with its successive rows 
of closely pitched blades is a much more recent develop- 
ment. Parsons started to build axial compressors in 
1901, the first having an efficiency of about 60 per cent,"! 
probably, according to Constant,'* because the blades 
were stalled, but it was not until the twenties that the 








development of this type of compressor to efficiencies 
high enough to use in an effective gas turbine really 
began. Whereas a quasi-one-dimensional approach was 
satisfactory in turbine design, the successful develop- 
ment of axial compressors has needed the application of 
theories of the isolated airfoil, and their modification to 
suit the different conditions in the cascade. 

Papers published in the late twenties (references are 
given by Keller,'® Tyler,'' and Howell”) record the 
application of airfoil theory to propellers and fans and 
the development of theories of the flow through two- 
dimensional cascades. Various methods were em- 
ployed for considering the effects of neighboring blades. 
Betz,'® for instance, calculated corrections by replacing 
neighboring blades by vortices. The work of Weinig™ 
tepresents an example of the theoretical development in 
the mid-thirties. 

More recent work has made it possible to solve the 
*It was explained by J. Reeman, then at the Royal Aircraft 
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potential flow problem in cascades, at least for incom- 
pressible flow, by a variety of methods for both the 
direct problem—i.e., when the profile is predeter- 
and the inverse problem when the pressure or 


Some approximate 


mined 
velocity distribution is specified. 
methods are available for compressible flows, al- 
though Lin’s method** is available for theoretically 
exact results. While mathematically, at least for in- 
compressible flows, the potential flow problem is solved, 
the results of the application of the techniques has been 
relatively disappointing. Sheer computational difficulty 
may partly account for this and for this reason there has 
been some use of the electrolytic tank. But the main 
difficulties appear to be due to factors which make the 
characteristics of compressor blades different from those 
of isolated airfoils. At reasonable Mach Numbers pro- 
file shape does not appear to have much influence on 
compressor performance, perhaps because its influence 
cannot easily be separated from other factors such as 
stage matching and three-dimensional and end effects. 
At the normal Reynolds Number and turbulence level 
in compressors, transition of the boundary layer on the 
profile probably occurs very near the leading e ge. 
Hence control of the boundary layer by altering the 
profile shape has fairly limited potentialities. Also, in 
a compressor cascade the blades are pitched so closely 
together that it is possible to consider them as the 
boundaries of bent diverging passages. The passage 
area ratio and overall diffusion is determined by the 
flow inlet and outlet angles. The “‘passage’’ approach, 
which has been used in turbine design for years, has not 
lost all usefulness and is complementary to the airfoil 
approach. Howell" showed in his summary of cascade 
test data in 1942 that losses began to increase when the 
diffusion exceeds a certain value dependent on the 
pitch-chord ratio and flow outlet angle and that the 
influence of camber and blade profile was small. 

The idea of stall has been transferred from the iso- 
lated airfoil to cascades, although in the latter case it is 
the increase in total pressure loss rather than the loss of 
lift which is significant. Howell’ defined the stall as 
occurring at the incidence at which the total pressure 
loss doubled its minimum value. A designed or nominal 
deflection was then defined as 0.8 of the stalling deflec- 
tion. This corresponds roughly to the condition of maxi- 
mum lift/drag ratio and to reasonably satisfactory 
pressure distributions. 

An axial compressor stage, however, is generally said 
to stall when the pressure rise reaches a maximum value. 
A similar criterion could also be adopted for cascade 
stalling. There is at present no completely accepted 
definition of stall which covers both cascades and com 
pressors, although the maximum pressure rise criterion 
has been most generally accepted for the compressor 
stage. 

The development of a method of predicting “stall” 
based on “‘passage’’ theory has been continued by 
+ Blade thickness and camber line shape affect the performance 


at high Mach Numbers 
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Fic. 6. Flow in annulus between two rows of blades; cu, 
tangential velocity; c:, axial velocity; 7, radius; x, axial coordi- 
nate. 


workers at the NACA. Meanwhile, the potential flow 
solutions have been used to confirm cascade pressure 
distributions and turning angles. The prediction of 
stall has also been studied by Goldstein and Mager*‘ and 
others. On the whole, however, potential flow theory 
has not yet shown the way to any great reduction in loss 
or improvement in operating range, and in spite of its 
help in predicting the onset of some compressibility 
effects, the theory has not yet been of as much help as 
it has to the external aerodynamicist. 


*THREE-DIMENSIONAL FLOW IN COMPRESSORS AND 
TURBINES 


The flow in an axial compressor or turbine is not a 
simple two-dimensional flow. The blades are disposed 
in an annulus with the object of producing large changes 
in the angular momentum of the fluid about the 
rotor axis. The radial pressure gradient is generally 
large enough to produce considerable differences in the 
Hence veloci- 
ties and flow angles vary with radial distance. Also, the 
alteration of the angular momentum of the fluid as it 


flow from one radial position to another. 


passes through a blade row introduces changes in the 
radial pressure gradient from one side of the blade row 
to the other which may cause radial flows within the 
blades, as well as in the annular spaces between the 
blade rows. 

These effects were not at all understood in early steam 
turbine practice in which attention was concentrated on 
the behaviour of single nozzles or small groups of blades 
as though they were discharging into free air. 

About 1930, a number of workers* unknown to each 
other began almost simultaneously to apply aerody- 
namic theory not only to two-dimensional cascades but 
also to the three-dimensional flow in the annular ar- 
rangement of the blades. 

The main development in this period was the use of 
constant circulation blading. Ina potential flow, blades 


with the circulation at radius shed no 
vorticity into the stream, so that the potential flow is 


It was thought 


same every 
retained through succeeding blade rows. 
that the absence of vorticity would minimize induced 
effects and yield higher efficiencies. 
If no circulation is shed, the circulation enclosed in 
\. A. Griffiths advocated the use of constant circulation 
blading in 1926 and built a small experimental unit in 1929 (see 
Howell®). O 
turbine blades with constant circulation in 1931. F 


Westinghouse designed long steam 
Whittle de- 
There 


Tietjens at 


veloped free-vortex blading for his jet engine in 1937.%° 


may have been others. 
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any circumferential circuit drawn in the annulus by 
tween two rows of blades and enclosing the rotor ayj 
remains constant. Hence, 27rc, = 
and the tangential velocity c, varies inversely as th 


constant (Fig, | 


radius. 
blade rows with each particle of fluid possessing th 


This gives a free-vortex flow between th 
same angular momentum. Hence, as each particle ; 
fluid passes through a rotating blade row, it suffers th; 
same change of angular momentum, and, consequently 
the work done by the blades is the same at every radius 
Since the input (or output) of work is uniform over th 
annulus, the stagnation pressure and enthalpy is als 
uniform, and a little analysis shows that the axial y¢ 
The flow wit! 
It is the 
basis for the frequently used free-vortex design. | 
leads to stator blades which are relatively untwisted 


locity, c,, is also invariant with radius. 
rc, and c, both invariant with radius is simple. 


but in the rotor blades a considerable twist is required 
from root to tip. The pressure change across the rotor 
blades is least at the hub or root and highest at the tip 

When the flow may be regarded as incompressible, or 
at least the variation of density p, with radius is negli 
gible, then, in the free vortex design with c, invariant 
with radius, continuity requires that the streamlines lic 
on cylindrical surfaces—i.e., if p and c, are invariant 
with radius, then pc, is invariant with radius and there 
is no radial flow. (It is assumed that the inner and 
outer radii of the annulus containing the blades is sub- 
stantially constant.) 

In these conditions Euler's equations for frictionless 


incompressible flow with axial symmetry reduce to 
(1/p)(Op/Or) = ¢,?/r 
1/p(Op/Ox) = c,(Oc,/Ox) 


Now in the more general case of frictionless, com 
pressible flow with axial symmetry, Eq. (1) becomes 


(1/p)(Op/oOr) = (c,°/r) — 
cC,(Oc,/Or) — c,(OcC,/Ox) (2 


which reduces to Eq. (1) when the radial velocity, « 
If the blade rows are far enough apart 
with radius, and 


dies out to zero. 
whatever the variation of p and c, 
whether the blading is free vortex or not, c, and Oc, 0 
will tend to zero between the blade rows and Eq. (! 
will apply. Eq. (1) is called the condition for ‘‘radial 
equilibrium,’ and it is commonly assumed in nonfret 
vortex and compressible flow designs that this condition 
obtains at the inlet and outlet from the blade rows 
i.e., that any radial velocities and displacements of the 
flow occur well within the blades so that at the leading 
and trailing edges c, and Oc,/Ox are negligible. 

The application of free-vortex design methods and the 
condition of radial equilibrium may in some cases hav 


interesting consequences. In a free-vortex high-speed 
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single-stage turbine, for example, owing to the large 
pressure gradient in the annular gap between nozzles 
and blades, the density, p, of the gas at the root is con- 
If the axial velocit) 
the mass 


siderably less than it is at the tip. 
is constant, the specific mass flow—i.e., pc;, 
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fow per unit area -between the nozzles and the blades 
isless at the root than it is at the tip. Consequently, 
radial velocities are developed in the flow through the 
nozzles and blades, as is shown approximately in Fig. 7. 
Owing to the radially inward flow in the blades, con- 
tinuity is sometimes theoretically satisfied at the root 
{ the moving blades by making the blade passage 
larger at the exit than at the inlet. Several turbines 
were built with blade passages which, if regarded as 
two-dimensional cascades, would be classed as com- 
pressor blades. The root sections had blade passages 
which diverged in the direction of flow (Fig. 7) with 
the minimum cross-sectional area at the inlet,* while, 
wing to the inward radial flow in the rotating blade 
passages, a small pressure drop across the blades was 
whieved. The application of this surprising result of 
the radial equilibrium assumption gave satisfactory re- 
sults in a number of cases. However, when the radial 
velocities are large, it appears to be impossible to 
achieve all the required radial displacement of the 


streamlines within the blade passages. 


AcTuATOR Disc THEORY APPLIED TO COMPRESSORS 


The radial equilibrium assumption, which has been 
extensively adopted in compressor and turbine design, 
has proved satisfactory when radial variations of density 
and circulation are not large. 

*I am indebted to F. Schwoerer of the A. G. T. Division, 
Westinghouse Electric Corporation, for the information that some 
large, low-pressure steam turbine blades are currently designed 


in this fashion 
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Fic.7. Flow through a free-vortex turbine showing meridional 
streamlines and blade passage and flow angles at the root (hub). 
Displacement of streamlines shown twice actual value. Veloci 
ties are multiples of critical velocity in upstream flow 
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Meridional streamlines in the flow through a single 


stage compressor.” 
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For aircraft engines the free vortex or constant circu- 
lation design in the compressor has the disadvantage 
that it gives somewhat large frontal areas, if losses in 
efficiency due to high Mach Number either at the tip of 
the rotor or root of the stator are to be avoided. It ap- 
pears to be more reasonable to equalize the Mach 
Number over the radius by avoiding the extreme varia- 
tion of c, with radius implied by the free-vortex law, 
To achieve this result, some compres- 


i. constant. 
sors have been built in which the circulation around 
some blades has varied six-fold or more between hub 
and tip. 

In some turbines, attempts to avoid radial displace- 
ment of the streamlines have been deliberately made by 
introducing the design requirement that pc, is inde- 
pendent of radius. In these ‘constant specific mass flow 
designs” there is also a radial variation in circulation. 

One of the earliest papers dealing with the effect of 
variable circulation on the flow in a compressor is that 
of Ruden.*! Merchant”? in 1946 and Marble** in 1948 
applied the actuator disc concept of propeller theory to 
a row of compressor or turbine blades. In the actuator 
disc method it is assumed, as in propeller theory, that 
the blades are replaced by a thin disc which may be re- 
garded as an infinite number of blades of vanishingly 
The disc is normally located at the center 
Across the dise continuity is 
The change in 


small chord. 
of pressure of the blades. 
satisfied but losses may be introduced. 
angular momentum across the dise is determined by 
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using the flow angles obtained from two-dimensional 
cascade measurements. The flow angles so obtained 
are assumed to be those at the trailing edges of the 
blades and are imposed on the flow at the appropriate 
distance downstream from the disc. 

More exact theoretical solutions for the flow in an- 





nular rows of blades have been attempted, notably by | © suction s 
Wu.*! But the more exact the solution, often the more q| Me rolls into 
tedious its application, and since for many purposes the “ > observer 
“radial equilibrium’? method appears satisfactory, it = 7 The ut 
remains to be seen whether anything more accurate z S sional be 
than the actuator disc method is necessary. S 2 yelopmet! 
The writer, Horlock,* and others have, therefore, 4 ~ by the t 
developed the method for use in multistage compres- ” ao was init: 
sors. For design calculations it gives some improve- cades an 
ment over the radial equilibrium method when radial cascades 
variations in circulation are large. It is perhaps most In the 
useful in the calculation of the flow in blading off the dynamic 
design point, for at some conditions of operation, even 
in the most carefully designed stage, the radial velocities 
can become appreciable. This point is illustrated in where V 
Fig. 8, which shows the measured meridional stream- vorticity 
lines in a single-stage axial compressor, consisting of a ind p th 
row of guide vanes, rotor blades and stator blades. : —_s . The A 
Fig. S(a) shows a section of the surfaces containing the ghtCjjll 4, Stk, wuaization of secondary, ean a] a of by 
streamlines at the design point (full lines) and ap- on bottom wall near inlet. It passes toward and up suction | becomes 


: . - : surface and, separating the wall, spirals mstreat 
proaching the stall (broken lines). Fig. S(b) shows the icc d, separating from le Wall, spira downstream 


flow just after the stall with the rotor tip and stator root 
stalled. It will be noticed that the flow ‘‘avoids’’ the , 
ment of the flow toward the root has unstalled the 


stalled regions at the rotor tip and stator root and that 
? stator root. curl V X 


Fig. 9 shows the observed axial velocity variation 


and the | 


the radial velocities are large in consequence. Fig. S(c) 
is at an even lower flow rate, and shows that the rotor 


ie ; i with radius at various stations for the badly stalled 
tip is now badly stalled and that the resulting displace- : 


case shown in Fig. S(c). The broken lines show the} The thir: 





predictions of the actuator disc theory already referred | ‘Url v= 
to. The actuator discs were placed at the center of | pressible 
== 100 pressure of the blade rows (for rotors and stators at the | ¢quation 


2a ax 
he hk. ee 
— 
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leading edge). Approximate blade performances wert 
obtained from published two-dimensional cascade tests 
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and all losses were assumed to occur across the actuator Suppose 







100 discs. The inviscid equations of motion between the 
GUIDE — discs were linearized, and the equations for each row oi | “iV VX 
blades solved by an iterative procedure. The method 
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SOF = is given in reference 25 and will shortly be published a 
XN ed ec 
elsewhere. 
ROTOR TRAILING \ fe iets ileat is " rod wit | otter 
OF EDGE 100 t seems like y cmat ¢ le actuator disc met 10¢ = ies in 
enable the location of stall to be determined, and 
er te consequently it may be of help in the design of com- 
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150 |~ > 50 pressors for optimum performance over a wide operating ie attend 
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x range. 
IOOL STATOR TRAILING NX ~~-4 9 
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50+ —e— OBSERVED It has been recognized for some time that the losses in : 
\ compressors and turbines are only partially accounted * grad 
---- PREDICTED . ne , glint the 
‘ for by the profile and skin-friction drag of the blades 
04 05 O6 O7 O8 O89 10 rhe clearance gaps which have to be maintained be | po. io, 
R= RADIUS tween rotating and stationary parts are another sourct eer 
TIP RADIUS of loss. The boundary layers which build up on the} gq. << a1, 
Fic. 9. Radial variation of axial velocity in the flow through walls of the annulus and in the region where the ends o/ 


a badly stalled single-stage axial compressor. Comparison be- a : 
tween theory and experiment .2 the blades meet the walls also are a source of loss and 
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distortion of the flow. An illustration of the complex 
three-dimensional flows occurring in a blade passage as 
aresult of the distortion of the boundary layer on the 
end walls by the action of the blades is given in Fig. 10. 
Smoke admitted through a hole on the bottom wall near 
the entrance to the passage sweeps toward and up the 
suction surface of the blade, then leaves the surface and 
rolls into a vortex which winds downstream (toward the 
observer ). 

The understanding of the flow in these three-dimen- 
sonal boundary layers has been assisted by the de- 
velopment of an inviscid theory, which was prompted 
by the turbine and compressor problem. The theory 
was initially developed by Squire and Winter” for cas- 
cades and has since been applied to bends, struts, and 
cascades by Hawthorne.** ** 

In the inviscid incompressible flow of a fluid hydro- 
dynamic theory shows that 


V X Q = grad po p 


where V is the velocity vector, and Q = curl V is the 
vorticity vector, po is the total or stagnation pressure 
and p the density. 

The Helmholtz equation is obtained by taking the 
curl of both sides of this equation. The right-hand side 


becomes identically zero—hence, 





curl V KX Q = 0 (3) 
and the left-hand side 
curl V K Q = (Q-grad)V — 

(V-grad)Q + Vidiv Q) — Q(div V) 


The third and fourth terms are zero since div Q = div 
curl V = 0 and, by continuity in the flow of an incom- 
pressible fluid, (div V) = 0. Hence, the Helmhotz 


equations are obtained, viz., 
(Q-grad)V = (V-grad)Q 
Suppose that we now examine the following identity: 


divV X (VX Q) = (VX Q)-curl V — 
V-curl (V X Q) 


we find that, since curl V = Q, the first triple product 
on the right-hand side is zero and by Eq. (3) the second 
term is also zero. Hence, since 
VX (VV X 2) = (a-V)V — (V-V)Q 
we obtain the result 
div (Q:-V)V = div (V-V)Q 

or 
V- grad)(Q-V) + (Q-V) div V = 

(Q-grad)(V-V) + (V-V)(div Q) 


For incompressible flow the second term on the left- 
hand side is zero and the second term on the right-hand 
side is also zero by the definition of Q. Hence, 


(V-grad)(Q-V) = (Q-grad)(V-V) (4) 
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The left-hand side of this equation expresses the 
growth of the product of velocity and the component of 
vorticity in the direction of flow along a streamline. 
The right-hand side is the product of the magnitude of 
Q and the derivative of (V-V) in the direction of Q. 
If we imagine a flow approaching a bend or a cascade, 
then near the wall the vorticity is perpendicular to the 
velocity. Hence (Q-V) is zero initially. As the flow 
passes round a bend, a gradient of (V-V) in the bend is 
developed normal to the direction of flow and hence in 
the direction of the vorticity. Consequently the right- 
hand side of Eq. (4) becomes other than zero and (Q-V) 
Therefore a component of vorticity in the 
Squire and Winter” have 


increases. 
direction of flow develops. 
shown that for a simple bend of deflection, €, the vor- 
ticity in the direction of flow becomes approximately 
2€Q; where Q) is the vorticity initially perpendicular to 
the flow. 

It may be shown that for three-dimensional inviscid 
flow in bends and cascades, the growth of a vorticity in 
the direction of flow is of great importance since it im- 
plies the growth of secondary velocities transverse to 
the main flow, which will, as the flow continues, ap- 
preciably distort and twist the original flow. 

For the flow through cascades it is possible to follow 
kinematically the distortion of the vortex filaments. 
The flow is assumed to approach the blades with a 
velocity gradient normal to the end wall, and viscous 
effects are neglected in the subsequent flow. Hence the 
approaching flow contains vortex filaments which are 
normal to the flow. As the flow carries each filament 
through the row of blades (Fig. 11) one portion will 
travel at a different speed from another so that the 
filament no longer remains normal to the flow. The 
distortion of the filament is shown in Fig. 11. Starting 
at ab it passes through the blades to def. From a to e 
along the suction surface of the blade the filament trav- 
els faster than it does from 6 to f along the pressure 
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Fic. 11. Flow of a vortex filament through a cascade 
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Fic. 12. Components of vorticity in direction of flow at out- 


let from cascade. 
distortion of vortex filaments in the approaching wall boundary 
layer. (B) Trailing vortex sheet containing the vorticity from 
the circulation shed off the blade and the stretched vortex fila- 
ments from the approaching wall boundary layer. (C) Trailing 
vortex sheet outside the boundary layer containing the vorticity 
from the circulation shed off the blade only. (D) Streamlines 
of the secondary flow induced by A in the passage between the 
blades and the trailing vortex sheets. (FE) Cross hatching shows 
the region in which separation frequently occurs. 


surface. Hence the portion ef is inclined somewhat as 
shown and develops a component of vorticity in the 
direction of flow, the “distributed secondary vorticity.”’ 
The portion of the filament near a which passes along 
the pressure surface of the upper blade in Fig. 11 only 
reaches d when the nearby portion has reached e by 
passing along the suction surface. Hence a section of 
the filament lies stretched from d to e in the streamline 
trailing from the blade. As may be seen from Fig. 11, 
it is of opposite sign to the distributed secondary vor- 
ticity. In addition to these streamwise components of 
vorticity, there will be a component due to the change in 
circulation along the blades also lying in the trailing 
streamline. The various components of vorticity are 
shown in Fig. 12 for an example in which the velocities 
induced by them are not sufficient to distort the flow 
The picture of the smoke filament (Fig. 
Experi- 


appreciably. 
i0) shows a considerable distortion of the flow. 
ments by Hawthorne and Armstrong”® show that for 
thick boundary layers, the main features of the flow may 
be determined by assuming that the distributed vortic- 
ity induces secondary velocities within a 
bounded by the blades and the trailing vortex sheets 


passage 
leaving them. Streamlines of the secondary velocities 
are sketched in Fig. 12. There is evidently a change in 
the direction of the spanwise secondary velocity across 
the trailing edge of each blade. The magnitude of this 
change gives the strength of the trailing vortex sheet. 
At spanwise positions outside the wall boundary layer, 


(A) Distributed secondary vorticity from the 
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the trailing vortex sheet contains shed circulation only 
so that the secondary flow induces a decrease in circula. 
tion and hence a reduction in lift and deflection even at 
positions on the blades which are not immersed in the 
wall boundary layer. Inside the wall boundary layer 
the trailing vortex sheet contains both the shed circula. 
tion and the trailing vorticity due to stretching of the 
vortex filaments in the approaching flow as they pass 
through the blades (portions such as de, Fig. 11). 
These effects give rise to a spanwise variation of lift 
coefficient and outlet angle along the blade such as js 


shown in Fig. 13. The lift coefficient at each spanwise 
position is based on the approaching velocity appro 
priate to that position. Outside the wall boundary 
layer, where the approaching velocity is uniform, the 
lift coefficient and deflection decrease as the wall is ap- 
proached. Inside the boundary layer the trend is re. 
versed until a large lift coefficient is found near the 
wall. The increase in lift coefficient near the wall found 
by Hawthorne and Armstrong’ may be regarded as due 
to an upwash; consequently, a negative induced drag 
may be expected. Measurements by Mair” on an 
isolated airfoil in a stream with a spanwise velocity 
gradient, but no wall, have shown that an increase in 
lift coefficient and reduction in drag take place more or 
less as predicted. The inviscid theory is, however, in- 
adequate to account for all the effects which occur, and 
considerable flow separation is normally observed in 
the corner between the end wall and the suction surface 
(Fig. 12). 

The theory of the development of the secondary flows 
based on the assumption that the flow in the cascade is 
inviscid is not satisfactory. An alternative approach 
suggested by Mager®® and developed by Herzig and 
Hansen’! is based on the use of the boundary-layer ap- 
proximations. It is possible that further development 
will include a combination of the inviscid theory with 
the boundary-layer theory developed by the latter 
authors. 
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Fic. 13. Effect of secondary flow on outlet angle and lift 





coefficient. Compressor cascade inlet angle = 50°, normal 
outlet angle = 30°, normal lift coefficient = 0.92, s/c = 1.0, 
boundary-layer thickness/spacing = 0.43. 
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AERODYNAMIC 


PROPAGATING STALL 


Aircraft axial compressors have to operate over such 
4 wide range of rotational speeds that stalling of some 
of the stages for appreciable periods is inevitable. A re- 
markable aerodynamic feature of a stalled compressor 
js the phenomenon of stall propagation or rotating stall. 
The proper description of this phenomenon was first 
given by Emmons, Pearson, and Grant* in 1953. Pre- 
viously, Whittle and Smith and Fletcher (see reference 
299) had observed patches of reversed flow which regu- 
larly changed their circumferential position during the 
stalled operation of a compressor. By the use of hot 
Emmons and his co-workers were 
“stall 


not only on a centrifugal and axial compressor 


wire anemometers, 
able to detect moving patches of stalled flow 
cells”’ 
but also on a two-dimensional cascade. 

The phenomenon may be explained qualitatively by 
relerence to Fig. 14, which shows a two-dimensional 
cascade operating at an incidence which is almost suf- 
ficient to cause stall. <A slight disturbance to the inlet 
flow may cause one blade A to stall which results in 
The 
partial blockage of the passage diverts air toward blades 
Cand D. 
quently increased, and that on blade C reduced. Blade 
B therefore tends to stall and, if the incidence changes 
are sufficient, stall will propagate in the direction ABD. 


separation of the flow in passage AB as shown. 


The incidence on blades B and D is conse- 


Eventually passage A will become unstalled so that a 
cell of several stalled blades will propagate in the direc- 
tion shown. If the blades are on a rotor, the propaga- 
tion will be in the opposite direction to the rotation. 
Observations on stalled axial compressors show that 
one row of blades may have several stall cells propagat- 
ing simultaneously, and that stall propagation on one 
row of blades may occur simultaneously with stall 
propagation on others. Fig. 15 shows an instantaneous 
picture of stall cells propagating simultaneously on a 
rotor and The 
shaded regions at the tip are two stall cells on the rotor, 


the following stator (see Wood’). 
propagating in this case at 62 r.p.s. The shaded regions 
near the root are two stall cells at the stator root propa- 
gating at 24 r.p.s. The rates of propagation quoted 
in both cases are relative to a stationary observer. 

A schlieren picture and an interferometer picture 
taken by the M.I.T. Gas Turbine Laboratory) of stall 


propagation in a cascade are shown in Fig. 16. It is 
possible to observe vortices of opposite sign being shed 
Irom the blades as the angle of incidence alters. Stall 


propagation in cascades has been observed over a four 
degree range of incidence. 

There is now a considerable literature on the experi- 
mental observation of stall propagation sources 
such as NACA, M.I.T., C.1I.T., N.G.T.E., and various 
On the theoretical side there are also a 


>of 
ez. 


companies. 
number of papers in which the problem has been tackled. 
It is perhaps premature to emphasize one particular 
author’s solution in favor of another. Emmons,” 
Sears, Marble,®! and Kriebel® have developed theories. 


More recently Rannie®* and Fabri® have developed 
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~! 
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OF AIRCRAFT ENGINES 
small perturbation theories which predict the speeds of 
propagation when the disturbances are not large. But 
the major difficulty is due to the difference in some 
cases between the velocities of propagation measured 
in a cascade wind tunnel and the performance in a com- 
pressor. 

Stall propagation is undoubtedly a phenomenon of the 
greatest importance to the engine aerodynamicist and 
is one in which our understanding is still in the pre- 
liminary stage. 


COMBUSTION 


In the early stages of the development of turbojet 
engines, the development of the combustion chamber 
posed the most formidable problem. It was required to 
operate at heat release rates an order of magnitude 
greater than those previously met with in industrial 
Complete utilization of the fuel and low 
The combustion had to 


practice. 
losses in pressure were desired. 
be stable and efficient over a much wider range of mix- 
ture strength and inlet and flow conditions than had 
hitherto been required. The products had to be well 
mixed so that a desirable temperature distribution 
would be obtained in the gases entering the turbine. The 
problems set by these requirements, the difficulties, 
and the development of suitable combustion systems 
have been described elsewhere (see Whittle,*’ Lloyd,” 
Watson and Clarke*'). The general features of the flow 
in a typical modern combustion chamber are shown in 
Fig. 17.°° The Figure shows a primary zone in which 
stabilization of the flame and preliminary combustion 
is achieved and a dilution zone for completion of com- 
bustion and mixing of the excess air with the products. 

In turbojet combustion chambers, the total air-fuel 
mixture is too lean to permit premixing of all the air and 
fuel, so that fuel is injected into a portion of the air- 
stream in which it mixes and burns. Dilution air is 
added later. Premixing of fuel and air before combus- 
tion is used in ram-jets (and afterburners) where the 
air-fuel ratio is within the limits of inflammability of the 
mixture. 

In the early development of combustion chambers 
there was considerable attention paid to the mixing 
process. It was realized that the rate at which combus- 
tion proceeded depended on the rate at which fuel and 
oxygen molecules were brought together. Attention 
was, in other words, directed toward the aerodynamic 
and convective features of the flow rather than to the 
chemical kinetics, which are microscopic rather than 
macroscopic phenomena. The development of this 
idea may be traced from the work of Groume-Grji- 
iailo*® who made liquid models of the gas flow in fur- 
naces, to the work of Burke Schumann**! who 
applied a molecular diffusion theory to the burning of 
laminar jets of fuel gas in air streams, and to the work 
at M.I.T. in 1937® on the combustion of turbulent jets 


and 


of gas. 
Numerous experiments on the means of mixing in 
order to achieve both rapid combustion and dilution 
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AERODYNAMIC 


were conducted at the Royal Aircraft Establishment 
during the war. Attempting to devise methods of in- 
creasing the rate of mixing, Minchin® found that when 
, circular cylinder was placed at the end of a baffle 
separating streams of hot and cold air flowing in a duct 
Fig. IS), an almost uniform temperature profile could 
be achieved in a short distance (profile A). Further 
iownstream, however, an apparent unmixing occurs, 
the temperature on the hot side becoming less and that 
on the cold side greater than the average (profile B, Fig. 
\s). Complete mixing had not in fact taken place at A, 
the temperature measured being the statistical average 
of evenly distributed eddies or vortices which were in 
motion from one side of the duct to the other. 

In flames intimate mixing between molecules is re- 
quired to complete the combustion. It is not achieved 
as soon as a Statistical average of the concentration of 
juel and air reaches the stoichiometric proportions. 

The idea that the aerodynamic flow pattern con- 
trolled the combustion rate, because it controlled the 
mixing of fuel and air, was very helpful in the early de- 
velopment of combustion systems. In more recent ex- 
perience it has been found in afterburners that the 
mixing process is equally if not more important since it 
tends to be more difficult to control. 


The actual flow in combustion chambers is so complex 








that it defies theoretical analysis except of the crudest 
type. There has, however, been a considerable study of 
nonburning air and water models and the development 
{ mixing devices, many of which have been adapted for 
use in practical combustion chambers. 

Another aspect of the combustion problem which has 
received great attention is flame stability. The main- 
tenance of a flame depends on the diffusion of heat and 
probably certain active radicals to the reactants enter- 
ing the flame zone. In gas turbine combustion cham- 
bers, flame propagation is controlled by heat and mass 
transfer in a complex flow. Broadly speaking, stability 
is achieved in practice by recirculating hot or burning 
gas so that it heats up the incoming reactant mixture. 
Suitable recirculatory flows have been obtained by 
hours of cut and try development testing in combustion 
chamber test rigs. 

We have obtained some understanding of the prob- 
lem of flame stability by studying relatively simple 
systems. One basic problem, which has been exten- 
sively studied, has been that of stabilizing a flame in a 
premixed mixture of fuel and air. Scurlock,* Long- 
well, Barrére and Mestre, and Spalding” have con 
sidered the problem of flame stabilization in the wake 
of rods and other baffles. The phenomenon is illus- 
trated in Fig. 19. In the Figure the broken lines show 
a boundary-layer region between the wake and _ the 
main stream. In the wake immediately behind the 
stabilizer rod there is a region of intense mixing in 


which no unburned gas appears.'! Further downstream 





a well defined recirculation is found which is supplied 
with burned or burning gas from the turbulent interface 
between the wake and the main stream. If the Reyn- 
olds Number based on rod diameter is low enough—.e., 
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lower than about 10‘—the boundary layer between the 
wake and main stream is laminar for a distance behind 
the baffle, becoming turbulent further downstream. 
Since the rates of molecular diffusion across the laminar 
boundary region are not the same for the different 
chemical species, the mixture in the immediate wake 
will have a different atomic composition from that of the 
approaching stream. For many fuel gas-air mixtures, 
the composition of the wake is thereby shifted toward a 
leaner mixture than that in the main stream.* 
mixtures are then required in the approaching stream 


Richer 


to obtain stability. 

As the Reynolds Number is increased, the laminar 
boundary region shortens, and eventually disappears 
altogether, probably when the critical Reynolds Num- 
ber for the body is reached. The mixing, which is then 
by turbulent diffusion, occurs at approximately the 
same rate for all species so that the atomic composition 
of the gases in the wake is about the same as that in the 
main stream. 

The presence of a laminar boundary rezion not only 
aifects the mixture range for stability, but also the 
velocity at which blowout occurs. Relow the critical 
Reynolds Number the blowout velocity tends to be 
proportional to the square root of the stabilizer width; 
at higher Reynolds Numbers it varies directly as the 
stabilizer width. The difference is attributed by Zu- 
koski and Marble* to the effect of the laminar boundary 
layer on the downstream length of the recirculation zone, 
The 
length of the recirculation zone is important because, if 
it is too short, the time available to complete the reac- 
If mixing in the turbulent region is 


and hence on the time available for combustion. 


tion is too small. 
toc rapid, the temperature may decrease and the reac- 
tion rate be slowed down. 

Research stimulated by the jet engine has covered 
many more aspects of the combustion process than the 
two briefly mentioned here. Fuel injection, droplet 
burning, the effect of turbulence on flame speeds and un- 
steady or oscillatory burning have been studied in- 
tensively. They are all phenomena controlled at least 
partially by the aerodynamic features of the flow. 


FLow WitHh Heat ADDITION 


While the behavior of flames is affected by the flow in 
which they burn, it is also true that the density changes 
occurring during combustion will alter the character of 
the flow. The study of the effect of heat addition in fast 
moving gases has been almost entirely initiated by the 
development of the new aircraft engines. Due recog 
nition should, however, be given to the work of Hugo 
niot, Chapman,** Jouguet,’* and others on the dy- 
namics of moving flame fronts, and to the work on the 
condensation shock when steam becomes supersaturated 
during expansion in a nozzle. During the development 
of the jet engine, interest in flows with heat addition 
started with the problem of combustion chamber pres- 


* For fuel gas of low molecular weight—e.g., methane—the 


shift is toward richer mixtures 
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frictionlessly caused by density changes in heating and cooling. 


sure loss. Aerodynamicists are (or were) normally ac- 
customed to regard losses in pitot or stagnation pres- 
sure as evidence of irreversibility*—e.g., skin friction, 
Part of the 
stagnation pressure drop in a combustion chamber is, 


Let us con- 


viscous dissipation in eddies or shocks. 


however, not due to irreversible effects. 
sider a one-dimensional frictionless flow in a duct of 
uniform cross section, Fig. 20, in which in region X the 
density is decreased, for instance, by adding heat to the 
fluid. Following Rayleigh*® and applying the equations 
of momentum and continuity across the surface X we 
find 


[pr + ('/2) Vi?) — [be + ('/2)p2Ve?] = 


& 2) prVi"7 | (py p2) — 1] 


which, when the velocities are not large, gives the change 
in pitot or stagnation pressure. Provided the flow is 
frictionless, the result is the same whether the heat is 
added reversibly or not, and we may imagine processes 
by which small amounts of heat are added reversibly in 
the region X and later removed reversibly in a region Y 
further downstream (Fig. 20). If the amount of heat 
removed at Y is the same as that added at X, the state 
at (3) will the same as at (1), Fig. 20, and the stagnation 
pressure change across Y will be equal and opposite to 
that at X. As the heat exchange is performed re- 
i.e., With infinitely small temperature dif- 
the surroundings are also restored to their 


versibly 
ferences 
initial state, so that the change from (1) to (2) or from 
(2) to (3) is a thermodynamically reversible one. 

The simplified analysis given above may be extended 
to the flow of a perfect gas with stagnation temperature 
7). It may be shown* that the change in stagnation 
pressure, fo, in frictionless flow along a duct of uniform 


cross-section is given by 
dpo/ po = —(k/2)M?(dTo/T») (5) 


where .\/ is the local Mach Number of the flow and k the 
ratio of the specific heats. This reversible change in 
stagnation pressure decreases to zero as the velocity 
and Mach Number are reduced. The change in Mach 
Number due to stagnation temperature change on heat- 


* Meteorologists were, however, aware of the implications of 
the effects of heat addition in low-speed flows in the atmosphere 
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ing or cooling is given by 


dM? (1+ M*)\1 + [(k — 1)/2]M?} dT, 
Me (1 — M2) y 


Owing to the presence of the (1 — J/?) term, the 
effect of heat addition (7» increase) is therefore to jp. 
crease the Mach Number in subsonic flow and reduce jt 
As M approaches unity, the Mach 
Number changes rapidly for a small change in stagna- 
The equations restrict the point at 


in supersonic flow. 


tion temperature. 
which a Mach Number of unity can be obtained to the 
downstream end of any heat addition process. This 
implies a choking phenomenon in which in subsonic 
flow the Mach Number upstream of any heat addition 
can never exceed the value which gives a Mach Number 
of unity downstream. The elucidation of the problem of 
the pressure loss and choking phenomenon? as a result 
of heat addition has been a significant contribution to 
the theory of compressible flow arising from study of 
the aircraft gas turbine. Its implications are similar to 
those of choking by area constriction or friction. The 
flow upstream of a region of heat addition adjusts its 
pressure or mass flow rate so that the Mach Number 
ahead of the region does not exceed the choking value. 
The temperature rise with hydrocarbon fuels is so large 
that choking Mach Numbers tend to be small—ie., 
from 0.2 to 0.3 
‘duct in which combustion occurs has to be large. This 


so that the cross-sectional area of the 


sets important practical limitations in engine design 
and determines the frontal area of ram-jets and after- 
burners. 

Any loss in stagnation pressure due to friction and 
viscous effects must be added to the reversible pressure 
drop, and the effect of the combined loss allowed for 
in determining the inlet condition for choking.*® 

The variation of temperature due to change in stagna- 
tion temperature by heat interchange is given by 


dT (1—kM*)}1+[(R—1)/2]/ M2} dT 
‘i (1 — AW?) To 


It was pointed out by A. A. Griffith that between VM = 
1//k and M = | the effect of heat addition was to re- 
duce the temperature of the gas. The fact that cooling 
of a gas can be achieved by adding heat to it or, as 
shown earlier, by frictional effects has not escaped most 
of us who value the use of the paradox in teaching as a 
technique for stimulating thought! 

The possibility shown by Eq. (5) of a gain in stagnation 
pressure by heat removal and density increase is of in- 
terest, although the reversible gain is inevitably re- 
duced or lost completely by frictional effects which 
accompany any method for increasing the density of the 
stream. Engineers have used one application of this 
namely, the steam jet injector—for a 
The injector 


phenomenon 
hundred years as a boiler feed pump. 
takes steam at boiler pressure, and by using the re- 


+ See Hawthorne and Cohen. In this reference the effect of 


injecting and mixing another gas stream with the flow was also 


studied 
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versible gain in stagnation pressure obtained by mixing 
cold water with a high velocity jet of steam, so that it 
condenses, raises the water to a pressure above that of 
the boiler. A gain of stagnation pressure of 3 to 4 per 
cent has recently been achieved by evaporating water 
injected in fine streams of droplets” into a fast-moving 
stream of hot gas. If such a device were ever to find 
application, it would be a perfect illustration of how 





understanding reached in the development of one in- 
vention gives rise to new inventions. 

When the flame front* is inclined at an angle to the 
direction of flow, the effect is somewhat similar to the 
fow through an oblique shock wave. Across the flame 
front, however, the deflection is toward rather than 
away from the normal. There are also upstream effects 
since, in the subsonic flow ahead of the flame front, 
disturbances can be transmitted upstream. Fig. 21 
shows the direction and magnitude of the deflection, 6, 
jor various angles of obliquity, 6, assuming that the 
i.e., the velocity normal to the flame 


fame speed 
For each value of the ratio of stagna- 


front—is small.* 
tion temperature of the burned gas, 7\», to the stagna- 
tion temperature of the unburned gas, 7\, there is a 
maximum deflection. When the flame front is normal 
to the flow, the Mach Number always increases. For 








an oblique flame front, only the normal component of 
velocity is altered, so that the more oblique the front 
the less the change in velocity. The velocity of sound, 
however, always increases with the temperature so that 
there is an obliquity at which no change of Mach Num- 
ber occurs across the flame. When the flame is more 
oblique than this, the Mach Number downstream is 
less than that upstream. The value of the angle @ at 
which the changeover occurs is the same as that for 
maximum deflection and for small flame speeds is shown 
in Fig. 21. 

* The term flame front is used to describe the region in which 


combustion occurs in a stream of fuel and air. In laminar flow 


the flame is thin. It is not so thick in turbulent flow that the 


concept loses validity. 
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Fic. 22. Photographs of streamlines passing in laminar flow 
through an oblique hydrogen flame 





Deflection across oblique flame front for small flame 
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exact solution by numerical calculation.** (b) Bottom 
velocity choking flow, linearized solution.” 


Fig. 21 shows that the deflections across oblique 
flame fronts can be large. A large deflection is seen in 
Fig. 22 which is a schlieren photograph of an oblique 
flame burning in a laminar stream of hydrogen and air. 
The streamlines are made visible by magnesium pow- 
der. 

When the flame front is curved, vorticity is created 
in the downstream flow as behind a curved shock. 
Tsien*® has treated this elegantly in a paper dealing 
with the problem of the behavior of a flame stabilized 
by a rod placed in the center of a duct. Flame fronts 
in such practical applications are oblique at angles less 
than the limiting angles shown in Fig. 21, so that the 
Mach Number generally decreases locally across the 
flame. Buta flame stabilized in a duct of uniform cross 
section, even though it is inclined to the flow, will still 
cause a Mach Number increase at the outlet from the 
duct and tend to choke the inlet flow. 

In an early paper, Scurlock® analyzed by approxi- 
mate means the flow through a flame front stabilized 
at the center of a two-dimensional channel. Fig. 23(a) 
shows an exact numerical solution by Emmons, Ball, 
and Maier* of the problem when the velocity is so low 
that the pressure changes do not affect the density. 
The velocity of the flame normal to the flow (the flame 
speed) is assumed constant, and the vorticity generation 
due to curvature of the flame front is included. Al 
though the unburned and burned gas streams both flow 
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in the axial direction far upstream and downstream of 
the flame, there must be a deflection across the flame 
front, so that the gas bends away from the axis just 
upstream of the flame and toward the axis immediately 
downstream. Since the streamlines are bent just up- 
stream of the flame, vorticity appears in the flow down- 
stream. The converging streamlines upstream and 
downstream of the flame indicate an increase of velocity. 
Hence, for a constant flame speed normal to the flame 
front, the flame must become more obliquely inclined 
to the axis as the flow proceeds. 

A linearized solution for higher velocities, from the 
work of Fabri, Siestrunck, Foure, “ 
density is affected by the pressure, is shown in Fig. 23(b) 


and when the 
and is generally similar to the incompressible flow solu- 
Over almost all of the flame front the Mach 
Number decreases and, in fact, the flame front is so 


tion. 


oblique to the flow just upstream that there is prac- 
tically no change in velocity and pressure across it. All 
pressure changes and velocity increases occur in the 
flow.up and downstream of the flame. 

The convergence of the streamlines in the unburned 
gas and the decrease in Mach Number across the flame 
imply that the unburned gas velocity may reach the 
velocity of sound just ahead of the flame front near the 
point where it reaches the wall of the duct. The linear- 
ized theoretical solutions by Tsien“ and Fabri, et al.‘ 
give the same choking inlet Mach Number as the simple 
one-dimensional treatment but show that supersonic 
velocities are attained in the unburned gases near the 
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Fic. 24. Flow through stovepipe ram-jet. (a) Top. With 
heat addition. (b) Center. With cooling. (c) Bottom. Effect 
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wall and in the burned gases near the axis. The Conse. 
quences of the supersonic flow field near the wall hay¢ 
not yet been resolved. 

The problems of heat addition are not separable from 
those of combustion. The spreading of a stabilized 
flame is both a mixing and flame propagation problem 
and a heat addition problem. It seems, however, that 
our understanding has grown by dealing with each 
aspect of the problem separately before considering the 
interaction between the separated parts. 


SOME REMARKS ON THRUS1 


Propulsion is achieved by so altering the pressure dis- 
tribution about an aircraft, or parts of it, that enough 
force in the direction of flight is developed to overcome 
its skin-friction drag and any component of its weight, 
In the past it was possible to some extent to separate the 
activities of the engine maker from those of the airplane 
manufacturer. We could distinguish between the posi- 
tive force which it was the responsibility of the propul- 
sion man to produce on the blades of the propeller with 
a minimum consumption of fuel, and the negative force 
on the rest of the aircraft, which the airplane man 
sought to reduce as much as possible. The new develop- 
ments in aircraft propulsion are making this distinction 
less easy, and even though there are still determined 
attempts to distinguish between thrust and drag,” it is 
likely that the distinction will eventually disappear al- 
together. The forced nature of the distinction may be 
illustrated by considering the flow through a simple 
“stovepipe”’ ram-jet consisting of a hollow thin tube of 
uniform cross section, Fig. 24(a). In this simple case 
we assume that heat is added frictionlessly in the ram- 
jet so that p, is less than p;. If we assume the flow is 
inviscid, the velocities are small, and the Kutta condi- 
tion applies at the trailing edge, the force acting on the 
thin-walled shell is a thrust proportional to the product 
of (V2 — V) and the mass flowing through it. 

Ahead of the region of heat addition, Bernouilli’s 
equation yields 


Pa + (1 2) pV? = Pr + (1 2) pV," 


where p, is the atmospheric pressure and the other 
nomenclature is shown in Fig. 24. 

Applying this result to the equations in Fig. 20 and 
noting that po = py, 


(1/2) pV? — (1/2) p.V.? = (1/2) p1Vi?[(p1/ p2) — 1] 


Using the continuity relation, Fig. 20, p,V; = p2V2, we 
eliminate V; and obtain the result 


(V' V2)? = (p2/ pi)[2 — (p2/ px) ] (6) 


shown graphically in Fig. 24(c) 

Hence we conclude that as ps is less than p;, V2 is 
greater than V and a thrust is obtained. The thrust 
must appear on the body as a pressure force and the 
only place it can act is at the leading edge, where, as a 
result of a negatively infinite pressure acting on an in- 
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fnitely thin body, a thrust is theoretically possible. 
The streamlines are shown approximately in Fig. 24(a). 

The satisfactory development of this thrust is nor- 
mally the responsibility of the propulsion team, who 
concern themselves primarily with the flow passing 
through the body. The thrust in this ideal case is, how- 
ever, developed by the infinite negative pressure at the 
infinitely thin leading edge. Any failure to develop the 
thrust due to improper rounding and cambering of the 
leading edge could perhaps be said to introduce an ex- 
ternal drag since it occurs in the external flow, which is 
normally the concern of the aircraft team. 

In the same ram-jet it is possible at least theoretically 
to develop a thrust when p, is greater than p), in other 
words by abstracting heat from the flow. If p»/ p; is, say, 
39 then (VV) = +/3 2, so that Vz is greater than V, 
and consequently a thrust is produced. As before, the 
thrust is developed by an infinite negative pressure at 
the infinitely thin leading edge, since nowhere else can 


the thrust be sustained. 





Now, from Eq. (6) and the continuity equation in 


Fig. 20, 
(V/V,)? = (V/V.)2(V2/V,)? = 2(p:/p2) — 1 (7) 


so that, as p»’ p; becomes greater than unity, so (V, V) 














becomes greater than unity, and the ram-jet takes in 
more air. Hence, for the ‘‘reverse’’ ramjet, 1.e., the one 
with po > p;, the streamlines must be as shown in Fig. 
24(b). 

In such a case, and following the same division of re- 
sponsibility, the propulsion team, concerning themselves 
with the internal flow only, would round and camber the 
intake to avoid any loss of thrust due to internal drag, 





such as might be caused by separation at the leading 
edge. This would naturally leave to the external 
aerodynamicist only the portion of the flow external 
to the stagnation streamline shown in Fig. 24(b). Eq. 
7) shows that the engine aerodynamicist could, by rais- 
ing p2/p, to 2, make V, approach infinity and conse- 
quently reduce the external flow to negligible propor- 
tions, and make the external aerodynamicist super- 
fluous! It is perhaps fortunate that the compressibility 
and imperfections of real fluids, together with the 
greater ease with which density decreases rather than 
increases are made possible by conventional fuels, save 
the external aerodynamicist from extinction and indeed 
give him, in terms of the flow shown in Fig. 24(a), a 
fairly formidable job. 

The difficulty of distinguishing between thrust and 
drag would become even greater if one side of the ram- 
jet were omitted and propulsion were achieved by alter- 
ing the density in the flow around a body. It has been 
shown theoretically by Hicks** and others that a thrust 
may be developed on a cylinder by adding heat to the 
flow around it. While it may not yet be possible to 





propel aircraft by external combustion, it is possible for 
us to avoid the inevitable rigidity of thought which 
stems from trying to departmentalize and compart- 


mentalize our ways of thinking. 
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CONCLUSION 


In conclusion I should like to return to my theme. 
The conception and development of a machine leads to 
more than merely the production of a device. It is 
true that many of us get satisfaction out of making 
things and making things work, but I think there is 
a greater satisfaction than this. It is the satisfaction 
which comes from the growth of our understanding. If 
I have done anything in this lecture to show that 
machines produce ideas just as surely as ideas produce 
machines, I shall have achieved my aim. It is this 
process that makes engineering a profession which offers 
to the best minds the highest of rewards—that of chal- 
lenging intellectual problems which are directly con- 
nected with the most dynamic aspects of life in the 
Twentieth Century, and yet have the enduring value of 


knowledge itself. 
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The Transient Downwash Resulting From 
the Encounter of an Airfoil With a Moving 
Gust Field 


NORMAN P. HOBBS* 
Massachusetts Institute of Technology and AviDyne Research, Inc. 


SUMMARY 


The transient downwash resulting when an airfoil enters a 
moving gust field at an arbitrary rate (not necessarily at the for- 
ward speed of the airfoil) is considered for incompressible non- 
viscous flow For the downwash 
ecurring at downstream positions in the plane of the airfoil is 


The effects of finite span and of vertical displace- 


these conditions, transient 


determined 
ment of the observation point from the plane of the wing are also 
considered for the case in which the airfoil instantaneously enters 
i gust field (or, correspondingly, for a sudden change in angle of 
attack). The theoretical downwash is compared with shock 
tube data obtained by observing the direction of the wake 
emanating from a wire placed downstream from an zirfoil. 
Finally, the downwash for the case of a sudden change in angle 
of attack is used to determine the resulting lift on a two-dimen 
sional airfoil, representing a horizontal tail, in the plane of the 
wake. 
the downwash in dynamically loading the horizontal tail 


This lift, in turn, is used to determine the importance of 


SYMBOLS 


( = airfoil chord 
C(k) = Theodorsen’s unsteady lift function 
Hy, Ay = Hankel functions of the second kind 
Jo, J = Bessel functions of the first kind 
k reduced frequency, we/2U’ 
K = spring constant 
[ = lift 
i = steady-state lift resulting from gust velocity, wp 
m = mass 
= nondimensional time parameter (distance tray 
eled in semichords), 2 l’t/« 
t = time 
l = airfoil velocity 
V = velocity of gust front 
= gust velocity 
= vertical velocity of airfoil relative to fluid 
Wa = magnitude of sinusoidal vertical velocity of ait 
foil 
= downwash velocity from two-dimensional airfoil 
= magnitude of sinusoidal downwash velocity 
= downwash velocity from three-dimensional wing 
= nondimensional distance from the mid-chord of 
the airfoil to the observation point 
y = displacement of single degree of freedom system 


= static displacement of single degree of freedom 
system in response to L 
nondimensional spanwise position of trailing vor 
tex 
= nondimensional airfoil coordinate in chordwis¢ 
direction, 22/¢ 
= chordwise airfoil coordinate 
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s* = nondimensional vertical displacement of obser- 
vation point from plane of wake 
= circulation per unit length on airfoil 


Va = magnitude of sinusoidal circulation per unit 
length on airfoil 

A = nondimensional coordinate in the downstream 
direction 

é = nondimensional coordinate in the downstream 
direction 

o = dummy variable of integration 

y = Kiissner unsteady lift function 

w = frequency 

Q = reduced circulation parameter 


INTRODUCTION 


I DETERMINING the structural response of horizontal 
tail surfaces to gust loadings, transient wing inter- 
ference effects are commonly neglected. Consideration 
of the sense of the vortex sheet shed from a two-di- 
mensional wing in incompressible flow undergoing a 
sudden change in angle of attack indicates that a hort- 
zontal tail located behind the wing will be subjected, at 
early times, to an effective increase in angle of attack 
because of the presence of the wing. The increase will 
last at least until the leading edge of the vortex sheet 
has passed the tail. Under such conditions, it is logical 
to inquire whether this effective increase in angle of 
attack might not lead to important increases in the 
stresses produced in the tail. 

The answer to this question obviously is dependent 
upon the magnitude and time variation of the interfer- 
ence effect, or transient downwash. Accordingly, the 
present paper is largely concerned with the determina- 
tion of the transient downwash produced when a two- 
dimensional airfoil suddently undergoes a change in 
angle of attack in incompressible, inviscid flow. This 
downwash is determined at observation 
points in the plane of the wake; however, the effect of 
vertical displacement of the observation point from the 
Attention is also given 


downstream 


plane of the wake is considered. 
to finite span effects of the wing. 
The transient mentioned 
which arises in the Wagner indicial lift problem,! and 
is accordingly termed the indicial downwash. The 
an airfoil 


downwash above is that 


downwash resulting from the encounter of 
with a sharp-edged gust, the Kiissner problem, is also 
of interest. As may be seen from Fig. |, the Wagner 
problem and the Kiissner problem may both be con- 


sidered to be specialized cases of the same general 
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Fic. 1. Airfoil encountering moving gust field 


problem. Referring to Fig. 1, it is seen that the gust 
boundary itself has a horizontal velocity, |’. When 
I” becomes infinite, the airfoil is instantaneously im- 
mersed in the vertical gust field, or, equivalently, sud- 
denly undergoes a change in angle of attack. Thus, 
the Wagner problem corresponds to an infinite value of 
l’. Similarly, when J is zero, the airfoil flies into the 
gust field at its own velocity, UW’. This, then, is the 
sharp-edged-gust problem considered by Kiissner.* 


It should be mentioned that the lift response to a 
traveling gust has recently been calculated by Miles.’ 
This lift response has also been determined independ- 
ently by Diederich and Drischler! and by the present 
author.’ In addition, in connection with the problem 
of gust simulation in a wind tunnel, Kuethe, Schetzer, 
Garby, and Roensch® have calculated the lift response 
to a moving gust having a given intensity distribution. 
The lift response to a moving sinusoidal gust has been 
determined by Kemp’ in connection with the effect of 
the wakes of stator blades on a rotor blade moving 


through them. 


In the present paper, the transient downwash in the 
plane of the wake is determined for values of |’ between 
positive and negative infinity. It should be noted that 
for | less than minus l’, the moving gust field overtakes 
the airfoil from behind. 


The downwash corresponding to the Wagner problem 
is, as mentioned before, termed the indicia! downwash, 
while the more general downwash corresponding to an 
arbitrary value of I’ is termed the generalized indicial 
downwash. 


The unsteady lift response to the indicial downwash 
of an airfoil representing the horizontal tail is deter- 
mined, and this lift is utilized in determining the im- 
portance of the indicial downwash in producing tail 


stress. 


GENERALIZED INDICIAL DOWNWASH 


The generalized indicial downwash will be developed 
Ly means of a Fourier integral over the known sinusoi- 
dal downwash responses. It is first necessary to spec- 
ify the houndary conditions in suitable forin. Fig. | 
demonstrates the velocities involved—namely, Ll’, the 
forward speed of the airfoil, |’, the velocity of the gust 
boundary; and w», the vertical gust velocity behind 
the gust boundary. The coordinate system used is 
defined in Fig. 2. From Fig. 2, the gust velocity to 


which the airfoil is exposed is given by 
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a. , ta Eas. ( 
7 j0, —(c/2)+ (U+ V)t <2 < (c/2) In Eqs. \” 
Wa = * es : ] 
‘ (—wm, —(c/2)<2< —(c/2)+(U+ Vyt “IL cay 
where w, is the vertical velocity of the airfoil relatiy 
to the fluid, positive up. It is convenient to defiy, 
dimensionless time and distance parameters as . 
g-= BUT ) () 


All dista 
Substitut 


With these definitions, Eq. (1) becomes 


0, —-1+ [1+ (V/U)|s<2<1 


0) 


a = , 
(—uw, —-l<2< —1 + [1 + (V/U))s 
(x* 

Eq. (4+) may be expressed in integral form as 

apy . — Pe 72 2—fi+(V/l 

ay ( 
UW, = - E " dk 

2mt . R 

Eq. (9) | 


where the syinbol ~ indicates that the path of integra 
is 7 present stu 


tion must take an infinitesimal loop below the origin “am 
This is the desired representation of the boundary con , 
ditions. 

Having established the boundary conditions, use o| 
Fourier synthesis to determine the generalized indicial 
downwash requires the sinusoidal downwash responses 
The sinusoidal downwash at observation points in the 





airfoil is given by* 
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sent, respectively, the contributions to the downwash 
made by the bound vortices and by the starting vor- 


tices. In Eq. (6), 
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In Eqs. (6)-(S i is the nondimensional distance from the 
: : ‘ ; midchord of the airfoil to the observa- 
l (x* is the magnitude of the two-dimensional inde: ancl 
; : 10 01 
. downwash at the observation point ' : = 
elatiy * FF z, Xr are nondimensional coordinates in the 
Xs of Sec. . ° 
defi, : 7 . : : downstream direction 
( is the magnitude of the circulation per : ‘ 
: tl , tl tl ‘foil. f k is the reduced frequency parameter 
unit length on the airioil, ft./sec. ‘ . > , 
ot ance HT, IT, are Hankel functions of the second kind 
aD is a circulation parameter used by Reiss , y= ; cae 
a ee ; f indicates the Cauchy principal value ot 
ner,” 1t./ sec. . : 
the integral 
All distances are made dimensionless by dividing by the semichord of the airfoil. 
Substitution of Eqs. (7) and (S) into Eq. (6) and integration of the simpler terms leads to 
x*—-14f' l+2 w,(z) 
™ \ , Ps dz + 
T r* + | \ ; l—sa*-—gsg 
? | » pw ee eh 
rll? + tHe? J_, J; l1—zs YWA-—1x*-—A y 
mais Eq. (9) gives the magnitude of the sinusoidal downwash in response to one frequency component, e* In the 
ne present study, it is desired to determine the downwash response to the step-function gust velocity defined by Eq 
gin 2 ‘ — , ss ; 
j). Combining Eqs. (5) and (9), 
con 
le | is* — 1 a) TOMO! Nee ly 
se of = 
+s i T Vv. + | \, l \ l_—g x* — ¢s 
licial 
ses | °1 ‘ Y jis I L l | > P tk( A+ 1+2 + (VL } 
| the \ \ = dk |dzd\¢ (10 
; rid it; i—2 FW — 1 x* ~All. kIT, ill 
an 
The bracketed hook integral may be simplified in the following way. The integral representation of the Kiissner 
function is given in reference S as 
l 4 1C(R)[Jo(Rk) — tJi(k)| + i(k) fF e* 
y(s) - dk (11 
(6 28 ik 
wre- | Where J) and J; are Bessel functions of the first kind 
ash | and C(k) is the Theodorsen function, 
‘Or- CU IT, (k) : , 
(R ' (12 Differentiating Eq. (14) with respect to the dimen 
2 (k) + iy 2k ; . 
sionless time variable, s, 
Combining Eqs. (11) and (12), and noting that! 
; ; 5 ; d ey ( 7 
Jo(k\I2(k) — Si(k)He®(k) = —(2/mik) (13) ys) = ™ = —~ 2 
ds x. R[IT,° (Rk UT y'*)(R) 
7) | the Kiissner function may be expressed as 
} 1¢ It is seen that the bracketed hook integral in Eq. (10 
| ¥(s } RUT OU 2) dk (14) can be expressed in terims of the slope of the Kiissner 
wT. a | (R) i! (R i‘ ; ; : 
\ | ) function. Accordingly, Eq. (10) becomes 
= . 7 
: | ti ] iv is l+ 2 dz 
T \ ae 1 | \ 1 — 3.x" 


in a een 


Consider the limits of integration on the double integral 


in Eq. (16). y’ is zero for (s + 1—A ((l+s This is valid only as long as the upper limit 1s greater 
+ (V/U)]}) < 0, so the upper limit on the inner than unity, so the upper limit on the outer integral 
integral becomes (s + 1 — ;(1 + 2)/{I (V/U)]}). must be | —1 [1 (V7/U)|s| or unity when this 
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Fic. 3. Indicial downwash. 


limit becomes greater than one. It is convenient to 


define 
A=s— j}(14+2)/[1 + (V/U)}} (17) 
1+1 , 
‘ nN l (1 1 — dA) 
ax Jt u dy (18) 
1 A-— 1 x* — 
Using these definitions, the expression for the downwash 
becomes 
we. ?((V/U), x*, s] | x* — | 
Wo ~ rg a 7 


f 1+[1+(V/U) |s} or 1 tl 1 
ae l1—silLx* -—2 
V 
1A, x*) |e, (1 + r)>° (19) 


Similarly, for rearward immersion of the airfoil, 


Wy ((V/U), x*, s] | x* — 1 
= - x 


f al | , 
(1+ [1+(V/U)|s} or —1 l1—siLx* — 2 
” V 
A, 2) | as, ( + r) <0 (20) 


A=s+}(l —3)/[l1+(V/U)]! (21) 


where 


The indicial downwash, or downwash in the Wagner 
problem, is a special case of either Eq. (19) or Eq. (20). 
When I is infinite, .1 is equal to s and loses its depend- 
ence upon z. The integral in Eq. 19 may then be 
evaluated analytically, and the downwash expressed 
in terms of the integral /. 


Wy'?( 0, x*, s) 


Wo 

I -_ | 
—1+ [I + Ts e")) Qe (22) 

a oe 
It is thus seen that the determination of the indicial 
downwash requires simply the evaluation of the in- 
tegral/. Also, the evaluation of / represents the major 
part of the solution for the generalized indicial down- 

wash, Eqs. (19) and (20). 
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The indicial downwash is plotted in Fig. 3. — It shoy| 
be noted that the time scale has been shifted an amoy 
equal to (x* — 1), so that the horizontal scale is (5 + 
— x*). This shift has been made to emphasize ¢} 
similarity of the downwash at different observati 
points, x*. The time at which the leading edge of th 
vortex sheet reaches the observation point is given | 
s=x* — l,ors+1-—x* = 0. Thus, the origin, 
Fig. 5 represents the time at which the leading edge 
the vortex sheet reaches the observation point. 
this time, the downwash has an infinite value, in q 
upward direction. When the leading edge of the vort 
sheet has just passed the observation point, the dow 
wash has a finite value in the downward direction. 

From Fig. 3 it is apparent that the downwash depen¢: 
largely upon the distance of the observation point fro 
the leading edge of the vortex sheet, particularly a 
short times. At later times, of course, the downwas| 
approaches its steady-state value which depends onh 
upon x*, the distance of the observation point from th 
airfoil. 

The asymptotic behavior of the indicial downwast 
as the leading edge of the vortex sheet approaches the 
observation point may be deduced from Eqs. (18) and 
(22). Defining a dimensionless distance between the 
observation point and the leading edge of the vortex 
sheet, 


| 
| 
| 


Eq. (18) becomes 


mer N+ 1 (x* —A— ©) 
| v : dr (24 


Te, x*) = 
1 A - l x* — i 


The abo 
or transien 
immersion 
ner pré yble 
determine 
sent down 
The result: 
jorward al 
tively. Ft 
the gust fie 
taken fron 

From F: 
tinuous fo! 
ity occurs 
[his repre 
of the vor 
the time ft 
in the gust 
forward 01 
sient dow! 
when the 
passage of 
ton point 
sponding 
slope of t] 
the down 

when the 


observatie 





Sears'' has developed an expansion for the Kiissner 
function for small values of time. The leading term 1s 


‘ / a 
¥V(s)~ V2s/n7, s<l 
The slope of the Kiissner function for extremely small 
values of time thus approaches 
y'(s) ~ 1/9V2s, s<1 (25 


Consideration of Eqs. (24) and (25) demonstrates that 
the main contribution to /(¢, x*) occurs in the region of 
4 =x* — ec. In this region, y’ may be replaced by Eq. 


(25) and the lower limit of the integral in Eq. (24) may 
* 


e 


be replaced by x* — 6, where 6 is small but is much 


larger than «. Accordingly, Eq. (24) becomes 


T(e, x*) ~ 


Introducing this result into Eq. (22), 
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t shoy| 3 — 12 4 —_ 
amc The above remarks apply to the indicial downwash, : | ] 
amour ; hi . ; 
ats or transient downwash corresponding to instantaneous a (14%) me H 
77 ‘ ° es . ‘ . “s ' 
size {h,| immersion of the airfoil in the gust field, as in the Wag- ', 7 =e 4____ | —__—_—}—____.—— 
S1Ze th ‘ 
ervatin| ner problem. Eqs. (19) and (20) were also used to ° (1+¥) = 3 iva 
ation} ” : ny 6 1 
-.. | determine the generalized indicial downwash, or tran- < ra i Tht” | 
e of th nd namie Sa ——— 
sv, sent downwash for immersion rates other than infinity. x (+f =I | Pad U 
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ian taken from behind by the gust field. oe eee dep Kept Mo y=.25 3 
€ - so ree he . ‘ te) | ! 
7 tig. S see ‘ = mwash is discon- 
yee From L ig i, it iS seen that the dow nw ish is liscon NON: DIMENSIONAL TIME PARAMETER, S 
n tinuous for positive immersion rates. The discontinu- 
epends| ity occurs when s = x* — 1 4+ /|2/[{1 + (V/U)] i 
= (Is . . . - . 
‘ ‘his represents the time required for the leading edge wT 
it fron This — os DIRECTION OF 
} of the vortex sheet to reach the observation point plus a POSITIVE V 
rly at 2 ak : zee - oe 
— the time for the airfoil to become completely immersed =F _ a 
Wasi a . = ; : 5 U 
s on! the gust field (regardless of whether the immersion 1s 2 K+ h)--3 _—_—— TT 
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corte: | ope of the downwash is discontinuous. In addition, NON- DIMENSIONAL TIME PARAMETER, s 
the downwash becomes infinite in an upward direction Fic. 4 (Top Generalized indicial downwash for positive 
he » leading edge of . ie chne a eee ; immersion rates (x* = 6) 
- when the leading edge of the Teves sheet reaches the Fic. 5 (Bottom). Generalized indicial downwash for negative 
(<9) | observation point in the case of rearward immersion. immersion rates (x* = 6 
EFFECT OF VERTICAL DISPLACEMENT OF OBSERVATION sidered here. The downwash will be determined by 
POINT using the Biot-Savart law* with the circulation per unit 
(2 ge . : : 
a4 ; length on the airfoil and in the wake. In this study, it 
Heretofore, the downwash considered has been that : : ; : 
: : : ia is convenient to separate the downwash into the com 
ssner| observed at points in the plane of the wake. The 
' : ponents contributed by the wing and the wake. 
mis, downwash at points displaced vertically from the plane ne PRG ; 2 a 
: ee oe SE rhe wing contribution will be considered first. The 
of the wake 1s also of interest. Only the indicial down . i ce 
‘ al circulation per unit length on the airfoil may be found 
wash (downwash in the Wagner problem) will be con Pia ie IE : 
by using Eq. (7) in a Fourier integral. 
nall 
4 ° | 
2w, l é ¢ ‘ti ae dz 
al g = 
25 ; Wigelt YN, sé g 
. > a | - 
hat | tts + 1 ° eta } 
“ : | | | dk | daddy (28 
rl J —1W 1 l-—g2 A-—1? AL RTT, t+ tf 9°? | 
Eq. ~ 
od Utilizing Eq. (15), and evaluating the first integral, 
ich 
| ‘ e+) — , 
fl —€j A+ 1 (s+ 1 d) / nn 
Y AE) . PW ) l . | dr (29) 
l+¢ Ji h- 1 é— dX f 
Using the Biot-Savart law, the airfoil contribution to the downwash at points displaced vertically from the plane of 
the wake becomes 
a : I — § ‘igisdia I + 1 p(s l1—A x* —¢ a 
(} (' : } Y 4 as v 2 = ant ~ dé (30) 
wy J, wT J, i+et J} d l é-— dX 1 (x* — £)? + 2*? 


where s* is the vertical displacement of the observation 
point from the plane of the wake, made dimensionless 
by dividing by the semichord of the airfoil. 

The wake contribution to the downwash may be de- 
termined by using the fact that the wake circulation per 


unit length is proportional to the slope of the Kiissner 


function. 


- d) 2rwmy'(s + 1 — XA) (3 
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os | | | | values of x* due to the increased importance of th: sonal dov 
12 2*-0 z*-2 bound vortex. tribution. 
° g | the variat 
E ks EFFECT OF FINITE ASPECT RATIO 
a 4 i i J } } Bis ; P 7 ote CoMPARIS 
x Finite aspect ratio of the wing affects the downwag 
7} Z=10 “a ja ; ' 
$ a at the tail in two ways. The first is that the distrjhy 
FS tion of the starting vortex system is different for The ex] 
2 infinite and a finite aspect ratio wing due to the addej | wash 15 1 
ae a is 20 25 presence of trailing vortices. However, for moderat. | attammmen 
aspect ratios, this effect s of primary in. | ficult. T 
NON- DIMENSIONAL TIME PARAMETER, s aspect ratios, this effect should not be of primary in ( | 
portance, since near zero time the additional trailing} img a SteP 
Fic. 6. Effect of vertical displacement of observation point : ; 7 hock 
“ . ae ee , Ly cycte , > OIC 7 mo F rn n. shock 1 
from plane of wake on indicial downwash (x* = 10) vortex system due to the gust is short in length and , > 
therefore have only a small influence on the flow. ;] studying 
, ' later times, this influence becomes more important, byt | vertlicatio 
Using Eq. (31) and the Biot-Savart law, the wake con- ae , See : ‘ ; : 
te Aisi uid the Biot-Sava , ' . it is the starting vortex distribution near zero time which | attemptes 
ribution becomes ek ae ie ‘ ; ty niet 
——— is Of primary importance in determining the downwash An atrt 
The second effect is the loss of the starting vorte the test s 
apy (2 es+] : i , me ; 
(7 ) <a v(s +1—\) xX system between the wing tips and infinity and the pres the tube 1 
| ‘ < - : e 
We J ence of the trailing vortex system. The starting and] object 0! 
cP — : trailing vortex systems assumed are shown in Fig, 7.| of the wa 
a — Gr (32) se ee, i oe : letermine 
(vn* — vr)? + 2*? As may be seen from Fig. 7, the trailing vortices are} “eter 
assumed to be grouped at the dimensionless spanwise | foil. 51m 
sf : ‘ : ms Stites. aah ye: niin: sree cco by dividine| behind th 
Phe results of evaluating Eqs. (30) and (32) are shown distance, y*. y* is made dimensionless by dividing} ™ 
by the semichord. The dotted lines in Fig. 7 then rep cools the 


in Fig. 6 for x* = 10. The greatest influence of the 
vertical displacement of the observation point is seen resent the vortex system which must be subtractel}| ¢Ve - 
to occur in the neighborhood of s = «* — 1. As pointed from the two-dimensional result to obtain the finite | result is 1 

; the wire | 


out before, this represents the time at which the leading 





aspect ratio downwash from the starting vortex sys 


edge of the vortex sheet reaches the observation point tem. The observation point is taken at mid-span later. | Mte-spar 
ally, so that the influence of the dotted vortex system is} the ebser 


(or, in the present case, reaches a point directly above 
. : . , ‘trical fror » two sides stream all 

or below the observation point). Physically, the reason symmetrical from the two sides. ’ sl 
é - he 7110 

be negligi 





that the main effect of vertical displacement of the ob- The finite aspect ratio correction for x* = 7 is show! 
servation point occurs near s = x* — | is roughly as in Fig. S for various values of y*. This correction [ ee 
follows. From Eqs. (31) and (25), it is seen that at the must be added to the two-dimensional results to obtai _ : ” 
leading edge of the vortex sheet y jumps’ discontinu- the finite aspect ratio downwash. The downwash at The ag 
ously from zero to infinity and then decreases rather v* = 7 is shown in Fig. 9 for y* = © (infinite aspect ented : 
rapidly. The vortex sheet is thus roughly similar to a ratio) and for y* = 6 (aspect ratio of about eight ontonng 
concentrated vortex at the leading edge of the sheet. For these conditions, it is seen that finite aspect rati cca 
vortex sh 


Consider, then, a concentrated vortex in a plane dis- affects mainly the steady-state downwash and is ol 


placed vertically from the observation point a dis- 
tance 3*. When the horizontal distance from the 


comparatively little importance near zero time. The resul 


It should be noted that Fig. S can be used to deter 
mine a finite aspect ratio correction to the two-dimen 


values bi 
: : ; shock fro 

vortex to the observation point is large compared 
reached 1 


with s*, the vortex produces essentially the same down- . 
conmlparist 


wash as if s* were zero, resulting in little change 
the shock 





in the downwash due to vertical displacement of the ob- | 






































ae , ; Coho. 
servation point from the plane of the wake. When the y ry a 
; : age | a | 
horizontal distance is small compared with 2*, however, | a 
the vortex can induce little vertical velocity at the ob- iW “slo ° 
servation point, but for the same horizontal distance, if t | | ‘| 
i. ; : : z 
2” were zero, a large vertical velocity would be induced. TRAILING os 
; P | - VORTEX 
Accordingly, the magnitude of the downwash should be Le iw 
greatly diminished when the horizontal distance from —+—--—-— — —i----- 8 
° . . ~ - 
the leading edge of the vortex sheet to the observation | x Sars ° 
ane : ; , } | vi i 
point is small compared with s*, and should be virtually Stl < 
unaffected when the horizontal distance is large com- OBSERVATION z 
? i : eae POIN 
pared with s*. These effects are shown in Fig. 6. a ac = SS = 
2 
s* > . - * = 
Although Fig. 6 applies only for x* = 10, the same a SS Srenn 8 
trends depicted therein would hold for other values of —t---—_ + ——--——— 
x* as well, the main difference being that the steady- 
; Fic. 8. | 


state downwash would be more affected for smaller Fic. 7. Vortex system for finite wing. 
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sonal downwash for any assumed trailing vortex dis- 


tribution. In addition, it may also be used to indicate 


the variation of downwash along the tail span. 


INDICIAL DOWNWASH WITH SHOCK 


TUBE RESULTS 


COMPARISON OF 


The experimental verification of the indicial down- 
wash is not a simple problem, because the experimental 
attainment of the necessary conditions is extremely dif 
The shock tube offers the possibility of simulat 
The MIT-WADC S- by ?4- 
was designed with the idea of 
Experimental 


ficult. 
ing a step-function gust. 
in. shock tube!” at M.I.T. 


flow phenomena. 


studying unsteady 
verification of the indicial downwash was accordingly 
ittempted in the MIT-WADC shock tube. 

An airfoil was placed in the shock tube upstream of 
the test section, and a thin wire was stretched across 


The 


ibject of the experiment was to measure the direction 


the tube between the windows at the test section. 


of the wake emanating from the wire, and thereby to 
determine the downwash at the wire caused by the air 
joil. Since the wire is at room temperature, and the air 
behind the shock front is heated by the shock, the wire 
At the same time, how 
ever, due to viscous effects, the air is heated. The net 
result is that the wake is heated, making the wake of 


cools the air passing over it. 


the wire visible on interferograms as distortions in the 
The temperature of the wake at 


finite-spaced fringes. 
above the free 


the observation point is about 3.5°F. 
stream air temperature, so that buoyancy effects should 
be negligible. 

A 10 per cent thick symmetrical double-wedge airfoil 
with a chord of 3.97 in. was used in the experiment. 
The airfoil was placed at an angle of attack of approxi 
about 10 semichord lengths upstream from 


mately 5 
Runs were made to determine the timewise 


the wire. 


position of the starting vortex (or leading edge of the 


vortex sheet), which was visible in the interferograms. 
The resulting positions were compared with calculated 


values based on the free-stream velocity behind the 


shock front using the time at which the shock front 
reached the airfoil trailing edge as a reference. 
comparison demonstrated that the starting vortex in 


the shock tube was approximately 1/4 of a semichord 
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Fic. 8. Finite span corrections to indicial downwash (x* = 7 
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length of flow late in arriving at downstream positions. 
The vortex appeared to be moving at free-stream veloc 
ity when it passed the test section, approximately 9 
semichord lengths behind the trailing edge of the air 
foil. The duplication of the position of the vortex in 
successive runs was excellent, the average streamwise 
variation being about 0.03 in. and the average vertical 
variation being about 0.02 in. 

On the basis of these experiments, which determined 
both the vertical and longitudinal positions of the vor 
tex, the airfoil was placed vertically in positions such 
that values of 2* of +0.5 and +0.2 could be obtained 
The Mach Number of the flow behind the shock front 
was 0.4, and the Reynolds Number was approximately 
One of the resulting inter 

In Fig. 10, the starting 


650,000 based on the chord. 

ferograms is shown in Fig. 10. 
vortex and the wake from the wire are both 
Although no reflected shocks (reflections of the main 
shock from the airfoil) or trailing shocks (weak shocks 


which follow the main shock) appear in this particular 


visible. 


interferogram, these shocks represent one of the main 
obstacles in the present experiments. The reflected 
shocks, in particular, can influence the vertical velocity 
at the wire quite strongly compared with the vertical 
velocity caused by the The first reflected 
shock in the present experiment arrived 
The trailing shocks arrive at random 


airfoil. 
at approxi 
mately s = 6. 
times after the arrival of the main shock. 





wire 


wake 


Fic. 10. Interferogram showing vortex and from 


s = 8.75, s* = O5 

















































738 JOURNAL OF THE AERONAUTICAL 
8 “——_ a | EXPERIMENTAL 
| POINTS 
Re = 650,000 | 
| x g¥e+5 
TT apo. TO 
a | | ' 
Mzi0 | | | | | 
2 me ; See ~ 
o | THEORETICAL 
Fd | INDICIAL 
a | DOWNWASH | 
= 2 <<<} 
= " | 
77) } | | 
< 
oO = J 
: — a a = 
° | 
a | | 
-2} —}+——+} 
— 2 4 6 8 16 
NON-DIMENSIONAL TIME PARAMETER, s 
6 : 
M+0.4 
Re * 650,000 
- 7 -2—-T HEORETICAL - 
“3,2 INDICIAL 
ai? DOWNWASH 
7s 7 ’ | | | 
‘a 
« Ln | 
= © T 
a | | 
<= | | 
= | * 
z 
zx -2 — - 
o 
a | | 
a 2 4 6 E 12 14 16 
NON-DIMENSIONAL TIME PARAMETER, s 
Fic. 11 (Zop). Comparison of indicial downwash with shock 


tube data (x* = 10, 2* = +5) 


Fic. 12 (Bottom). Comparison of indicial downwash with shock 


tube data (x* = 10, s* 


= ().2) 


Aside from the problem of reflected and _ trailing 
shocks, the accuracy of reduction of the interferograms 
This corresponds to an error 
+0.1. An- 


is only about +0.5°. 


in the downwash ratio, w,,‘?)/wy, of about 


other error is introduced by the fact that the direction’ 


of flow produced by the main shock can apparently 
vary by as much as +0.5 

Since the shock-tube conditions do not conform to 
those assumed in determining any of the generalized 
indicial downwash functions, the question arises as to 
which of the downwash functions can most reasonably 
be compared with the experimental results. It is ap- 
parent that the airfoil in the shock tube has a finite 
immersion rate given by the speed of the shock wave. 
rhis consideration suggests the use of the generalized 
indicial downwash function for the appropriate immer- 
rate. However, in the theoretical analysis for 
finite immersion rate, it is assumed that the entire air- 


sion 


foil, and notably the trailing edge, is instantaneously 
influenced by a disturbance at any point on the airfoil 
flow). Thus, the 
begin to form when the gust front reaches the leading 


incompressible starting vortices 
a 


edge of the airfoil. This, of course, is far from the case 
in the shock tube. Since the shock itself is supersonic 
relative to the flow ahead of the shock, the trailing edge 
is unaware of the disturbance until the shock reaches 
the trailing edge. At that time, the lateral separation 
between the shock passing over the lower surface of the 
airfoil and that over the top of the airfoil causes flow 


upward past the trailing edge, which is the beginning of 
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the formation of the vortex sheet. Since the speed of 
sound behind the shock is greater than the shock veloc. 
ity, the consequences of the disturbance over the entire 
airfoil begin to be felt at the trailing edge when the 
shock reaches the trailing edge. This means that the 
noncirculatory flow past the trailing edge corresponding 
to the Wagner problem is established over a very short 
period of time after the shock reaches the trailing edge. 
It is believed, therefore, that the indicial downwash 
function, which corresponds to the Wagner problem, 
beginning at the time when the shock reaches the trail. 
ing edge, should correspond reasonably with the results 
of the shock-tube experiments. 

The comparison of the indicial downwash with the 
shock-tube results is shown in Figs. 11 and 12 for z* = 
+0.5 and +0.2, respectively. The dimensionless time 
parameter, s, was calculated from the free-stream veloc- 
ity and the airfoil chord, using the time at which the 
shock reaches the trailing edge as zero time. The lag 
of 1/4 of a semichord length of flow in the position of the 
starting vortex mentioned previously was not taken into 
account in plotting the data. This means that if it 
were desired to compare the experimental data with 
the theory on the basis of distance of the leading edge 
of the vortex sheet from the observation point, the time, 
s, for each of the experimental points should be reduced 
by 1/4. 
or minus 1/S of a semichord length of flow. 


The time resolution is believed to be about plus 


It is apparent from Figs. 11 and 12 that the experi- 
iments neither confirm nor deny the theoretical results, 
the 
the theoretical and experimental results. 


although same general trends are evident in 
Perhaps the 
greatest value of the experiments lies in pointing out 
ways in which the experimental techniques could be 
improved. First, it would be wise to use higher angles 
of attack of the airfoil in order to produce larger and 
more easily measurable angles of flow of the wake from 
In doing this, an airfoil with a rounded lead 


an airfoil of 


the wire. 
ing edge should be employed. Second, 
smaller chord could be used, delaying the arrival of re 
flected shocks at the observation point till a later value 
of dimensionless time, s, for the same value of dimen 
sionless position of the wire. Using a chord of 2 in. 
(approximately half that used in the present expert 
ment), the first reflected shock would arrive at about 
= 10 instead of s = 6. In using a smaller chord, car¢ 
would have to be taken with the Reynolds Number, 
since Reynolds Number effects may be different in 
unsteady flow than in steady flow experiments. Fi 
nally, attempts should be made to improve the basic 
flow by eliminating the trailing shock waves insofar 
as possible and by making the direction of flow mor¢ 


repeatable. 


HORIZONTAL TAIL LirT RESPONSE TO 
INDICIAL DOWNWASH 


In order to determine the structural response of the 
horizontal tail to the downwash, it is first necessary to 


determine the lift response of the tail. Attention 1s 
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restricted here to the lift response of a two-dimensional 
airfoil in the plane of the wake to the indicial downwash. 
The effects of vertical displacement of the tail from the 
plane of the wake, of finite wing span, and of finite im- 
mersion rate of the wing on tail lift are not considered. 
[he qualitative effects of these parameters are appar- 
ent, however, from consideration of Figs. 4-9. 

From Fig. 3, it is apparent that, over a small range of 
downstream positions, the downwash is essentially a 
function of the distance of the leading edge of the vor- 
tex sheet from the observation point. Thus, essentially 
the same downwash sweeps over the tail at a rate equal 
to the forward speed of the tail. The Kiissner func- 
tion may thus be utilized to express the lift response of 
the tail. Use of the Kiissner function is not strictly 
correct, especially near zero time when use of the Kiiss 
ner function inherently assumes that the disturbance 
exists only at the leading edge of the tail. This 
course is not the case according to the incompressible 


of 


downwash predicted herein; however, the downwash, 
and hence the lift, is small at small times, so that this 
discrepancy is unimportant. 

In determining the tail lift response to the indicial 
downwash, it is convenient to base the nondimensional 
time parameter, s, on the horizontal tail chord, c¢, 


that is, 


i= ZUG = (6/6). (33) 


where s, is the nondimensional time parameter based 
on wing chord, previously called simply s. 

The lift response to the indicial downwash is shown 
in Fig. 13 for a chord ratio, ¢,/c,;, of two and for the 
downstream position of the tail based on the leading 
edge, x*,z, equal to seven. The quasi-steady lift 
is also shown for comparison. 
quasi-steady lifts are made dimensionless by dividing 
by L.., the steady-state lift resulting from the gust 
The infinite lift which occurs when the 


The unsteady and 


velocity, wy. 
leading edge of the vortex sheet reaches the leading 
edge of the tail is of apparent mass origin. 

THE HORIZONTAL TAIL 


STRUCTURAL RESPONSE OF 


In considering the structural response of the hori 
zontal tail to the indicital downwash, the tail will be 
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Tail lift due to indicial downwash (x* = 7, Cw/cr = 2) 
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Fic. 14. Structural response of horizontal tail to gust and 
indicial downwash (x* = 7, k = 0.757) 


represented by a single degree of freedom system having 
a frequency, w, which is representative of the first 
natural frequency of a horizontal tail. The frequencies 
considered also include fuselage first vertical bending 
frequencies. The equation of motion for the response 
of an undamped single degree of freedom system repre 


senting the horizontal tail is 


m(d*y/dt?) + Ky = Lit) (34) 


where y is the deflection, m is the mass, and A is the 
spring constant. Since the lift is expressed in terms of 
the dimensionless time parameter, s,;, it is more conven 


ient to rewrite Eq. (34) in terms of s,. 


m(2U/c,)*(d*y/d?) + Ky = L(s,) (35) 
The solution of Eq. (35) is 
Ky =k | L(o,) sin k(s, — o,)do, (36) 
where 
i — oe? ae 
k = VK/m (c,/2U) = we,/2l (37) 


w is the frequency of the single degree of freedom system 
and k is the reduced frequency used previously. The 
static response to L«, which will be designated yo, 1s 
given by 


3S) 


Ky = L 
Dividing Eq. (36) by Eq. (38), 


y/yo =k 39) 


| [L(o,)/L@] sin k(s; — a,)do, | 


The gust lift and the downwash lift were used sepa 
rately as excitations in Eq. (39). A characteristic 
response is shown in Fig. 14. From Fig. 14, it is pos- 
sible to determine the ratio of the peak response to the 
gust and downwash combined to the peak response to 
the gust alone. This ratio will be called the response 
ratio. It should also be noted from Fig. 14 that the 
peak response to the gust and downwash combined 
occurs at the fifth response peak for the particular 
conditions applicable to Fig. 14. The response ratio 1s 
plotted in Fig. 15 for x* = 7. 
response peak at which the peak response to the gust 
Fig. 15 


The number of the 


and downwash combined occurs 1s indicated. 
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Fic. 15 Respouse ratio (x* = 7 


indicates that the importance of the indicial downwash 
can be considerable in producing deflections of, and 
therefore stresses in, the horizontal tail. 


assumed conditions and the conditions which will ac- 
tually exist before conclusions can be drawn as to the 
actual structural importance of the indicial downwash. 
The first, and probably the most important, factor is 
that aerodynamic and structural damping have been 
neglected in the present analysis. Actually, the re- 
sponse to the gust will decay considerably in a few cycles 
due to damping. This means that for cases in which 
the peak response to the combined gust and downwash 
occurs after several oscillations, the base from which 
the response to the downwash operates (that is, the 
response to the gust) is considerably smaller than that 
assumed. It is seen from Fig. 15 that the highest re- 
sponse ratios occur when the downwash is phased such 
that the peak response occurs after several oscillations. 
The inclusion of damping would normally considerably 
When the 


peak occurs on the first or second oscillation, the re- 


diminish the response ratios in this area. 


sponse ratio would be diminished somewhat, but not 
to a very large extent. 

A second factor to be considered is motion of the air- 
craft. The airplane tends to ride with the gust, and 
the resulting vertical translation decreases the total 
lift on the tail. 
mally tends to decrease the lift. 
dependent, and will be most important for light air 


In addition, the pitching motion nor- 
This effect is time 
planes. Assuming constant mass and inertia of the 
airplane, the response ratio will be decreased most by 
aircraft motion for large values of tail distance, x* 
since, for large values of x*, the airplane has had more 
time to accumulate velocity in the direction of the gust 
before the large upwash reaches the tail. 

An allied effect is the combination of motion of the 
vortex sheet with the gust and the motion of the air- 
plane. These motions mean that the tail cannot al- 
ways be in the plane of the wake. However, in con- 
sidering the maximum possible structural importance 
of the downwash, if it is assumed that the tail is located 
in such a way that the portion of the vortex system 
responsible for the large upwash passes very near the 


RONAUTICAL 


‘stresses since use of 
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horizontal tail, this effect will be of little consequence 
For other vertical displacements of the tail from the 
plane of the wake, the qualitative effect may be as. 
sessed from Fig. 6. 

The influence of finite wing span must also be con 
sidered. From Fig. 8, it is possible to obtain an order 
of magnitude of this effect. The response to the finit, 
span correction in the period of interest for moderat; 
aspect ratios will be on the order of 5 per cent of th 
(This js 


obviously only a rough order of magnitude, since it does 


response to the gust, in a negative direction. 


not involve variation with aspect ratio or downstrean 
position of the observation point.) A decrease of 
per cent (of the gust response) in the response to th 
downwash will decrease the response ratio at most by 
about The influence of finite wing spa 
is thus seen not to be of primary importance. 


The effect of finite tail span on the response rati 


) per cent. 


should be of even less importance, since it will affect 
the tail lift due to the gust and due to the downwash 
similarly. 

One final point on the influence of finite span of the 
wing should be made. The analysis presented hereir 
demonstrates conclusively that use of the [1 — (de/da 
steady-state downwash correction would lead to ver 
unconservative results in predictions of tail loads and 
this factor would reduce the r 
sponse to the gust alone by 30 or 40 per cent. As 
minimum, the || — (de/da)| factor should be omitted 
in sharp-edged-gust studies of tail response. 

Finally, the effect of immersion rates other than in 
finity must be considered. For negative immersioi 
rates, the tail will be subjected to the gust before th 
so that there will be a lag between the gust re 
sponse and the response to the downwash. This lag 
will aiford time for the gust response to damp out to 
fairly large extent before the tail is subjected to the im 


wing 


Ss? 


portant upwash. The case of rearward immersion 1s 
accordingly probably less critical than is the case | 
infinite immersion rate. 

For forward immersion of the wing, it is possible t 
have immersion rates such that the peak response t 
the downwash will phase in with the first or second peah 
For such the effect 0 


of the gust response. cases, 


damping will be considerably lessened. The curves ¢ 
response ratio in Fig. 15 should give reasonable est 
mates of the possible importance of the downwash fo 
forward immersion rates. The case of forward immer 
sion rate with proper phasing between the gust respons 
and the downwash response may well be the most criti 


eal. 


CONCLUDING REMARKS 


In summary, it is estimated that the peak respons 


ratio for infinite immersion rate, including the eifects 
of damping, will be on the order of 1.4 and will occu 


for the lower values of reduced frequency. For nega 


tive immersion rates, the peak value of the response 


Continued on page 754 
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Cross-Coupling Dynamics and the Problems of 
Automatic Control in Rapid Rolls 


J. D. WELCH* anno R. E. WILSON* 
General Electric Company 


SUMMARY 


Some of the fundamental properties of inertial cross coupling 
e reviewed, with specia) emphasis on the inertia coupling terms 
vhich may cause apparently similar aireraft to exhibit different 

racteristics. Typical time histories of all the components of 
ircraft accelerations are presented to give a detailed insight into 
the dynamies resulting from inertial cross coupling 


Phe problems of automatic control in minimizing objectionable 


haracteristics during rapid rolls are discussed. A_ siniplified 
iethod for analyzing the contro] problem is presented as well as 

discussion of the practical problems such as measurement 
techniques and reliability 

SYMBOLS 

] = roll moment of inertia 

| = pitch moment of inertia 

| vaw moment of inertia 


cross product of inertia 
= roll acceleration due to unit sideslip angle 
= roll acceleration due to unit rudder deflection 


4413 


L roll acceleration due to unit spoiler deflection 


= roll acceleration due to unit roll rate 


roll acceleration due to unit yaw rate 


. pitch acceleration due to unit angle of attack 
Va pitch acceleration due to unit rate of change, angle of 
attack 
V5, = pitch acceleration due to unit tail deflection 
V, = pitch acceleration due to unit pitch rate 
Vz = yaw acceleration due to unit angle of sideslip 
V5, = yaw acceleration due to unit rudder displacement 
= yaw acceleration due to unit spoil r deflection 
V, = yaw acceleration due to unit roll rate 
VY. = yaw acceleration due to unit yaw rate 
p = roll velocity 
pitch rate 
Vaw rate 
Laplace operator 
I = forward velocity 
Vg = angle of sideslip rate due to unit angle of sideslip 
} angle of sideslip rate due to unit angle of rudder ce 
flection 
Ze = angle of attack rate due to unit angle of attack 
Zsr = angle of attack rate due to unit angle of horizontal tail 
deflection 
a ingle of attack 
a = initial angle of attack and also angle between stability 
axes and principal axes 
3 angle of sideslip 
6 rudder deflection 
6 = aileron deflection 


Presented at the Stability and Control Problems During Large 


Maneuvers Session, Twenty-Fifth Annual Meeting, IAS, New 
York, January 28-31, 1957 
* Aircraft Control Systems Section, Light Military Electron 


Equipment Department 


spoiler deflection 


6; horizontal tail deflection 

o roll angle 

/., = moment of inertia of engine rotor 
p engine speed 


quantity in stability axes system 


CROSS-COUPLING DYNAMICS IN 


RAPID ROLLS 


PART (I) 


A) INTRODUCTION 


Mr RN HIGH-PERFORMANCE aircraft may have 
complex cross-coupled response problems during 
rapid rolling maneuvers. This is due primarily to these 
aircraft having their mass concentrated in elongated 
fuselages and, in addition, being capable of high roll 
rates. The dynamics of the aircraft motion including 
cross coupling have been the subject of considerable 
investigation during recent years, with most of the 
studies emphasizing the dynamics of a specific aircraft. 
In this paper we shall attempt to generalize on the 
dynamics as they would apply to various aircraft and 
clarify some of the apparent inconsistencies that appear 
in the dynamics of various aircraft. 

The prime objective of the first part of this paper is 
to present a physical interpretation of the major factors 
involved in the cross-coupling phenomena. This di 
cussion of the basic dynamics of the problem will serve 
as a background for Part (II), which will discuss the 
important topic of automatic control to eliminate the 
objectionable characteristics of cross coupling. Our 
analysis considers not only the theoretical aspects of 
the control problem but considers many of the prac 
tical ramifications as well. 

CRross-COUPLING 


REVIEW OF THE (GENERAL 


PROBLEM 


(B) A 


The cross-coupling problem during fast rolling maneu- 
vers can be roughly divided into divergent or near 
divergent response and convergent response but with 
undesirable characteristics. 

Most of the effort in this country appears to have 
been directed toward the first of these problems. 
Phillips, in his well known report,? has shown that an 
“critical” roll velocity may 
experience a divergence in one of its modes. The 
sustained 


aircraft rolling above its 


Phillip’s criterion predicts divergence or 


t Phillips, W. H., tudinal and 
Directional Stability, NACA TN 1627, 


Effect of Steady Rolling on Long 
1948 
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oscillation during a continuous roll if the roll rate ex- 
ceeds the lower of either the yaw or pitch nonrolling 
natural frequency. The situation is analogous to the 
classical stability problem of a rotating shaft. The 
range of roll rates between the two nonrolling natural 
frequencies is the “‘critical’’ range of roll velocity where 
divergence may occur. Several investigators have 
pointed out this divergence by techniques such as root 
locus plots. Furthermore, flight tests have substan- 
tiated the response predicted by the Phillip’s criterion. 

Quite apart from this divergent response, the prob- 
lem of response for convergent motion may be most 
important for many aircraft during practical rolling 
maneuvers. In fact, some aircraft may not be ca- 
pable of rolling fast enough (relative to their nonrolling 
natural frequencies) to get into the critical divergent 
region. Yet these same aircraft may experience severe, 
perhaps even hazardous, response characteristics at 
those rates below the critical value. 

Whereas much of our presentation will be applicable 
to both of these problem areas, we shall place emphasis 
on the convergent response which we believe has been 
insufficiently emphasized in the past. As a result, all 
of our examples of aircraft response in this paper will 
demonstrate nondivergent response characteristics. 


(C) THe Various Cross-COUPLED COMPONENTS OF 
THE AIRCRAFT FORCES AND MOMENTS 


General 

Two types of cross coupling occur in the dynamic 
response of the aircraft, inertial and aerodynamic. The 
inertial cross coupling is the cause of disturbances 
which result in aerodynamic cross coupling and is 
are taken to 
and 


generally more significant. If means 
minimize objectionable excursions in 
angle of attack, the aerodynamic coupling may be 
In this report most of our attention will 
In the next 


sideslip 


negligible. 
be focused on the inertial coupling terms. 
section, however, certain manifestations of the aero- 
dynamic coupling will be demonstrated. 


The Choice of Axes Systems 

One difficulty in discussing the coupling problem is 
the fact that some investigators have used the stability 
axes coordinate system* while others have used the 
principal inertial axes system. 

The principal axes system is preferable for the con- 
trols engineer because, for all flight conditions, it re- 
mains in fixed orientation with respect to his sensing 
instruments and because it avoids the added complica- 
tion of the product of inertia terms in the analysis. 
We believe, however, that certain physical descriptions 
of the various factors involved in the cross coupling 
can be accomplished best by referring to the equations 
of motion in the stability axes system. 

Before we begin a physical description of the cross 


* Stability axes referred to in the paper are body axes aligned 


initially with the velocity vector. 
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coupling, we would like to clarify the relationship of 
the two systems. The equations describing the fiy, 


degrees of treedom of the aircraft in stability axes are 


Bs = [Lys — Tes)/TrslQsts + 
(Teze/Lns) (%s + Psds) + UM ry (acm 
G = [Le, — Trg) /Tys Pst — (resus) Pes ts + 
(Lzes/Tye) (727 — Ps®) + TM yy caer 
t= [Ur, — Lys)/Tus|PsGs + (Cres /Isz) Pes lds + 


(Lr25/L5) (Ps — Gs%s) + TM, 
B, = asp, — r, + IF,, 
— Bp, + IF, 


(aero and gravity 


as = qs aero and gravity 


where the subscript s refers to quantities measured in 


the stability axes system. 
In the principal axes system these equations become 





p= (WC, LS/1 Ne A SMe tace 
J [7, — 12)/Ty|pr — (Lrehe/Ty)r + 2M, (acre 
* = [U, — 1,)/T.)pq + (Lecbe/T2)9 + 2Msz acto 


B i ap waste 3 =p OP stew end gravity) 
sate, eae Bp ar is tents wel gravity) 


~ 
= 


The only way in which the equations differ explicitly is 
the inclusions of the product of inertia terms in the 
stability axes system. Not explicitly shown is a shift 
in the angle of attack zero reference by the angle ay be- 
tween the principal axis and the flight path axis when 
we go from one system to the other. A simple coor- 
dinate transformation shows that the only other major 
change is in the yaw degree of freedom. The yaw rates 
differ in the two systems approximately by the incre- 
ment pay and the yaw accelerations differ approxi- 
mately by the increment pay. A quick check of prac- 
tical values of a, roll rate and roll acceleration shows 
that, for a given maneuver, yaw rates and yaw ac- 
the vastly dif 





celerations in two systems may be 


ferent. 


Inertial Coupling 


The inertial coupling terms appearing explicitly in 


the force equations are the terms ap and Bp. These 


are best explained by considering the principal axes 
system, momentarily neglecting all other forces in the 
side and lift degrees of freedom and assuming no yaw 
or pitch moments. Fig. | gives a simplified concept of | 


the mechanics involved. The equations become | 


B ap, @ —Bp 


If these equations are solved for a and 8 assuming a 
constant roll rate, a sinusoidal component of angle of 
attack and angle of sideslip results with a radian fre- 
quency equal to the roll rate. In the typical case with 
an initial trim angle of attack ap, and no initial sideslip, 
the response would begin with a decrease in angle of 


attack and an increase in sideslip angle. It would con- 
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AIRCRAFT IN INITIAL TRIM AIRCRAFT AFTER RAPID 
90 DEGREE ROTATION 


Simplified dynamies of lift and sideslip considering only 
the inertia coupling terms in the force equations 


tinue as a cyclic interchange of relative wind angle in 
the lift and side directions. If there were no other 
forces or moments involved, the aircraft would tend to 
roll about its minimal inertial axis with no yawing or 
pitching motion. The lift force and side force would 
vary sinusoidally, and be relatively displaced in time 
by 90”. 

Prior to discussing the inertia coupling terms in the 
moment equations let us consider the major factor which 
causes the initial disturbing accelerations namely, 
the inclination of the principal inertial axis relative to 
the flight path or stability axes. For a wide class of 
practical flight conditions and maneuvers, the initial 
angle of attack relative to the zero lift line for typical 
fighters may range from +12° to —4°.* Typical 
fighter aircraft may have an inclination of the principal 
inertial axis ranging from 0° to 4° below the zero lift 
axis. Considering these facts, we find that the prin- 
cipal inertial axes may be oriented as much as 8° below 
the stability v-axis, giving maximum negative values 
of a9 to as much as 12° above the stability .-axis, giving 
maximum positive values of ap. These relationships 
are shown in Fig. 2. 

We will now discuss the moment equations. 
make some assumptions concerning realistic inertias 


Let us 


so that we can determine the effect of the inertia mo- 
ment terms. Assume that /, iP 10/,. These are 
typical relationships for a modern fighter. 
that /,. (J, — /,) sin ay, the pitch and yaw moment 


equations in the stability axes system become 


Recalling 


= 0.9p,7, — 0.9a0(7,* — £7) — 


(ft leet Migs Gor 


gd 
Us 


-, = O.9p = O.9ao( ps — qrs) ob 
ep. I )q + J/ “TO 


Where sin ay & a for the angles in which we are inter- 
ested. The inertial coupling terms in the roll equation 


are very small compared to the aerodynamic rolling 


* Angle of attack for the landing condition may lie outside the 
This is of little interest here since roll rates will 
Also just what constitutes the practical 


region shown. 
be low during landing 
limit of negative load factor maneuvers may admit of consider 


able argumentation 


~I 
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moments and may be neglected in most cases. Also, 
the apr,” component in the pitch acceleration response, 
the a (—gr,) component of the yaw acceleration re- 
sponse, and the two coupling terms due to the engine 
will be small compared to the other components and 
can be neglected in our first analysis. 

The pitch and yaw moment equations then become 


Gs = O.9pP.rs + OYaop.? + Mys (acro 
és = 0.9p.q — O.9a0p. + Mex (oro 


We have called the first term on the right-hand side 
of each moment equation the gyroscopic term, becau:e 
the action of these components is similar to that which 
occurs during the precession of a conventional gyro 
scope —1.e., a rigid body having angular velocities about 
two mutually perpendicular axes will experience an 
acceleration or torque about a third axis perpendicular 
to the other two. Thus, for example, if we simul 
taneously roll and yaw the aircraft, it will experience a 
pitch acceleration equal to [(/, — /,) J,|pr. 

The other inertial terms appearing in the moment 
equations are due to the cross product of inertia which 
is consequent upon our choice of axes system. The 
most significant effect here is the ‘centrifugal’ pitch 
ing moment due to high roll rates. If the principal 
inertial .-axis is below the stability axis, this component 
will pitch the nose down and if above, it will tend to 
pitch the nose up. If we were describing the dynamics 
in the principal axis, these ‘‘centrifugal” pitching effects 
Actually the effect would 
This 


would tend to be obscured. 
be described by a modified gyroscopic term. 
becomes clear when we recall that the yaw rate differs 
in the two systems by the term ap. 

The other significant component due to cross product 
of inertia is the yaw acceleration described by — ap, 
term. This yaw acceleration component due to roll 
acceleration gives a favorable yaw component if the 
principal inertial x-axis is below the stability .-axis and 
an unfavorable yaw component if the principal inertial 
axis is above the stability x-axis. This effect on the 
yaw acceleration would have been obscured if we had 
used principal axes. We will recall that the yaw ac- 
celeration differs in the two coordinate system approxi- 
mately by the component ayp. The sideslip response, 
however, must be the same regardless of the axes system 
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selected. The ayp component of yaw rate, which no 
longer exists in principal axes, is compensated for by 
the fact that ap and a,p, in the side acceleration equa- 
tion differ by the same ay component. 

There remains one other type of inertial cross-cou- 
pling moment —that due to engine momentum. These 
engine moments are more closely analogous to ‘“‘gyro 
the 


above and cause a fundamental dissymmetry of aircraft 


scopic’ moments than inertial terms described 


response with respect to direction of roll. 
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(D) Tue Dynamic RESPONSE OF THE AIRCRAFT I) 


ABRUPT ROLLING MANEUVERS 


General 


In the preceding section of this paper we described 
the 
ponents of forces of moments. 


various individual inertial cross coupling con 


In this section we wil 
describe the interaction of these various components 
which produce the rather complex dynamical respons 
during high roll velocity maneuvers. 

Our investigation was made primarily by means 0 
an extensive analog computer facility. Thirty-sj 
channels of recording equipment were used in order 
record each velocity and acceleration term as well 
all of their important components. 

Two significant difficulties present themselves to th 
The first 


is the large number of variations of aircraft configu 


vestigation of the cross-coupling problem. 


rations, flight conditions, and maneuvers which ma\ 


exist. During our investigation, we considered th 
effects of several aircraft configurations as well as many 
flight 


must limit our presentation here to certain representa 


conditions and maneuvers. Of necessity, w 


tive cases. Thus, for example, all the time histories oj 
the 
principal inertial .-axis is aligned with its zero lift axis 


aircralt responses will be for an aircraft whose 


‘In a qualitative way we shall allude to the consequences 
Also, We 


will restrict our presentations to a discussion and com- 


of a change in the principal axis orientation. 
parison of 360° maneuvers performed in low altitude 
high subsonic flight, and initially trimmed at low and 
high angle of attack, respectively. 

The second difliculty is a certain descriptive and cot 
ceptual difficulty which always presents itself when 
ever there exists an intercoupled system of several 

The 
Our description, of necessity, tends to bi 
sequential. With due regard to this handicap we shall 
attempt to explain in turn the response of each of th 


variables. interaction is intrinsically simul 


taneous. 


dynamic variables. 

The response of a hypothetical high-speed fighter is 
first described for a spoiler-controlled rolling maneuver 
with fixed horizontal tail and rudder. To show the 
effect of flight path curvature and initial angle of at 
tack, two rolls, one initiated from straight and level 
flight and the other from a rolling pull out, will be ex 
amined. In addition, in order to show the effects ol 
stability augmentation, a rolling pull-out maneuver 
at the same flight condition is repeated later with yaw 
and pitch damping included on the aircraft. 

Since, for the maneuvers selected here, the effect 
of engine momentum is small and, in addition, rather 
straightforward, it will be neglected in the following 


discussion. 


Response Without Stability Augmentation 
The components of roll acceleration 


Roll Response 


appear in the equation 
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These terms (except the inertial component, which is 


negligible) are shown in Fig. 3 for the roll initiated 


from a high g pull out and in Fig. 4 for the roll initiated 
from straight and level flight. 

Unlike the simple one degree of freedom solution 
which would indicate a roll development until a balance 
is achieved between roll damping and roll command 
moment, the roll response is complicated by roll mo 
ment due to dihedral 1,3. For the case of the high ¢g 
pull-up this term can be as large as the roll damping 
nd roll command moments. In addition, a signifi 
cant effect on roll acceleration is caused by the roll 
due to yaw and the changes.in roll command and roll 
damping moments, resulting from changes in angle of 
ittack during the maneuver. 

Dihedral coupling has a very significant effect in 
either decreasing or increasing roll rate depending on 
the direction of the developed sideslip. The direction 
of the sideslip development, as described in the next 
section, depends primarily on the inclination of the 
principal inertial axis relative to the flight path and on 
the yaw due to roll control. Initial sideslip will be 
positive for positive roll maneuvers when initiated from 
thigh angle of attack attitude. It is generally nega 
tive for positive rolls when initiated from low or nega 
tive angle of attack condition. Fig. 3 shows the posi 


tive sideslip with its corresponding large moment 


opposing roll acceleration and Fig. 4 shows a negative 
sideslip case with a small moment aiding roll accelera 
tion. * 

Some effects of aerodynamic coupling can also be 
seen in the roll acceleration response of Fig. 5. Both 
the roll control effectiveness L;,.6, and the roll damping 
component L,p decrease in magnitude soon after the 
roll commences even though the spoiler deflection and 
roll rate remains constant. The decreases are due to 
the increase in angle of attack. No significant change 
in roll rate occurs for this case since the decreases in 
magnitude of roll control effectiveness and roll damp 
ing approximately cancel each other. 

Side Acceleration Response -Next we examine the @ 
components in Figs. 3 and 4 as indicated in the equa 
tion 


] Pa ie | V3 + V5.5, + (g 1”) sin @ cos @ 


We can see that the total side acceleration as given 
by (V,6 + Y,6,) or [8 + 7 ap + (g 
is small in comparison to the ap and 7 components of 
the total Furthermore, since the 
gravity term is also small we conclude that, if we wish 
to keep 6 small, then r must approximately equal ap. 

Due to the large initial angle of attack in the rela- 
tively high g trim condition, the ap term in the sideslip 


1”) sin @ cos 0 


side acceleration. 


"If the pilot should initiate an abrupt roll while in a zero or 


negative ¢ maneuver, the roll moment, 1g8, may be so great as 


to cause a roll divergence even though the pilot may reverse his 


r, 


ollcontrol 
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equation is initially larger than the yaw rate com 
As a result, a positive sideslip is built up as 
This con 


ponent. 
the aircraft rolls in the positive direction. 
dition is shown in Fig. 3. 

For the low g trim conditions in Fig. 4, however, the 
lower initial angle of attack results in an ap term which 
is relatively smaller than the yaw rate component. 
(Later we shall see that the large yaw rate is mainly 
a consequence of yaw moment due to spoiler deflection.) 
Consequently, a negative sideslip develops. 

Normal Response We now 


the change in angle of attack by examining the com- 
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ponents of 
a=q-— pBp+Z,a + 256, + (g/V) cos ¢ cos 6 


Actually the total normal acceleration is given by 
the sum of the aerodynamic components Z,a + Z;,6; 
(or by the sum of all the pure acceleration components 
including gravity). The change in aerodynamic lift 
term is much more significant in this case than the side 
force in the side acceleration equation, but still less 
than the inertial components of the normal acceleration, 
Bp and gq. Changes in 8p and g must nearly balance 
each other or else large changes in angle of attack will 
occur. 

We note that z/ it were possible to keep sideslip (8) 
to zero, then in order for Aa to be zero, we should keep 
Aq approximately zero. This contrasts with the sideslip 
case where a steady-state yaw velocity is required to 
keep 6 zero. 

Also note that in both Fig. 3 and Fig. 4, the angle of 
attack increases soon after the roll starts. This is due 
to the pitch rate component described later which is 
largely due to inertial coupling, |(7, — /,) J/,|pr. 

The sideslip-roll rate inertia component, however, 
will aid or oppose the angle of attack build-up depend- 
ing on the direction of sideslip. As described in the 
previous section, the initial sideslip during positive 
rolls tends to be positive for the high angle of attack 
maneuvers and negative for the low angle of attack 
maneuvers. Hence, the sideslip-roll rate coupling 
component will tend to oppose a positive angle of 
attack development for the first case and aid a positive 
angle of attack development for the latter case. This 
effect is evident by a comparison of Figs. 3 and 4. 

The term Z,a, being the primary aerodynamic force, 
produces a force on the airplane tending to translate 
it in the pitch plane in such a direction as to reduce 
It is noted that this 
aircraft, 


the change in angle of attack. 
term exerts its effect by translating the 
whereas the two inertial terms g and 8p come in to play 
through rotation of the aircraft about its three axes. 
Since the aircrafts rotational response is more rapid 
than its translational response, the rotation effects 
predominate. 
Pitch Response 


pitch velocity g, we examine the acceleration compo- 


In order to examine the change in 


nents of the pitch equation. 


g = Uf, — 1,)/1,|pr + Mga + Mog + Mua + M; 6, 

In this case the aerodynamic damping term \/,@ is 
negligible, but all the other terms have a major effect 
and are plotted in Figs. 3 and 4. 

The largest term is the inertia coupling term [(/. — 
I,) I,\pr. As described earlier, for most high-per- 
formance fighter aircraft, the inertia ratio is positive 
and approaches unity. Thus, for a large positive roll 
rate and a resulting positive yaw rate, this yaw-roll 
coupling causes an upward pitch acceleration com- 
ponent of the airplane. The positive pitch acceler- 
ation component is the primary cause of the positive 
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pitch rate, which, in turn, is the principal cause ¢ 
angle of attack build-up. 

We note that if we had chosen a condition wher, 
the principal inertial axis was below the zero lift ayjc 
there would have resulted a considerable _pitchiny 
down tendency and a consequent decrease in angle oj 
attack. 

It is significant to notice in both Figs. 3 and 4 thay 
as the roll control reversed on rolling out of the maney 
ver, the yaw rate changes direction at about the sam 
time as does the roll rate. Thus, the inertial coupling 
in the pitch response remains positive throughout th 
maneuver. 

The usual aerodynamic theory indicates that th 
two major aerodynamic components, /,@ and J, 
initially tend to oppose changes of the airplane respong 
from the trim conditions. However, these correctiy 
components lag considerably the disturbing inertial 
component. As a consequence, the large aircraft pitch 
up develops. 

Yaw Dynamics We next study the yaw acceleration 


components, 


r= (JU, —1,) Ilbg + No8 + 


Nw + Nop + Ny8, + Ng 


For the fixed tail case, the main dynamic effect her 
‘is the unbalance among the inertial term, the weather 
cock stability, and the yaw moment due to spoiler. 

As in the usual theory, the weathercock moment is 
always in the direction to cause a yaw acceleration com- 
ponent to reduce the sideslip angle. However, the 
yaw moment due to spoiler, N;,6,, is positive for spoiler 
deflections causing positive roll, but may increase or 
decrease sideslip, depending on other factors. N,4 
causes a component of yaw rate which is in the opposite 
direction to the ap disturbance when the principal axis 
is above the flight path. The effect, for the high a 
case, is to help compensate for the ap disturbance and 
thus help reduce the resulting sideslip. For the low 
a cases, however, the N;,6, component causes a yaw 
rate the dis 


turbance and as a result a sideslip develops in the op 


component which ‘“‘overcompensates” 
posite direction. 

The third significant factor in the yaw response 1s 
the inertia coupling term, [(/, — /,) /.|pg, which for 
a positive roll maneuver and a positive pitch rate pro 
duces a negative component of yaw acceleration. Fig 
3 shows that, for the high trim angle of attack maneu 
vers, the inertial component approximately equals the 
This inertial 


component is the chief factor which prevents the weath 


weathercock moment of the airplane. 


ercocking of the airplane into the relative wind for 
the aircraft with fixed tail controls. 

For the example described the airplane has a positive 
yaw for a positive roll, in both the high and low a 
cases. However, there is a fundamental difference 
between the responses in the high ay and low a» maneu- 
vers. For the high a maneuver, the positive r helps 


compensate for the positive sideslip rate component, 
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». In contrast, for the low a) maneuvers, the posi- 


tive yaw rate is larger than the ap “driving force.” 
that. is, the component 7 ‘‘overcompensates”’ for the 

» component. This “overcompensating” r is 
srincipally to N;5;. 

As a consequence of the different effects of yaw rate 
in the different maneuvers, it is found that, from the 
standpoint of sideslip alleviation for the case shown 
here, yaw damping (Vr), which restricts yaw rate, is 


due 


undesirable for high ay cases but desirable for the low 
, cases. The magnitude of the natural yaw damping 
may be small (as shown in Figs. 3 and 4), but the di- 
rection of the effect of damping is significant when con- 
jdering an airplane with artificially augmented damp- 
More will be said concerning artificial stability 


ing. 


wgmentation in the next part of this paper. 


Response With Artificial Stability Augmentation 


Fig. 5 shows the same high trim maneuvers for the 
same flight condition as in Fig. 3, except that now yaw 
nd pitch stability augmentation has been added. 

The total effect of the dampers is complex and will 
detail here. There are certain 
For the high angle of 


not be discussed in 
ibvious highlights, however. 
attack cases, the yaw dampers oppose the yawing of the 
airplane which is required to minimize sideslip build-up. 
\s a consequence the sideslip is increased in the posi 








tive direction. This represents a major distinction 
between the character of the maneuver at roll rates 
below the ‘‘critical’’ values as compared to those di- 
vergent maneuvers with roll rates greater than the 
critical’’ value. In the latter case yaw dampers help 
prevent divergence in a manner analogous to that of 
the nonrolling aircraft. This has been demonstrated 
by several investigators. One of the consequences of 
the resulting larger sideslip is a decreased roll rate due 
to the larger dihedral rolling moment. 

The primary effect of the pitch dampers is to restrict 
the excessive pitching of the aircraft, caused by inertia 
coupling. Restricting the pitch rate prevents unde 
sirable changes in angle of attack during most of the 
maneuver. 

One of our main purposes in showing this particular 
case with dampers is to demonstrate that caution 


must be observed when one attempts to predict the 





elect of dampers. While pitch dampers usually alle- 
viate the cross-coupling problem, yaw dampers can 


either alleviate it or make it worse. 


AUTOMATIC CONTROL DURING 
ROLLING MANEUVERS 


PART (II) 


ABRUPT 


(A) INTRODUCTION 


Many of the undesirable manifestations of cross cou- 
pling may be minimized by including appropriate de- 
sign features on the original air-frame configuration. 
However, we may penalize the overall aircraft per 
lormance by restricting its maneuvering capabilities 
or by adding the extra weight and drag of enlarged 
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fixed surfaces which may not otherwise be required. 
In some cases the fixed surfaces will not give us the 
best reponse for ail rolling maneuvers and will not 
represent a desirable solution. 

In this part of the paper we will present some basic 
background for another approach, the control of cross 
coupling by automatic means. We shall consider not 
only the functional requirements of such a control, but 
will also consider practical problems such as control de- 
flection limits, measurement difficulties, and reliability. 

A simple technique will be used for determining con- 
trol requirements without considering the stability re- 
quirements. This approach is one way of satisfying 
the need for a simple method to examine control re- 
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quirements and is in contrast to the well known Phillips 
technique which approximately specifies the stability 
requirements without specifying the control require- 
ments. Our approach in determining the required 
control functions is first to define an idealized desired 
roll then to the 
control required to achieve this response. 


maneuver and determine idealized 


This method shows that such controls will have 
rather complex functional requirements in order to 
satisfy our ideal maneuver and at first glance its useful- 
ness may be questioned. However, this method does 
provide an insight into the selection of practical con 
trols which may allow some tolerable degradation ot 
performance as compared to the ideal. In addition, 
it also indicates approximate control surface require- 
ments and provides some insight into the magnitude of 
the major aerodynamic and inertia parameters which 


cause undesirable aircraft response. 


(B) SIMPLIFIED CONTROL CRITERIA 


Ideal Maneuver 


First we define an idealized roll maneuver as one 
whose sideslip, sideslip rate, change in angle of attack, 
(That is, 
This maneuver 


and angle of attack rate are all zero. we de- 
sire to roll about the velocity vector.) 
will not necessarily keep the aircraft headed in a 
straight path because of the gravity effects or, in the 
case of a rolling pull-out, because of the initial curvature 


of the path as well as gravity. 


We will show that, as a consequence of this definition, 
the five complex nonlinear simultaneous equations are 
reduced to five simple equations which are nearly linear 
and independent. In order to further simplify the 
two-force equations we neglect the lmear accelerations 
due to control surface deflections and gravity. These 
terms are small compared to the remaining terms. * 
We then find that the side force equation independently 
defines the yaw requirement and the lift equation inde 
pendently defines the pitch requirement. 

We use each of the three moment equations as fol 


lows: 
pitch equation defines the horizontal tail deflection 


the roll equation defines rolling velocity; the 


requirements; and the yaw equation defines the rudder 
requirements. Each of these equations will now be 
examined using the principal axes system. 


The side force equation reduces from 


V8 + Y;6, + g 1”) sin @ cos 6 


B ap —r+ 


where a ay + Aa, to 


r auP 


* Section (E) points out that the linear accelerations due to con 
trol surface deflections may have major importance when con 
sidering the stability requirements cf the flight control feedback 
loops. However, as a first approach in the ideal case above, only 


the control requirements are considered, 
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and we see the major requirement for zero sideslip 
that we yaw the airplane proportional to the prody 
of rolling velocity and ay (the angle between the prit 
cipal inertial axes and the velocity vector). 7) 
completely defines our yaw requirement. 

The lift equation reduces from 


a Gg = pr + ZO Z 5,5, + (g 1") cos & cos f 


where g = g + Ag, a a + Aa, and 6, = 6, + Ad 


constant 


q LM 


We see that in order to have zero change in angk 
attack we must keep the pitching velocity appro; 
mately constant. This defines our pitch requirement 
Next we will examine the roll equation which 
our ideal maneuver, is approximately linear and ind 
pendent of any other equation. We find the compk 


form is reduced from 


p (J, — 7.) I,-\gr 


where the underlined coefficients vary significantly wit 
angle of attack, to 


p Lpp + Ls,5a 


The term L;,6, is negligible except for low angl 
attack cases where it tends to be significant becaus 
it may change the total rolling velocity as much as? 
per cent. This equation can be solved by inspection t 
give rolling velocity or roll acceleration as a function 
aileron deflection. The transfer function relating r 
rate to aileron deflection becomes the familiar for 


—((L,./L,){1/{1 — (s/L,)]}) 


P(s) 6,,(s) 


or, if we wish to examine the most severe case of a ste 
aileron command, we can express the roll accelerati 
and roll rate as functions of time, 


/ 


p (L;,.52) ” ‘# 


aud 
p (Ls, 5a/Ly) (1 — e”") 


Control Requirements for Ideal Maneuver 


The two force equations and the roll equation cot 
pletely define the ideal maneuver. The pitch and y 


equations will be used to define the control requir 
ments in order to satisfy this ideal maneuver. 


The lift equation gave the criterion that pitch 


velocity had to remain constant in order to keep ze! 


change in angle of attack. The pitch equation defin 


the required horizontal tail deflection to satisfy t! 


pitch requirement. The pitch equation reduces fro! 


q (7, —T1[,) I, |rp + Mya + Maa 4 


Mig + M;,6, — (rep. I 


where a ao + Aa, q qo T+ Ag, and 46, On 7 
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CROSS-COUPLING 


Ad, 6; required 
[d: — LM; 1,|rp + 


The major moment produced by the tail control must 
equal the disturbing moment due to the so-called 
eyroscopic effect which is proportional to the product 
of yawing velocity and roll rate. We note, as explained 
previously, that the so-called ‘‘centrifugal” pitch mo- 
ment terms appearing in the stability axis system do 
not appear here. 

We can also express this required change in hori 
zontal tail deflection as a function of roll rate by using 


the zero sideslip criterion (7 ayp) to give 


Sere [(/ [.)/I, M;,| ach? 4 


(/ P. i. A/ aynp 


ind, in addition, using the solution to the roll equation 
to give the tail deflection required as a function of roll 


control. 
Bicsitha (ay M;,) (Ud: — 1) 1,| 


(IrePe/Ly) | Lieda(s)/Lp[1 — (s/Lp)]}) 


r Z.. 6,(5) Bett 2 
L»)|\? - 


These equations show that the horizontal tail deflec 
tion requirements to satisfy this ideal maneuver are 
reasonably simple. Generally speaking, it 1s found 
that sufficient horizontal tail positions can be obtained 
to eliminate completely this gyroscopic coupling. 

We now use the yaw equation to determine the re 
quired rudder control, remembering that the side force 
equation determined the required yaw rate 1.e., 7 
ap and implies that * = avp. 

The yaw equation reduces from 


de = i, }/i | pq + NB + N,r 4 


N»p ++ N; Oa + N; 6 


+ (Lrebe I:Ns,)qo 


where the underlined coefficients may vary appre 


ciably with angle of attack in the general case, to 
(1/N;.) (aup — N;,6.— | N» + 
(J, —1,) T:lgo + awN,| p) — (Unde T2Ns,)q 


When the values of the various terms in the equation 
above are examined, we find that for a general treat- 
ment none of the terms can be neglected. Referring 
to the rudder requirements above, the following typical 
situations may exist for an aircraft with spoilers for 
roll control: (1) the first term (ap) is the largest 
initially, but the remaining terms predominate as the 
roll velocity increases; (2) for low g rolls at high dy- 
namic pressures the yaw acceleration due to roll rate and 
spoiler deflection predominates; (3) for high g rolls at 
high dynamic pressure the inertia term [(/, — /,) 
I.\qop, the yaw damping term a.V,p, and the yaw due 
to spoiler are all large; and (4) for low dynamic pressures 
the yaw damping term and the yaw due to spoiler pre- 
dominate. 

We thus see that the desired rudder control is a func- 
tion of three separate but related variables— namely, 
roll acceleration, roll rate, and roll control deflection. 


(I eb-/M;, I,)r 
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If we desire to relate rudder to a single variable, we 
can use the simple solution to the roll equation to ex 
press the desired rudder as a function of anyone of the 
three. For example, to express rudder deflection en 
tirely as a function of roll control deflection (assuming 
ailerons were used) we can use the roll equation to 


eliminate p and p to give 


Dy 5 Jrcauivedl 6,(s) [—.N;, Av + Tis) N 1 + 72s) 


where 
Ky —(Ls, N5Lp») (Np + Cao + avr 
i, la (N, Cy au.) 
t= =i 
c Ud, —TI,) 1 


or, if we assume a step roll control command, we can 
write the required rudder deflection as a function of 


time and get 


Two examples of the solution of this equation are shown 
in Fig. 6. 

From the above equation we can see that a step 
aileron input will require an initial step rudder deflection 
ol 


Typical values of this equation are shown in Fig. 7. 
The steady-state rudder requirement is 
ee -6.[+(N5,/Ns,) + (L5,/LpN5,) X 

(N, +4 Cd + aqyN,)|, as t— 


I 


Typical values of this equation are plotted in Fig. 5. 
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Fic. 6. Typical rudder required to hold zero sideslip during 
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Fic. 7. Typical initial rudder required for step roll command. 


The major equations presented in this section are 


summarized in Table 1 and will be evaluated below. 


Interpretation of the Ideal Controls 

The elevator and rudder required to hold zero side- 
slip and zero change in angle of attack may best be 
interpreted by examining the pitch and yaw equations 


namely, 
5; required — [U2 — I,)/1,M5, |aop? 
Di misaibied (1 N5,) Coop — N50 
(Ny + [Ue — 1,)/Telqo + aoN,} p) 
Pitch Control The required horizontal tail deflection 


the 
If we specify 


the 
principal axis and the velocity vector. 


is largely determined by ay, angle between 


TABLE 1 
Summary of Requirements for Ideal Maneuver* 


EQUATICN CONDITION REQUIREMENT 
r useful equation) 
F Sideslip angle ap 
Equals zer 
ange 4 i 
a Atta 
, L-1, 1 Pot | 
- M, i les 
1 
w r Y 1 lap-n, i+a i 
) a I 
2 
Mee Pe % } 
I 
* For the above to have meaning a roll rate must exist. This 


roll rate is determined by the roll equation 


p = Léa + »p 


i 


Cc 
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our rolling velocity and also recognize that we may be 


unable to change the inertia significantly, we see that 
the only way to reduce the horizontal tail deflection 
requirement is to minimize ap. Appropriate distriby. 
tion of mass in the aircraft and utilization of variable 
incidence wings for trim flight are two methods which 
are in the realm of possibility for reducing a» and could 
be considered in the early design stages of an aircraft 
The most practical way of completely satisfying the 
pitch requirement is to control the elevator automat. 
cally. Usually it will be found that sufficient elevator 


position will exist. However, the control required 


varies as an inverse Hence 


in the case of high-altitude flight at the lower speeds 


of airspeed and density. 


the available elevator deflection might be insufficient 


for the ideal control unless some effort were made to 
minimize the angle between the principal axis and the 
velocity vector. 

Yaw Control——In order to see how the rudder should 
be controlled during the maneuvers, we will examine 
raw equation which specifies the re- 
First let 


rudder requirement due to roll acceleration. 


each term in the y 


quired rudder deflection. us examine the 
Fig. 7 
shows typical rudder requirements due to the roll ac 
The 
about primarily because of the sideslip 
that the y 


ts to hold continuously upon entering a roll, the yaw 


celeration component alone. component comes 
requirement 
raw rate must equal app. If the relationship 
the rudder control must also 
We 


rudder deflections required by reducing a as discussed 


acceleration caused by 


equal ay times roll acceleration. can reduce the 




















under the elevator control requirements above, or by 
restricting roll acceleration. 
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Typical steady-state rudder requirement to hold zero 


sideslip during a steady roll 
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CROSS- 


This component due to roll acceleration may call for 
excessive rudder deflections for short periods of time. 
We can slam the rudder against the position limits and 
ve with the sideslip or we can reduce the roll accelera- 
tion. In actual practice the roll acceleration is ap- 
preciably reduced from that shown in Fig. 7 because 
{ the finite rate at which the roll control can be de- 
fected. For example, a deflection rate of 100 deg. sec. 
ata high-speed altitude case may reduce the peak roll 





iceleration by as much as 40 per cent of that required 
ff a step aileron were possible. However, even this 


consideration does not preclude the possibility of 


excessive rudder deflections. 


The major function of the roll control of any aircraft 
isto obtain a given bank angle in as short a time as 1s 
practical. It is possible that roll accelerations can be 
increased in order to produce the desired roll rates. 
However, as we have just seen, large roll accelerations 
may demand large yaw acceleration, hence large rudder 
requirements if sideslip is to be minimized. An ideal 
roll control would command no more roll acceleration 
than that which corresponds to the maximum yaw 


acceleration available from the rudder. 


The importance of roll rate as a criterion for deter- 
mining the severity of cross coupling cannot be under- 
estimated, particularly relating to stability. How- 
ever, considering the problem of controlling the air- 
craft to minimize the disturbances in a stable rolling 
maneuver with limited rudder deflections, the magni- 
tude of the roll acceleration also becomes important. 


Next, we examine the rudder required to oppose the 
moment caused by the roll control deflection, as shown 
mn Figs. 7 and 8. 
quired, hence this important parameter must be con- 
Figs. 7 and 8 show that the yaw 


Large rudder deflections are re- 


sidered carefully. 
due to spoiler is helpful in certain cases and harmful in 
others. The yaw control problem would be greatly 
lleviated if yaw due to the roll control could be varied 


so as to be helpful in all cases. 


The last term to be examined is the rudder required 
to counteract the contribution due to roll rate. The 
various coefficients of roll rate in the required rudder 
equation are shown in Fig. 8. The first term J, is 
primarily a function of airplane configuration and 
varies appreciably with angle of attack. It is one of 
the major disturbances at the high Mach conditions. 
The second coefficient of roll rate ((7, — /,)//-|qo, the 
:0-called gyroscopic inertia term, is primarily a function 
of both the number of g’s the pilot is pulling and the 
true air speed. The third coefficient, a~N,, can be 
varied primarily by the inclination of the principal 
ixis. This last term exists because of the zero sideslip 
condition which requires a steady yaw rate proportional 


to roll rate. 


Fig. S shows that each of these three terms can be 
significant at different times. In the next section, 
methods of mechanizing the ideal control requirements 


just discussed will be briefly presented. 
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(C) MECHANIZATION OF IDEAL CONTROL 


Elevator Control 


It is a fairly straightforward matter to idealize on 
the possible control mechanizations by examining the 
various forms of the equations which define the required 
elevator deflections. We might multiply the output 
of a yaw rate gyro and a roll rate gyro; we might use a 
signal from a roll rate gyro alone, or we might use a 
signal only from the roll control deflection. These 
signals would be combined in a manner which would be 
a function of the source of signals as indicated by the 
different equations. One common problem is that, in 
order to eliminate the disturbance completely, the 
control gain would have to be varied as a function of 
(.\/;,), the moment due to horizontal tail deflection. 
This moment varies primarily as a function of dynamic 
pressure and Mach Number, and gain change com- 
promises would have to be made in a practical mecha- 


nization. 


Rudder Control 


As with the elevator control above, several kinds of 
idealized control schemes can be listed by inspecting 
the various equations which define the required rudder 
deflections. All of the schemes which attempt to op- 
pose the disturbances would obtain signals in one of the 


four methods listed below: 


(a) Utilize roll acceleration, roll rate, and roll control 
deflection. 

(b) Utilize roll control deflection alone plus a net- 
work. 

(c) Utilize roll rate alone plus network. 

(d) Utilize roll acceleration alone plus network. 


All the coefficients vary differently as functions of 
Mach Number, dynamic pressure, and angle of attack. 
The first scheme allows the greatest flexibility since 
each coefficient can be varied independently. Even 
here, in a practical case, compromises in the gain change 
requirements would have to be made. The remaining 
schemes could duplicate the first but the time constants 
of the different networks would have to be varied. 
Generally, it would be desirable that these time con- 
stants be fixed values and hence more compromises 
would be made. 

It is noted that in automatic flight control systems 
it is a common practice to move the rudder proportional 
to roll rate in order to improve transient coordination 
i.e., to minimize sideslip. The equations defining the 
rudder deflection show that if the control is moved 
proportional to roll rate alone, the initial transients due 
to roll acceleration and spoiler deflection as well as the 
effects of angle of attack would not be counteracted. 


(D) MECHANIZATION OF PRACTICAL CONTROL SCHEMES 


We now leave the ideal control schemes and discuss 
more practical approaches to the problem. The best 


practical scheme uses the typical feedback approach 
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of operating on an error signal. This technique par- 
tially circumvents all the many variables listed above 
and attempts to minimize angle of attack and sideslip 
directly rather than counteract each one of the many 
disturbances. We have emphasized the fact that, in 
order to minimize sideslip, efforts to control the rudder 
as a function of any other one variable might only par- 
tially compensate for the disturbances and would re- 
quire complicated gain change schemes. 

The difficult problem of measuring angle of attack 
The 
most common alternative scheme today is to use lateral 
and normal The outputs 


of these accelerometers can be used to control the 


and sideslip has been discussed in many papers. 
acceleration information. 


elevator and rudder so as to minimize any undesirable 


disturbances in angle of attack and sideslip. High 
gain accelerometer feedback schemes might easily 


worsen the flying qualities of the aircraft. Tight ac- 
celeration loops give very large directional and _ pitch 
weathercocking tendency so that some objectionable 
flying qualities might exist in rough air. This problem 
can be avoided by having tight loops only when the 
roll command or roll rate exists. 

In pitch, the use of normal acceleration feedback to 
reduce angle of attack changes may require special 
commands because it is intended that angles of attack 
vary during the maneuver. In this case it is quite 
possible that the best solution is to program the hori 
zontal tail as a function of the disturbing source as 
described in the preceding section. However, in yaw, 
it is desired that the sideslip angle remain zero. In 
addition, the rudder would have to be programmed as a 
function of so many variables that the lateral ac 
celerometer feedback appears desirable. 

A second problem with high gain accelerometer feed 
back loops is that of stability. Poor stability resulting 
from high constant feedback gains can be partially 
overcome by varying the feedback gains with the com 
mand information, such as roll command, roll acceler 
ation, or roll rate. However, it is noted that stability 
in the nonrolling sense and stability in the rolling sense 
must both be considered. In the rolling sense we find 
that the normal and lateral accelerometer feedbacks 
improve stability by raising the roll divergence bound 
ary in accordance with the Phillips’ criteria. In the 
nonrolling sense the conventional dampers used to 
improve stability in nonrolling flight improve stability 
in the rolling sense as well but with certain problems, 
as discussed below. 

In pitch, the conventional pitch damper, which al 
ways opposes any pitch rate, may be helpful during the 
roll. However, a pitch damper canceler, or washout 
circuit, is generally included in order that the damper 
output does not oppose intentional steady pitch rate 
commands. This canceler may appreciably reduce the 
benefit which might be obtained from the pitch rate 
signal in the steady rolling maneuvers unless careful 
attention is paid to the canceler time constant. 

In yaw, the conventional yaw gyro opposes any yaw 


rate. This may significantly aggravate the sideslip 
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when the principal .-axis is at large positive angle 
with respect to the velocity vector since, in this cay 
a steady yaw rate is desired. Careful attention to th 
gain and canceler time constants can minimize this 
effect, but, generally, a significant improvement 

found by eliminating the yaw rate signal completely 
In place of the yaw rate signal we may use the rate oj 
change of sideslip for stability. This can be obtaine 
by differentiating the output of the lateral accelerom 
eter, but extreme care must be taken in the practical 


case to obtain a suitable noise-free signal. 


Considerations involving stability, rough air flying 
qualities, and methods of introducing commands int 
the control loops limit the maximum reduction in sic 
slip and angle of attack which can be obtained by th 
methods mentioned above. Since this is the case, a 


examination of the major disturbing functions 


deseribed in Section (C) will show which of several 
signal sources should be used to aid the direct accel 


erometer feedback schemes. 
(E) PROBLEMS WITH MECHANIZATION 


Control Deflection Limits 


As has previously been discussed in Section (C 
there will seldom be enough rudder deflection availabl 
to keep side-slip zero in the most abrupt rolls wher 
cross coupling is severe. Momentary hitting of the 
rudder deflection limits, such as might occur during the 
initial entry into the roll, will not be as objectionabk 
as when the controls hit the limits in steady-state rolls 
In steady rolls with the control surface against the 
stops, we lose stability augmentation as well as the 


ability to prevent sideslip build-up. 


The major action which can be taken to reduce the 
required rudder deflection, besides the obvious way of 
reducing the rolling performance, is to minimize the 
inclination of the principal axis to the flight path and 
in some instances, minimize the vaw due to roll control 


Gain Changing 


It has already been noted that if controls are pro- 
grammed to oppose the cross-coupling disturbances 
careful attention must be given to the gain of the signal 
or the control may introduce as much of a disturbance 
in some Cases as it Opposes in others. Furthermore 
the variations of the different gains cannot be simph 
programmed as a function of any one variable. The 
tremendous complexity which would be introduced in 
order to compensate totally for all the effects guar- 
antees that, if present gain change techniques are em- 
ployed, controls will only partially compensate for the 


cross-coupling disturbances. 


Measurement of Aircraft Motions 


Sideslip and Angle of Attack--lfi a lateral acceler- 
ometer or pendulum is used as a measure of sideslip 
during abrupt roll maneuvers, a number of quantities 
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not important become very significant. 


mally 
hese are aS fe lows. 
1) A signal proportional to the side force due to 
dder deflection. This appears as a positive feedback 
mal and may cause instability for high gain feedback 
systellls. 

2) A signal which is a function of yaw acceleration, 
a product of pitching velocity and rolling velocity, 


yd the distance the accelerometer is located fore or 








it of the center of gravity. (Distances on the order 
3ft. are important. ) 

3) A signal proportional to the product of roll 
veleration and the distance above or below the roll 
“is. (Distance as small as | ft. may have a significant 
fect.) 

t) A signal proportional to the product of roll rate 
quared and the distance the accelerometer is mounted 
erpendicular to the aircraft plane of symmetry. (Dis 
tances on the order of | ft. may be important.) 

In addition, in the case of a pendulum, the gain of 
the signal 1S proportional to a number of components 
vhich make up the resultant normal acceleration. 

By proper the location, 
these additional components can be used to improve 
By im 


selection of accelerometet 
the stability of the lateral accelerometer loop. 
proper location these same components can result in 
the discard of a potentially useful control scheme. 

A similar situation is true in the case of a normal 
celerometer if it is used to minimize changes in angle 
f attack during these abrupt rolls. The normal ac 
elerometer signal is proportional to elevator deflection 
ind, in addition, may be proportional to roll and pitch 
wcelerations and roll rate squared depending on the 
location of the accelerometer. 

When these signals are meas 
their sensitive axis becomes 


Yaw and Pitch Rate 
ured, the orientation of 
more critical than in most applications because of the 
large roll rates. A misalignment of as much as one 
degree could result in a false yaw or pitch signal due 
to the component of roll velocity and may be as large 
is normal pitch and yaw signals. Proper orientation 
of a rate gyro can be used advantageously to combine 


vaw, pitch, and roll signals in some designs. 
5 s 


F) Fartt-SAFE CONSIDERATIONS 


There appears to be little question that, theoretically, 
the automatic control systems can remove the danger 
ous aspects of the inertia cross-coupling problem. 
There appears to be some question as to whether they 
can completely eliminate undesirable changes in side- 


slip and angle of attack during these abrupt rolls be 














cause of the limited control surface deflections. The 
major point which must be demonstrated is that the 
1utomatic systems can perform this control reliably with 
positive fail-safe features. This can and must be done. 

The subject of fail-safe considerations could well be a 
paper in itself and, hence, no attempt will be made here 
to cover this vital subject completely. However, in 
view of the large control surface deflection requirements 
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during these maneuvers and in view of the reluctance 
by many air-frame manufacturers to accept full au 
thority automatic controls, a few comments should be 
made. 

What is the basie problem ? 
control system which can move the control surface to 
faster 


We need an automatic 


large deflections at certain desirable instances 
than humanly possible, but vet never do this at the 
Actually a finite possibility of malfunc 


any 


wrong time. 


tion for system, automatic or manual. 


Hence, the problem really is to reduce the probability of 


occurs 


malfunction to an acceptably low value. 

This can be done in a number of ways: one practical 
way is combining limited authority devices to obtain 
full authority. The probability of failure can be re 
duced appreciably without 
ample, assume two completely independent channels 


a large penalty. For ex 
acting in series, each controlling one half of the allow 
able deflection, and with the probability of a danger 
ous malfunction as high as once in 1,000 hours for each 
channel. The probability of failure of two channels 
simultaneously would only be once 1n a million (1,000 
Furthermore, if we cousider that, during the 
a condition 


hours. 
total operation, the aircraft would be in 
where a hard-over signal would cause failure of the air 
craft less than 10 per cent of the time, then we get 
only one failure in 10 million hours (100 X&_ 1,000°). 
Practically, the two channels cannot be considered 
completely independent because of the common en 
vironment. The important point here is not the num 
bers, but the principle that with proper attention the 
probability of failure can be made insignificant with 
out a severe penalty. 

This technique is demonstrated by the dual channel 
hydraulic control systems used by most of the fighter 
The reliability is well known 

The electric control channels 
size 


aircraft flying today. 
to the aircraft industry. 


multichannel smaller 


could be made with even 
and weight penalty than the hydraulic systems. 

In addition to reducing the probability of malfune 
tion to as low a value as desired, the concept of using 


limited authority devices in series has additional bene 


fits. The comparison of signals in different channels 
can be used to indicate degradation in performance 
which, while not serious, is undesirable. The error 


signal can be used to warn the pilot or automatically 


take appropriate action. 


(SG) THe ADVANTAGES OF AUTOMATIC CONTROL 


The rapid growth of automatic control systems has 
made the problems associated with reliability, com 
plexity, and cost out of proportion to what they can 
and will be. Under these conditions it is with under 
standable reluctance that the air-frame industry ac 
cepts the automatic control systems of yesterday and 
today. Certainly remarkable improvement in_ reli- 
ability is possible if the effort is so directed and the 
weight and cost penalties are re-evaluated. 


There are considerable advantages to be gained by 
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solving the reliability problem and accepting what the 
automatic controls have to offer. For example, to 
double directional stability of an aireraft by the in- 
clusion of fixed surfaces, we might easily add several 
hundred pounds of weight to the aircraft, not to men- 
tion drag increases which might result. An automatic 
control might do the same job for 20 Ibs. with negligible 
drag increase. If proper emphasis were placed on re 
liability, we could carry two or more parallel control 
channels and their associated fail-safe circuitry plus an 
emergency power supply and still have a savings in 
weight. In addition, we would have the advantage of 
using the controls for other automatic flight control 


functions. 


In the longitudinal case, a similar situation exists. 
For example, if we accept the capabilities of automatic 
control, we can reduce static margins and obtain re 
duction in control power requirements, or we can ob 
tain significant advantages in the types and location of 
control surfaces which can be used. These advantages 
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can be realized if the proper emphasis is placed on re. 
liability and automatic controls are accepted on th 
basis of their potential advantages. 


CONCLUDING REMARKS 


In this paper we have reviewed some of the funda 
mental properties of inertial cross coupling. Speci 
emphasis has been placed on the inertial coupling 
terms which may cause apparently similar aircrafts { 
exhibit different characteristics. Typical time his 
tories of all the components of the aircraft accelerations 
have been shown to give a detailed insight into the dy 
lamiecs resulting from inertial cross-coupling. \\i 
have emphasized the problems of automatic control ij 
muininuzing objectionable characteristics during rapi 
rolls and have shown a simplified method for analyzing 
the control problems. We have shown the advantages 
which the automatic control can offer and have indi 
cated that satisfactory reliability and fail-safe designs 


are possible. 


The Transient Downwash 


Continued from page 740 


ratio should be somewhat less, and for positive immer- 
sion rates, it is believed that the response ratio may be 
somewhat, but not much, larger for cases in which the 
proper phase relationships exist between the response 
to the gust and the response to the downwash. 

It has been demonstrated that the transient down- 
wash occurring when an airplane encounters a moving 
gust field can be of considerable importance in produc- 
Use of the [1 — 
(de/da) | steady-state downwash factor would definitely 


ing stresses in the horizontal tail. 


lead to unconservative results in considering tail re- 
sponse to a sharp-edged gust, and may lead to uncon- 
servative results in the case of graded gusts. 


REFERENCES 


Wagner, H., Uber die Entstehung des dynamischen auftriebes 
von Tragfligeln, Z. Angew. Math. Mech., Bd. 5, Heft 1, February, 
1925 

? Kussner, H. G., Zusammenfassender Bericht tibe 
von Fiigeln, Luftfahrtforschung, Bd. 13, Nr 
12, December, 1936 

Miles, J. W., The 
\oving Gust, Journal of the Aeronautical Sciences, Vol. 23, No 
11, pp. 1044-1050, November, 1956 


den insta 


fionaren auftrieb 


Aerodynamic Force on an Airfoil in a 


Diederica, F. W., and Drischler, J. A., Unsteady-Lift F 
tions for Penetration of Traveling Gusts and Oblique Blast Wi 
paper presented at the Twenty-Fifth Annual Meeting, IAS 
New York, January 28-31, 1957. Preprint No. 714 
tical Engineering Review, Vol. 16, No. 8, p. 47, August, 1957 

5 Hobbs, Norman P., /ndicial Downwash and its Structural I 
fect on the Horizontal Tail, WADC TR 56-164, August, 1956 

® Kuethe, A. M., Schetzer, J. D., Garby, L. E 
R. L., Preliminary Experiments on the Simulation of Gusts 
Wind Tunnel, WADC TN 55-64, Mareh, 1955 

’ Kemp, N. H., On the Lift and Circulation of Airfoils in Si 
Unsteady Flow Problems, Journal of the Aeronautical Sciences 
Vol. 19, No. 10, pp. 718-714, October, 1952 

’ Bisplinghoff, R. L., Ashley, H., and Halfman, R. L., A 
elasticity, Addison-Wesley Publishing Company , Inc., Cambridge, 
Mass., 1955 

’ Reissner, E., Effect of Finite Span on the Airload Distr 
tions for Oscillating Wings, IT, Aerodynamic Theory of Oscillatu 
Wings of Finite Span, NACA TN 1194, 1947 

” Jahnke, E Tables of Functions with Formu 
and Curves, Dover Publications, Inc., New York, 1945 

'! Sears, W. R., Operational Methods in the Theory of A 
in Non-Uniform Motion, Journal of the Franklin Institute, Vol 
230, pp. 95-111, 1940. 

2 Witmer, E., Ruetenik, J. R., Herrmann, W., and Beals, 
V.L., Jr., The Design, Construction, and Operation of the M.1.1 
W.A.D.C. Shock Tube Facility, M.1.T. Aeroelastic and Structures 
Research Laboratory, WADC TR 58-228 (to be published 


Aeronau- 


. and Rox Hs 


, and Emde, F., 








The 


Measuremi 


ed to ¢ ilet 
fach 9.7. FF 


js found tl 
jent is neal 


A’ UT 
New 


Institute 
ided a m 
elocities 
{mong tl 
ameter; 
It./sec. 
spheres wi 
taming sp 
ur. The: 


it./sec. 


| tended at 


spheres fr 
covered t 
The value 

Received 


+ The au 


is helpful 


»Dr. E. J 
Bureau of ¢ 
under Con 
results 

* For 
Now, Phys 


“©! The Drag Coefficient of Very High Velocity 


on the — 


Spheres' 


funda A. J. HODGES 


Spe cial 
upling 
alts t 
e his 
‘ations SUMMARY 
he dy na . P 
Measurements of velocities of two sizes of spheres have been 
We ed to calculate the drag coefficient in the range Mach 2.2 to 


rol in PMach 9.7. Firings were made from a light-gas gun into open air 
jsfound that, at Mach Numbers from 4 to 10, the drag coef 





rapid 
: jent is nearly constant and equal to 0.359 
lvzing 
tages 
ee SYMBOLS 
indi 
‘Signs Le, © coeflicients of quadratic equation 
‘‘ = wire grid in Figs. 1 and 2 
acceleration of sphere 
velocity of sound in ait 
) drag force 
sphere diameter 
Ky drag coefticient 
} mass of sphere 
m initial mass of sphere 
” mass of recovered sphere 
§ transmission line shield in Fig. 1 
grid spacing 
time of arrival of sphere at velocity grid 
sphere velocity 
coordinate of velocity grid 
air density 
IAS 
Nal 
(1) INTRODUCTION 
E 
ge 1948 a light gas gun' was developed at the 
6: New Mexico School of Mines (now New Mexico 
} , . - . . rin ce he 
Institute of Mining and Technology). The gun pro- 
vided a means for accelerating projectiles to very high 
ces elocities for studies of exterior and terminal ballistics. 
{mong the projectiles fired were spheres of 3 S-in. 
liameter; the velocities attained ranged up to 12,000 
ge, . - ; ' 
It./Sec. Drag measurements of 3 S-1in. diameter 
spheres were made, as the gun furnished a means of at- 
F taming sphere velocities of very high Mach Numbers 1n 
ur. These spheres had velocities ranging from 3,250 





it. see. to 10,840 ft. see. The velocity range was ex- 
| tended at the lower end by shooting 5 16-in. diameter 
spheres from a smooth-bored 30-cal. service rifle. This 
covered the range from 2,410 ft. sec. to 4,510 ft. ‘sec. 
The values of the drag coefficient for spheres at the lower 
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velocity ranges agree with values previously published 
by Charters and Thomas? and by Clark and Harris.* 


(2) THe DRAG COEFFICIENT 


The equation for the drag force acting on a sphere 
moving in air may be written as 
D = Kyppd’v’ (1) 
where A, is the dimensionless drag coefficient of the 
sphere, d the sphere diameter, v the sphere velocity, and 
p the air density. With the drag force expressed as the 
product of the mass m and the acceleration a of the 
sphere, a differential equation for the drag coefficient 
may be written as 


(2) 





Kp — (ma / pd*v*) —(m/ pd*)-|(dv/dx)/v 


As shown in Section (4), this equation is useful in 
evaluating the drag coefficient from experimental data. 

There have been many investigations of the drag co- 
efficient of spheres as a function of various variables. 
Among the variables studied have been Mach Number 
and Reynolds Number, and some of the investigations 
of these variables are referenced here. Information for 
sphere drag forces at subsonic velocities has been sum- 
marized in reference 4, in which sphere drag data are 
preseuted for Reynolds Numbers varying from about 
0.04 to 10°. Sphere drag forces at subsonic, transonic, 
and supersonic velocities were studied by Charters and 
Thomas.” The Mach Number range was 0.29 to 3.96 
and the corresponding Reynolds Numbers varied from 
9.3 & 10% to 1.3 & 108 For Reynolds Numbers of the 
order of magnitude 10° they found no Reynolds Number 
effect at supersonic velocities, and they demonstrated 
that the drag coefficient could be correlated as a func- 
tion of Mach Number alone. At low Reynolds Num- 
bers, ranging from 15 to SOO, Kane® has presented data 
showing an increase of drag coefficient by a factor of 
2'/, over a Mach Number range of 2.1 to 2.8. 

In the work reported here the drag coeflicient of 
spheres, as defined by Eq. (2), was measured over a 
range of Mach Numbers from 2.2 to 9.7, using spheres of 


two diameters, 5/16 and 3/8 in. 


(3) EXPERIMENTAL PROCEDURE 


At an open-air firing range steel spheres were fired 
through fine-wire grids spaced at 40-ft. intervals along 
the trajectory. Breaking of a circuit by passage of a 
sphere through a grid produced a signal from a pulser, 
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A photo- 


graphic record was made of the oscilloscope trace by a 


and this signal was fed into an oscilloscope. 
rotating-drum camera. Timing pips on the trace were 
provided by a pulse generator. By measurement of the 
photographic trace, the time intervals of the sphere 
By 


the use of six grids, four evaluations of the drag co- 


flight between successive grids were determined. 


efficient could be made from each round. 
Two different pulser circuits were used, the first of 
which is shown in Fig. 1. When one of the fine-wire 
grids ‘‘a’’ is broken, a timing signal results principally 
from the discharge of the capacitance between one of 
its leads and the grounded transmission-line shield .S. 
The Thyratron fires also, resulting in stepping the trace 
to avoid retracing during successive revolutions of the 
drum in the drum camera. In the second arrangement, 
a sweep circuit was used with the cathode-ray tube, 
thus making the Thyratrons unnecessary and making it 
possible to use the simpler pulser circuit of Fig. 2. 
Before firing, readings were made of the wet- and 
dry-bulb temperatures and of the barometer, in order 
that the air density might be calculated. In the caleu- 
lations of the drag coefficient, the mass and the maxi 
mum diameter (very nearly the original diameter) of 
the recovered projectile were used. The projectiles 
were caught in a backstop consisting of Celotex sheets. 
Drag coefficients were determined from six rounds of 
5 S-in. spheres at Mach Numbers as high as 9.7 and as 
low as 2.9, and from three rounds of 5 16-in. spheres 


from Mach 4.0 to Mach 2.2. 


(4) EVALUATION OF THE DRAG COEFFICIENT 
The time of arrival ¢ of the projectile at some position 
v along the trajectory is determined from experimental 


* The sphere was not recovered, but since the mass loss had been very low for 


m,, for calculation of A, 
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data. It is assumed that the time ¢ can be expressed} 


a power series in the distance + —i.e., 


t F(x) (A + Bx + Cr’) + Dx? +4 


Precision of the measurements warrants use of the ter 


through the quadratic only. Therefore, 


v (dx dt) ~ [1 (B + 2Cyr)]) 
(dv/dx) ~ [—2C (B + 2Cyv)?]} 


Hence, Eq. (2) becomes 


Kp & (m/ pd*)-[2C (B + 2Cy) 


The coeflicients A, B, and C of Eq. (3) can be dete 
mined from experimental measurements by the solut; 


of three simultaneous equations: 


ty A Bx, T Cx," 
to A Bro + Cre 
fs A Bx i Cx3" 


Let distances be measured relative to the center oj 


range of three successive grids spaced at equal interval 


s. Then 
V1 S No 0) \ 5 
B l(t; — ty) (fg — t)] 2s 
( l(ts — te) — (te — t))]/2s 
At the center of the range, Eq. (5) for the drag coefficient 
reduces to 
: 2C om m 2[(ts — tb) — (b - 
A D —_ : a : 
B pil pd s[ (03 — fo) + (bo — I 


and Eq. (4) for the velocity reduces to 


PABLE | 
Round m 7) Loss of Velocity Mach 
No. (gm gm m, % Ft./Sec.) Number 
1S 3.866 3.844 0.57 10,510 9 5S 0.302 
9 400 8.56 Oot 
50 3.8615 3. S610 0.013 8,010 7 .ae 0.388 
7,180 6.56 0). 35 
6,420 5.87 0) 
5,750 5.26 0.39 
51 >. S64 3.857 0.18 7,830 7 a7 0.354 
6,980 6.39 ().ob4 
6,230 5.70 () 358 
5,550 5.08 () 3t4 
55 3. 867 5,410 $. 84 ().357 
$ S50 $ 34 0) 3 
$340 3.88 (). 362 
56 3.865 3.863 0.05 $530 $09 () 
$050 3.66 0. 3t 
3,620 o.24 () oft 
3,230 2 92 0.37 
59 3.864 3.800 1.6 10,840 9.74 0) 
Q 700 8.71 () of 
8,680 7.80 0 
7,700 6.98 ot 
5/16-in.-1 2.036 2.036 1,510 $.()2 0). of 
3,930 3.51 0 3t 
3,450 3.06 0.37 
2,980 > 66 ().37 
5/16-in.—2 2.037 2.037 3,640 3.25 ().378 
3,160 2 82 0.37 
5/16-in.-3 2.037 2.037 2,780 2.49 0.58 
2,410 2.16 (). 3M 


similar rounds, an estimate was 


_ 





made of the sphere mass 
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vy™~ 1/B ~ 2s/[(te — 4) + (tb — te) ] (7) 


The essential data are listed in Table 1, and values of 
Kp are plotted in Fig. 3. For Mach Numbers larger 
than 4, these data may be represented by 


Kp = 0.359 + 0.004 (5) 


where 0.004 is the standard deviation of the values of 
Ky. Assuming that Kp is constant, they may be ade- 


quately described by its average value: 


Kp = 0.359 + 0.001 (9) 


(5) GENERAL REMARKS AND CONCLUSIONS 


Two rounds were fired in order to establish that the 
inertia of the fine-wire grids was negligible. In one 
round the number of wires was doubled; in the other 
round the number of wires was increased tenfold. There 
was no significant effect on calculated values of Ap. 

As indicated earlier, the mass and maximum diameter 
of each recovered sphere were used in the calculations. 
The data presented herein are for spheres which lost no 
more than 1.6 per cent of their mass during their 
traverse of the gun barrel. In general, those rounds in 
which the barrel was wiped clean of oil suffered less loss 
of mass and less roughening of the sphere surface than 
those fired when the barrel was coated with an oil film. 
For spheres whose per cent loss of mass ranged from 1.9 
to 4.7, 


Kp = 0.360 + 0.0085 


One sphere of each category is shown in Fig. 4. 
Fig. 5 shows one of several spark shadowgraphs which 
were made in studying the general aerodynamic per- 


formance of the spheres. 
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The data, as presented in Table 1, Fig. 5, and Eqs, (§ 
and (9), are subject to several uncertainties which ar, 
believed to be of minor importance, but for which py 
strict accounting has been made. Analysis of the ex 
perimental setup and of the magnitudes of the uncer 
tainties in the various measurements to be made make: 
it evident that the greatest single source of error is j; 
the measurements of time. The magnitude of the larg 
est error in A, resulting from inaccuracies in the tiny 
measurements may well be in excess of one per cent 
This error is believed to account for nearly all the devia 
tions of calculated individual values of Ay from th 
mean. 

Drag coefficients of spheres in air were determine 
at Mach Numbers ranging from 2.2 to 9.7. 
indicate that at Mach Numbers ranging from 4 to 10 the 
drag coefficient of spheres in air does not vary signifi 


cantly. 
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A Theoretical Investigation of Laminar 


Separation in Supersonic Flow 


G. E. GADD* 
National Physical Laboratory, Teddington, England 


SUMMARY 
The effects of Mach Number, Reynolds Number, wall tem 
erature, and surface curvature on laminar separation in super 


sonic flow are investigated by an extension of Stratford's method 
Two-dimensional flows only are considered; where the surface is 
The wall 
temperature is taken to be uniform, although this restriction is 


urved the generators are perpendicular to the flow 
not essential. It is assumed that separation is of the usual prac 
tical type for supersonic flow, taking place well upstream of the 

ck wave or other agency which provokes separation. In these 
ircumstances the flow in the neighborhood of separation is en 
tirely governed by the factors considered in this paper-——namely, 
the Mach Revnolds the 
urvature of the wall. The results of the analysis for characteris 


ind Numbers and temperature and 
tics of the flow such as the pressure distribution at the wall agre« 


well in most respects with experiment. 


SYMBOLS 


= distance from the leading edge measured 
along the wall 
= distance from the wall measured normal to 
the wall 
= distance along a streamline 
velocity component in the x direction 
= velocity component in the < direction 
= static pressure 
= p between the leading edge and the position 
O of Fig. 1 or Fig. 2 
u = viscosity 
= density 
enthalpy (or ‘total heat” 
= Prandtl Number 


/ = | (u?/2 
ez 
y = muss flux puds between the wall and the 
/7 0 
point considered 
y = y at the join between the timer and outer 
profiles 
= ratio of the specific heats 
VW = Mach Number 
R = Reynolds Number ppg, Up, X/p,; 
T = viscous stress u(O0u/Odz) 
= nondimensional measure of the viscous stress 
in a constant-pressure boundary layer—set 
Eq. (12 
= velocity ratio u./u,, , in a constant-pressure 
boundary layer 
A,B,C, G defined by Eqs.(14) and (15 
m,0,.S defined by Eqs. (18) 
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be = coefficients.in the assumed relation (16) for 
the inner profil 
€ = (p/p ] 
x defined by Eq. (18’), Section (2.2 
y = (AK 1 ay ipproximately equal to the 
thickening of the boundary layer 
L, 4 defined by Eq. (26 
defined by Eq. (28 
G = pressure coefficient (2/7M,. .*)|(p/p l 
a = exponential factor in the upstream variation 
of C, with x 
W = (K — 1)Dj 
m = surface curvature parameter, approximately 


the reciprocal of the radius of curvature 

Subscripts 
( denotes conditions that would obtain at the general dis 
tance x from the leading edge if the pressure remained 


constant everywhere, equal to p 


E denotes conditions at the edge of the boundary laver 

W denotes conditions at the wall 

y, ¥ denote that the mass flux between the wall and the point 
considered is y, ¥ 

0 denotes conditions at the position O of Fig. 1 or Fig. 2 

s denotes conditions at the separation point $ 


(1) INTRODUCTION 


a SEPARATION in supersonic flow can be pro 
voked in a variety of ways, for instance by an 
oblique shock wave incident upon the boundary layer 
or by a step in the wall. In most practical cases, how- 
ever, Separation takes place well upstream of the agency 
provoking it. In these circumstances, the pressure dis 
tribution in the region of separation is of the ‘‘simple 
wave’’ type, the external flow being deflected from its 
original free-stream direction by the thickening of the 
boundary layer. This thickening is in turn governed by 
the pressure distribution, and the equilibrium between 
these processes therefore determines the local flow. 
The present paper considers the effects on such sepa 
rating flow of the main relevant variables--namely, 
Mach Number, Reynolds Number, wall temperature, 
and surface curvature. It is assumed that the flow is 
two-dimensional, and that where the surface is curved 
the generators are perpendicular to the flow, as with a 
two-dimensional airfoil. The wall temperature is as- 
sumed to be uniform, though in principle the method 
can still be applied if the wall temperature varies from 
point to point, provided that its gradient in the region 
of separation is not too great. Similarly, although it is 
assumed for simplicity that the Prandtl Number ¢ is 
equal to | and that viscosity is proportional to absolute 


temperature, these assumptions are not essential. 
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Flat wall~ Agency, such as a step in the 
wall, provoking separation at S 
Xo 
of 
S 
Pe - 
—_—_ 
x 
Fic. 1 The flow pattern and general shape of the pressure 


distribution for separation from a flat wall 


The method of analysis is an extension of Stratford's 
method! for incompressible flow. A similar approach 
has been used before by the present author,’~! but the 
treatment given here is more complete and makes fewer 
approximations. 

The fundamental principle of Stratford’s method is 
that the boundary layer is divided into two regions, one 
adjacent to the wall, and the other the remainder of the 
boundary layer. The flow in this latter region can be 
shown to be essentially inviscid when the adverse pres 
sure gradients are fairly sharp. Thus the shape of the 
outer part of the boundary-layer mass-flow profile can 
be predicted approximately by what are, in effect, in 
viscid flow considerations. The inner profile is deter- 
mined by the conditions that it must join smoothly on 
to the outer profile, that continuity must be satisfied, 
and that, at the wall, the rate of change of viscous stress 
normal to the wall must balance the pressure gradient. 
In the supersonic case the further condition is imposed 
that the pressure is related to the external-flow deflec- 
tion brought about by the thickening of the boundary 
layer. By applying all these conditions a complete ap- 
proximate solution to the problem may be obtained. 

Section (2.1) presents the analysis in a general form 
for the flat wall case. This analysis gives quite simple 
results, as in Section (2.2), when the inner profile is as- 
sumed to be only a small fraction of the total. Such an 
assumption is theoretically justifiable only at very high 
Reynolds Numbers, although the results agree sur- 
prisingly well in many respects with experimental data 
An account is 


given in Section (2.3) of some results computed with 


for much lower Reynolds Numbers. 


the inner profile considered as forming more than a 
small fraction of the total. In Section (2.4) the analysis 
of Section (2.2) is extended in a relatively crude and 
simple way to deal with the effects of surface curvature. 
Finally, in Section (3) all the results are compared with 
experiment, and the advantages and limitations of the 
method are discussed. 

ANALYSIS FOR A FLAT WALL 


(2.1) GENERAI 


2.1.1) The Outer Solution 


The boundary-layer momentum and energy equa- 


tions in two-dimensional flow are 
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pu(Ou Ox) + pw(Ou/0z) = — dp dx) + radi€ 
O Oz)|u(Ou Oz eloci 
siti 
pu (O//Ox) + pw (O//0z) = u (dp dx) 4 would 

(0 /0s)[(u/o)(O/ Oz)] + u(Ou dz)? (9 - 
If the Prandtl Number o is assumed to be equal to 
multiplying Eq. (1) by w and adding it to Eq. (2) giv ccep! 

pu(Oh Ox) + pw(dh Oz) 0 Oz)|u(Oh Oz a 
sharp 
where right- 
site 
h I + (u?/2) 
uter 
Eqs. (1) and (3) can be written spect 
i ne ve 
Ou?/Os = —(2/p)(dp/ds) (2 p)(O Os)[u(Ou Oz n ord 
Oh Os = (1 pu)(O/Oz)[u(Oh Oz) | ” 
ure 1 
where O Os denotes the rate of change with distan uter 
along a streamline. uch 
For separation from a flat wall the pressure distriby n apy 
tion has the general shape shown in Fig. 1. If the in-  presst 
creases in pressure near the leading edge* due to tl Fro 
rapid rate of growth of the boundary layer there ar 
ignored, the pressure can be treated as constant, equi 
to p., over a length x, from the leading edge to som vhere 
point O downstream of which the pressure begins t 
rise. (In theory, the pressure rises exponentially, s 
that there is no precise point O, but in practice O is de | But, 
fined well enough.) Consider the change of / along 
streamline between O and some general point P at 
distance x from the leading edge, where x > x,._ Ii th 
outer velocity and temperature profiles do not chang 
too markedly in shape between the two positions, th 
right-hand side of Eq. (5’) in the outer flow at positior 
P will be of the same order of magnitude as the value it By de 
would take if the pressure between O and P remaine 
constant at p,. Thus, at P in the outer flow, 

Oh, Os = Ofdhy ./Os] \lso, | 
where y is the mass flux between the streamline con- | 
sidered and the wall, and the subscript c denotes con _— 
ditions that would obtain at the general distance x fron ey 
the leading edge if the pressure remained constant ever) 
where, equal to p,. Also of course at O 

Oh, Is = Ohy ,/ds 
Hence, integrating along the streamline gives rhe « 
hy» — hy =h —hy . 11) i 
pro 
hy —hy x i 7 
oe O- 
Probably |h, , — hy will be small compared with Pa 
h, , and hence, to a close approximation, , 7 
‘qs 
hy, =h, . 
2.1.2 
Now consider Eq. (1’), which expresses the fact that Lay 
uw” decreases partly because of the adverse pressufe hey 
distin; 


* These become important only at hypersonic speeds 


Oz 


jual to 


hy as 
2) give 


dz) ] 


distanc 


distribu 
f the it 
e to tl 
here ar 
It, equi 
to som 
egins t 
ially, s 
O is de 
along 
P at 
If the 
chang 
ms, the 
positior 
value it 


maine 


t that 
essure 
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radient and partly because of friction. If the outer 


elocity and temperature profiles at the general 


osition P are not very different in shape from what they 
would be 


nates té 


at constant pressure, the equation approxti- 


Ou, OS —(2/p)(dp ds) 4 On Os) 6G 


xcept near the wall. At the wall the viscous forces 


st balance the pressure gradient forces, however 
sharp the latter may be—1.e., the two terms on the 


right-hand side of Eq. (1°) must be equal and oppo 


site so that the approximation (6) clearly fails. In the 
uter flow, however, it is quite a good approximation, 
especially when the pressure gradient is sharp. Then 


the velocity profile must be sharply curved at the wall 
in order that the balance of viscous and pressure forces 
in be satisfied there. The region of very sharp curva 
ture must, however, be limited in extent, so that in the 
uter flow the pressure gradient term in Eq. (1°) will be 
uch bigger than the friction term, and any errors made 
approximating to the friction term by its constant- 
pressure value will be unimportant. 


From Eq. (6), since p = [y (¥ 


j A {a 


] o ft 
ad log p) as y y 


But, from Eqs. (4) and (5), 


I, / + (Vy, 2/2) 


F 6| pl), 
O( pl)» pit) ae 
(pl), ‘i. = 


By definition 


éy = pu) 62, 


IN SUPERSONIC FLOW 761 


If the wall temperature is uniform, /, 


a small amount with x between the position O and, say 


will vary only 


the separation point, so that it can be treated as inde 


pendent of x or s. Hence 

d (log Pp) ds ¥ y¥—!1 ad ds } log # t | y 2); 
PP ; | + Uy fy rt ‘ 

Similarly, 

P Py « Pf ‘ 
so that the outer flow behaves in a quasi-inviscid man 
ner. Since / u | (y — 1).1/°, it follows from Eqs 

7) and (S) that 
9 ' 2 
1) bf 4 — a 
plu), = p p = x 
( + | /, { 
yj 


in the outer flow. 
Consider now the flow on an adjacent streamline, be 


tween which and the wall the mass flux is Y + 6y, where 
6y < wy. For this streamline write 
pit), y piljy 6 pu 
pl), pu O| pul 
uM, y Vy + 6\/, 
S 62 
z g 6 


Then, from Eq. (9) 


(pu), 62, 


\lso, for the constant-pressure boundary layer, from Eqs. (1), (3), and (4), 


I. + (u.?/2) = In 


where subscript Il’ denotes conditions at the wall, and subscript / conditions at the edge of the laver 


some algebraic manipulation, it follows that 
bad ( p 
| @ 3 p 
| 


‘he expression in braces on the right-hand side of Eq. 


l)M, 


is near to unity eveiywhere in the outer flow, so, 


ipproximately, 


O-( pit) (pl), O( pit) Oz], wa} 12 
| = (12) 
L Oz (pu), - [O(pu) Oz], Oz* Jy 

Eqs. (9), (11), and (12) constitute the outer solution. 


. . nae ; is 
" 2) Conditions in the Constant-Pressure Boundary 
ayer 


In the outer solution we need to evaluate the terms 
distinguished by subscript c, corresponding to condi- 


(he — hw)(te uy 


*(dwily ) + My — 1) 24M 


10) 


Hence, after 


tions which would obtain at the general distance x from 


the leading edge if the pressure remained constant 
everywhere, equal to p,. 

The solution for a constant-pressure boundary layer 
is relatively simple if the wall temperature is assumed 
to be uniform, the Prandtl Number o to be equal to 1, 
and the viscosity to be everywhere proportional to the 
absolute temperature. The latter two assumptions are 
not greatly in error, while that regarding wall tempera 
ture can well be approximately true under experimental 


conditions. Hence all three simplifying assumptions 
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TABLE | 
Nondimensional Viscous Stress ¢(n 


n 0 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 5 
t(n-) 0.332 0.332 0.332 0.331 0.330 0.328 0.325 0.321 0.316 0.309 0.300 () 99 
n 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 . 
tn 0.276 0.260 0.241 0.218 0.191 0.159 0.119 0.070 


are made here. In these circumstances, for all wall temperatures and Mach Numbers, the viscous stress ; 


in the constant-pressure layer is given’ by 


9 1/2 
To = pte(Ou/OS)- = (pg ce. c*/Rr ’“) (ne) 
where Ne = Ue/ UE 
R = PE cl} V ME 


and ¢ is a known function of n,, given in Table 1. Also, from Eq. (10), 
(ue ww) = U./Iw) = 1 — [(hw — he) Iwilne — U2 Tw)((y — 1)/2).M; 7 


Hence, after some manipulation, we can write 


(pu); = pw Up cA(ne) 
[O(pu)/ Oz]. = pw. ctte. Cle -/ Iw)(R, = x)B(n-) 
[O*( pi) Oz? |. = pw, UE. lie. e/tn *)(R, x*)C(n ) 


zo = (wile, (x R.’ *)D(n-) 


y = [(ou)atz = pw. clr, &l(n-) 
J0 
M. = Me le. Ty)" *F(n-) 
Iw/I. = Gln.) 
where 
A = 7G 
B = G*thl + Ue -/Iw)[(y — 1)/2) Me, <n} 


c = (yr L)(/; ? Iw) My ols *t*n, — BG(dt dn-) . 3G°Bty [(hu aes hz) Tw | - (I; Ty is = 1); net 


p= | >] = (hy — hy) Twn — (Jer, Iw)lQ — 1) 2).\; =n (dn f) 1: 


E = (1/D) | (n dn, t) 
~/7V 


PF = ae” 


Se ;1 — (hn = hy) ly In. = (Ie. . IwlQ — 1) 2)M, on 2 


(2.1.3) The Inner Solution 


The outer solution given above ceases to be valid near the wall. Denote by % the mass flux between the wall and 


the innermost point at which the outer solution is valid. Then for y < ¥ a different solution, the inner solution, must 
be used. The junction between the inner and outer solutions will not itself be a streamline, so that ¥ will vary with 


For the inner solution, assume 


(pu), = pw. clr. -la(z/2z) + b(2?/257) + c(23/2g*)] + (pu). = 2/K.c | 
? (16 . 
where 23 = Kzj.. 


Y in pu, O(pu) Oz, and O°(pu)/ dz”, and the 


ee) 
— 


The inner solution is made to be continuous with the outer solution < 


further condition is imposed that 


* Thus, when s = sj, the inner profile relation (16) reduces to 


(pu)y = pw. Up, (a +b +c) + (pur 


It foll 


where 


From 


where 
Se] 


since 


\ssut 
begin 
bect I 


the a 


Th 
terna 
flow 1 
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vary 
edge 
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angle 
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thick 
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y = } " pu ds = | , (pu) dz, 
~ WV 0 
) 0 55 
| 0.2% § [t follows from Eqs. (9), (11), (12), (14), and (15) that 
ar b+ i Aj Ny 
a+ 26+ 3c + (BD); = AQ; 
17) 
Stress - 2b + 6e + (CD?)5 = K?S; 
(a/2)+ (63) + (cC4 + hj = kVA 
where 
y (£) ), 21(p pe) ' — Hl , 
p.) (y—1)Me le F? § ~ 
2 l _ ‘ - : = as 27) 2 ] rt 1S) 
0 (*) (1 + (p’ pe) 1})1 = [(+ 1)/2|Me .? Ud; re) BD 
p \¥Y — 1)\M,; 2; Iw)F2;G+ [(y - 1) /2).M;, (I; Tw) F* | 


S = NQ(CD AB) 


From Eqs. (17) 


a = 3(N — A) — 2(AQ — BD) + [(K?*S — CD*)/2] 

b = —3(N — A) + 3(KQ — BD) — (K°S — CD?) 

c = (N — A) — (AQ — BD) + [(K2S — CD*)/2] - 
| (1 24)K3S — (1/4)K2O0 + [(3/4)(N — A) + (1/4)BD — (1/24)CD? + EJK — E =0 


where all the terms are evaluated at y. 
Separation occurs when O(pu) Oz = 0 at the wall—.e., 


when 


a = —0.332D; 20) 


since at the wall B = ¢ = 0.332. 


(Cf. Eqs. (15) and Table 1.) 


2.1.4) The Pressure Gradient Condition and the Final Equations 


From Eq. (1) 


dp dx = (O Oz) [u( Ou Oz) |n = (uw/pw, c) (Pe p)) [O?( pi) Oz" | + (3/I)(o/ Oz) |O( pi) Oz]i un 


\ssume that (O//0zs)w = (O//0z)n 


begins to rise (cf. Fig. 1), but it becomes progressively less accurate downstream. 


This assumption is of course accurate at the position O where the pressure 


However, at separation O( pu) / Oz 


becomes zero at the wall, so that errors made in (O/ 02) become progressively less important downstream of O, and 


the assumption for (O/ (Oz) is therefore good enough. 
d ( p ) YM, Lr.c Pej 26 
dx \ p. x Iw p (K2D3? 


The thickening of the boundary layer deflects the ex- 
ternal flow and it is because of this deflection of the 
flow that the pressure rises above its upstream value /.. 


0.996(hw — hela 0.331 (hw — he) (K — 1)) 
KDjlIy ly Kk? { 


Hence from Eqs. (10) and (13) to (16), and Table 1, 


the thickening of the boundary layer in a precise way, it 
hardly seems justifiable to express the relationship by a 


ill and 
ial Where the pressure varies along the wall, it must also very complicated equation in the analysis. At first 
vith x. vary along lines normal to the wall, because toward the sight the obvious relation to use would be the one based 
edge of the boundary layer the pressure tends to be con- on equating the deflection of the external flow to dé, dx, 
stant along the inclined Mach lines. Similarly, the flow where 6; is the displacement thickness. However, this 
Pe angles are not constant along normals to the wall. All leads to complicated results, and a much simpler ex- 
(16) this is, strictly speaking, incompatible with the preced- pression can be obtained by considering the inner flow, 
ing analysis based on the boundary-layer equations where most of the increase in boundary-layer thickness 
id the which assume that the pressure does not vary across the takes place. The distance from the wall, or ‘ordinate, ’’ 


thickness of the boundary layer. However, since the 
longitudinal variation of pressure is spread over many 
boundary-layer thicknesses, the assumptions are not 
seriously violated at moderate Mach Numbers. All 
the same, because (p — p-) cannot easily be related to 


of the outermost point of the inner flow is 23, and 27, , 1s 
what the ordinate of the point with the same mass flux 
Y between itself and the wall would be if the pressure 
Thus the thickening of the 
(K — 1)zj 


remained constant at p,. 


inner flow is y = 23 — 23. = by Eq. (16). 


“\ 
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For s > 23, the outer solution applies, and in this region 
the streamlines diverge slightly where the local Mach 
Number is less than | and converge slightly where it is 
greater than 1. The net contribution to the thickening 
of the boundary layer is small, and the angle of deflec- 
tion of the external flow due to the thickening of the 
boundary layer can accordingly be taken to be dy dy. 
Since this angle is small the simple wave flow relations 


give 
(p/p) — 1 = ly; (Mp 2 — 1)'?](dy/dx) (22) 
Hence, by Eqs. (1+), 


(d/dx)\ ((K — R,'’"|D;} = 


[(.\MW, -— 1)? /y Me 2, Iw) |(p i) oe | | 


or approximately, since the term x Rx!” in the dif- 
ferential varies only by a relatively small amount be- 
tween the position O of Fig. 1 and, say, the separation 
point, 


(d/dx){(K — 1)D3] = 2 —1)'"/yMe, 2] X 


Jr 


[( Mr 


(23) 


Iw)(R,? x) [(p/p.) — 1] 


From Eqs (14), (15), and (18) to (25), the complete 
approximate solution for the flow in the vicinity of 


separation may be obtained. 


AERONAUTI 


re 
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(Cf. Table 1.) Similarly, if the proportional increase 
pressure, €« = (p p.) — 1, is assumed to be very sma 
of order 7,”, relations (1S) become 
N = 7-(1 — x), 
where x elw/yMe. Le. ene’, R 
QY=RH=n, S = 0 
Hence the equation for A in Eqs. (19) becomes 
—(1/4)AK? + (38/4)A(1 — xz) — (1/2) = 0 
so that XJ (A 1)\(2 — A)/3A 
and from Eqs. (19 
a = —([(A — 1)(K + 2)/K]nj a 
b = [2(K — 1)(K + 1)/K Jn 


Hence, from Eq. (20) separation occurs when 


K = Kk, = v2 95 
Eq. (21) becomes 
(d/dx)\ |((K — 1)(2 — A)/K]nz. 2} = 


(1.323 /x)[(A 


while Eq. (23) becomes 


(d/dx){((K — 1)nz. -] = 0.11101, — ]) 


99 Tk > . _— r SAY 7 - ‘LENE ew > , : ) . - F 
(2.2) THE SOLUTION WHEN THE THICKNESS OF THE (I; Iw?)(R?/x) [(K — 1)(2 K)/K |n 28 
INNER PROFILE IS Mucu Less THAN THE TOTAL 
BOUNDARY-LAYER THICKNESS The two unknowns, A and nj, ,, can be determined 
. , er - from the two equations (21’) and (23) as follows. Fron 
For small values of ., relations (15) simplily to each equation we can obtain expressions for dA dx it 
1=7,, B = 0.332, terms of A, nz ., and dnj ./dx, and for dny, ./dx it 
as ” terms of A, nz ,, and dA /dx. Hence, by elimination, | 
C = 0.331 [(hw — he) Tw], (15°) expressions can be obtained for dA dx and dnj, , dx 1 
D = 7./0.332, E=7./2, F=n, G=! terms of A and nj, only. Then, by dividing one rela 
=n 32, = 7./2 = 7,. = — ‘“ ; 
tion into the other we find 
dnj nz. [1.323(K + 1) — 0.111(Mez 2 — 1)? Up 2/Iw?) R,'7(2 — K)(2 - K*)nj. .*] 
dK (K — 1) [—1.3823(K + 1) + 0.222(Me 2 — 1)'? Ue. 2/Tw2)R'? K(2 — K)2nz, 4] 
Hence, if 
nyc = (Me 2 — 1)7°/R,~ (Iw /Te. -)'/*L 
where LA 
dy 1V[1.323(K + 1) — 0.111(2 — A)(2 — K?)Y] = 
dK (A — 1)[0.222A(2 — A)?Y — 1.323(AK + 1)] Z 


Upstream, where p/p. —> 1, A —~ 1 by Eqs. (18’) and 
(24). Hence the numerator of the right-hand side of 
Eq. (27) must become zero, so that Y ~ 23.82 upstream. 
Similarly, by differentiation of the top and bottom of 
the right-hand side of Eq. (27), it follows that dV dk > 
66.70 upstream. With these starting conditions it is 
possible to solve Eq. (27) for Y and L as functions of A 
by numerical integration—by, for instance, repeatedly 
applying the approximate relation 


S ies Ki 2 | x 
(dV /dK)x, + 


VKnos = YR, + IK, 


(dV dA )k 
and solving Eq. (27) at each stage for (dV/dA)& 
The approximate solution, so determined, is given 1 
Table 2. 


From Egs. (23’) and (26), 
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TABLE 2 
rhe Solution of Eq. (27 
0 1.12 1.22 1.30 1.37 1.48 1.48 1.52 1.56 1.59 1.62 1.65 ay 1.69 1.71 1.73 1.75 
21 2.41 ».61 2.80 3.00 >.20 3.38 3.00 3.79 3.91 $10 $59) $40 $63 +S] » OL 5.26 
441.80 2.17 2.55 2.99 >. 48 >. 49 $50 5.13 B. 7a 6.37 7.31 7.99 8.87 9S] 11.02 12.26 
PABLE 3 
( V | tR.1 4 asa Function of 
16 14 12 10) S 6 } 2 () y } 6 S 10 
1 4% t () 09 0.13 0.20 0.29 (). 42 56 0.74 0.94 : 39 1.32 150 1.68 S4 1.98 
K-1)\(dL dK) + LI(dA, dx) = O111CM; “we which by Eq. (25) takes place when A 1.41. From 
; Eq. (32), the pressure gradient increases as the wall tem 
R,°'?/x)(1 Iw)’ “(AK — 1)22 — ADL? A] | : mo 
perature is reduced, and since the pressure coethicient at 
Hence, 1 separation remains unaltered, this means that cooling 
COV 1)3 8p 387 : 2 19° in effect makes it more diflicult to separate the boundary 
Tae) i) af l C iw) =-) _ bap ‘ “ ‘ : 
layer. For if the pressure distribution could be speci 
so that fied independently of the rate of thickening of the 
/ 8 R38 /y)(] : : boundary layer, an increased pressure gradient after 
ay \ V)ULy iy i i . 
; an unchanged upstream region at constant pressure p 
it follows that would lead to separation at a lower pressure if other 
conditions were unchanged. 
dk O11 A — 1)(2 — A)L ‘ : , ; 
; ; di F 29) Upstream, in the region of the point O of Fig. 1, 
dr A((K - IL di - i ; , ee a 
. ie head K > 1 and Y — 23.82, as stated above, so that L > 2.21 
Since L and dL dK are known as functions of A from and Eq. (29) approximates to 
Table 2, A can be determined by numerical integration 1» ide 0.245 
. a if oa OR, 
isa function of 7 -., Where 7, is the value of 7 at the 
. ° ° - a , I ~oen 5 0.245 
separation point at which A v/ 2 by Eg. (25). Also, where k = A 1. Hence k«« , where r, 1S 
vl F F the value of r at the position O, and from Eq. (30), C 
— = Y | F { { 245 1 
: ‘ , . w) X e” also. Thus, upstream, 
(AK — 1)(2 — A) 3A] 73 . . 
| 7 ( ( 
by Eqs. (1S’) and (24). Hence by Eqs. (26), the pres- where 
sure coeflicient . — 
a = 0.245[(r — r.)xo/(x — X)] 
C 2e 2L°~AK — 1)(2 — A) - 0.245(Me 2 — 1)*/R,2/*(Ig -/In 
= OU) 
yMe 2 3(Me,2 — 1)'*R,"7K 
; It has been assumed that the thickness of the inner 
; l/4p 1/4 “es , 
so that C,(.1/; 1)’ *R can be determined as a profile is much less than that of the total boundary 


function of r — r,. Results are presented in Table 3. 


\t the separation point the pressure c¢ efficient 1s 


‘| Medea! "alii 31) 


c = L.13(M; = 


and is independent of wall temperature. 
The gradient of C,, dC, dx, is given by 
\ dC, dx) 


1) SRS . Tw)? 7*(dK dr) X 


(I dK)(L°(A — 1)(2 — A)/K] — (C, 4) 


The maximum value (dA dr) X 


by Eqs. (28) and (30). 

ddk){(L7(K — 1)(2 — AK) K] occurs when K = 
1.37 and C 1.02(.M/; 1 tie! 
equal to 0.151. Hence, the maximum value of 


aC dy 


— 1S 


4 
, and 


is approximately 


0.101(Me 2 — 1)'*R,"(Ie ./Iw)*’* (32) 


This maximum occurs just upstream of separation, 


layer. This means effectively that nz. < 1 4, say. 


Hence by Eq. (26), 
| ing 1M; 1/8) (7 270 


1) I; 


At the upstream position O, Z is equal to 2.21 and at the 
separation point it is 3.15. Hence, if the wall tempera- 
ture is of the same order as the zero heat-transfer value, 


so that 

Iw/T = O}1 +4 y¥ — 1)/2) 1; ‘ 
and if J, . isin the range 1.2 to 5, the Reynolds Num 
ber R, ought strictly to be very large indeed, 10° at 
least, for the results of Table 3 and Eqs. (31) to (33) to 


be valid. 


THICKNESS OI 
VERY SMALI 


(3.3) 


THE INNER PROFILE NOT ASSUMED TO BI 


COMPUTED SOLUTIONS WITH THEI 


For inner profile thicknesses that are not small, the 
full Eqs. (15) and (18) should, in theory, be used rather 


than the approximations (15’) and (18’). The solution 
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becomes much more complicated, and at an early stage (AK — 1) (dDjz dnj, .)(dnz, - dx) + Djz(OK / Onj x feren 
in the calculations particular values must be assumed . scatt 
. i I : (dng. ./dx) + Dj [OK /O(p/ p-) |[d( p/p.) /d = 
for the Mach Number, Reynolds Number, and enthalpy ' ; 
ratio Jy Jp. Cases have been computed with since from Eqs. (15) Dj is a function only of nj . whik | crept 
My . = 3, with R,,, the Reynolds Number at the posi- from Eqs. (18) and (19) A is a function both of nz . and) been 
tion O of Fig. 1, equal to 0.35 & 10°, and with wall of p p.. Upstream, in the region of the position ( IS) 1 
temperatures of 245°K., 275°K., and 345°K., the stag- where p = p,, A as given by Eqs. (19) takes the limiting © of 
nation temperature being 275°K. These conditions value | for all assumed values of nj... Hence, 0A On; rh 
are close to those at which experiments have been done.° — (0, while 0K /0(p/p.) -|> 0 upstream. The terms) ponet 
The method of solution is as follows.. The left-hand dnj . dx and d(p/p.)/dx both tend to zero upstream, — weré 
side of Eq. (23), (d/dx)[(K — 1)Djz], becomes, when but both are of the same order, so that perat 
expanded, vA Eq 
(d/‘dx)((X\ — 1)Dz|—> 
. small 
Dj|OK O(p p.) |; i[d(p p.)/dx| & Thus, 
were 
Hence, Eq. (23) becomes the w 
) — spec ti 
d(p po) | (Me, 2 — 1)" In, R,' | p | aon 
=e _— 18 SMa 
dx Dz|OK O(p pe) |p p yMer Iw «Xx p comp! 
= right-hand side of Eq. (21) p tions. 
1.08 a 
yMer Pi I; P j 2 Ob 0.996(hy = hy) Ou (hy ms hy ) OA ( Dp ] wall t 
nas ae _ — 0.331 ad 
Xx Ty (Dz? Ol P pe) Djzly rel P pe) Ty rel P Pp. 4 t | p | fF while 
5: : , a Q.111. 
since, from Eqs. (19), a, Oa Onj, ., 6, and Ob Onzj. . > O upstream while Oa O(p p,) and Ob O(p p-) -|> 0. Thus temy 
empe 
(Vp 2 - 1)! 7 . tions, 
yMr 3 For tl 
: : re il S( 
D; jy OA 2 0b ——0.996(hw — he) Oa 0.83 (hw — he) OA j ay) fa litth 
‘ \O(p Pc) D;? O(p p-) Djln O(p p-) Ty O(p p-) f F 
The upstream limiting value of nj, , is the value for which Eq. (54) is satisfied, and can be determined from Eqs. (15 
(18), and (19). 
At the general position P downstream of the point O, Eqs. (21) and (25) give 
dW (Me 2 — 1) 7R 7/2 Me A] (bp ped (Cp pe) 1] ) ' 
= - - - om lust 
d(p/p-)  }\(2b K2D;z*?) — [0.996(hw — he)a/KDjzlw] — [0.331 (hw — he)(K — 1) lwK2]{9 
where W == (K — 1)D; 
. m 1/2- . a 1/2 ‘4: ° Now ¢ 
It simplifies matters to replace R, ~ in Eqs. (54) and (35) by R,, ~, the value at the position O. This does not in 
troduce serious errors, and it enables Eq. (35) to be integrated for II’ as a function of p /p, by a relatively simple ste] 
by-step process. Thus we can write, approximately, 
- - vet wer 
Wati = Wi + (1/2) [(/ de), — (p podali ldWV d(p pol, + [dW dip pd|n +15 it 
If IW, and [dW /d(p p.)|, have been determined, Hence 
\dW d(p) pe) |, 41 can be determined from Eq. (35) for a 
: ; ; , 1/2 
range of assumed values of nj , at the pressure ratio " dC, (MV; a. I; RV: C 
ie ; ; : = X eres 
(Pp Pe)n+1- The correct value of nz, . will be the one for dx yMer ly [dV d(p p p -iC sep 
which Eq. (36) is satisfied. Hence, IV’,,,, and [dIV- he ' elie 
: ‘ rhus, in the upstream flow, in the region of the pomt' yy" 
d(p p-)\, 4, can be determined and another step in the e 
integration performed. The starting values for IV’ and ( e 
aw d(p p.) upstream are Il, = 0 by Eq. (35), and by aii . 
Eq. (54) 31 
( Vy -— | )! ss I; R Vitl 
- P a = se oF 
dV d(p P| = Dz|OA O(p p ) | i My P ly dif ad(p Pp )] | nev b 
where this is evaluated at the limiting upstream value of This method of solution, though straightforward in | “4 
SZ, cc principle, is exceedingly tedious in practice. In the | “0! US 
The solution of Eq. (35) for nz, as a function of p p cases that were computed, the functions a, ), A, et Eq 
can be substituted into Eq. (21) to give the pressure dis- were evaluated at closely spaced intervals of 77, ., bu! Ji 
" ider t 


tribution. Eq. (21) becomes, approximately since the Eqs. (19) for these quantities involve the dil 


xX 
P/ Pe) /d 
NZ. Whi 
f NZ. - an 
osition ( 
e limiting 
OK On; 
‘he terms 
Ipstream 


S not in 


ple ster 


ward 11 
In the 
kK, ete 

i as but 
the dil 


; 
j 


i 


LAMINAR SEPARATION 


ferences of nearly equal quantities, there was some 
scatter in the estimated values when plotted against 

Thus, rounding was necessary and errors probably 
crept in. The only way of avoiding this would have 
been to evaluate the initial functions of Eqs. (15) and 
iS) toa very high accuracy at closely spaced intervals 


[he results obtained from Eq. (37) for athe ex- 


nonential rate of increase of the pressure upstream 


were 19.8, 16.8, and 12.3 for the cases with wall tem- 
neratures 245°K., 275°K., and 345°K.., respectively. 
Eq. (33), which presupposes that the inner profile is 


small in thickness, gives the values 16.9, 14.2, and 10.1. 
hus, although the computed upstream values of 7; 

were far from small, being 0.451, 0.474, and 0.499 for 
the wall temperatures 245°K., 275°K., and 345°K., re- 
spectively, the results based on the assumption that nj 

is small are not \ ery different from those of the lengthier 
computations. This is true also at downstream posi- 
tions. For instance, it was computed that when pp, = 
1.08 and C, = 0.0127, x(dC,/dx) = 0.126 and 0.102 at 
wall temperatures of 245° K. and 275°K., respectively, 
while the values corresponding to Table 3 are 0.132 and 
The computation process for the 345°K. wall- 
temperature case broke down, indicating unreal solu- 
tions, for a pressure ratio p p, slightly less than 1.08. 
For the 245°K. 


real solutions were obtainable up to pressure ratios p/p 


wall temperature, on the other hand, 


a little greater than 1.12. 


IN 


SUPERSONIC FLOW 


(2.4) THe EFFECTS OF SURFACE CURVATURI 


When boundary-layer separat:on takes place on a 
surface with convex curvature, such as a typical airfoil 
there will usually be favorable pressure gradients acting 
on the boundary layer well upstream of the region of 
separation. These will affect the state of the boundary 
layer upstream, and will thus affect the separation con 
There will also be a direct effect due 
the 


ditions indirectly. 
to the 


separating layer. 


surface curving away from underneath 
This direct effect can be investigated 
in isolation from the indirect one by an arrangement such 
as that shown in Fig. 2.. Here the surface is flat up to 
the point O, and curved, with uniform curvature, down 
stream of it, and the separation-provoking agency is so 
arranged that the pressure begins to increase at O. Ex 
periments of this sort are reported in reference 4. A 


peg is possi 


theoretical solution on the lines of Section 


ble for this case, but since it is somewhat artificial, 
further approximations are made to simplify matters, as 
follows. 

The analysis leading to Eq. (21) remains unaffected 


by the curvature. If, as in Section (2.2), the inner pro 
file is assumed to be small in thickness, and the pressure 
increases, (p/p.) — 1, are assumed to be small, relations 
(24) apply. Hence, of the terms in the braces on the 
right-hand side of Eq. (21), only the first is significant, 


and the equation approximates to 


(d/dx)(p/p-) = (yM, v) Tr, e/Tw)(2b/K*Dz*) (21”) 


Suppose that downstream of O the inctination of the 
surface to the flat wa'l upstream is given by (a Ra 


Then in place of Eq. (22), the relation 


p/p.) —1 = [yM; (CM; — 1)'*][(dy dx) — m(a v,) | 
must be used. Hence, Eq. (23) must be replaced by 
1 ‘dx)((K — 1)Dz] Mp, 2 — 1)""/yMe, 2) Ur, -/Tw)(Re'"/x) [(p pd — 1] + Ue, « Iw)(Re7/x)m(x — x) 
Now crudely, since both (A — 1)Dj and [(p/p.) — 1] become zero upstream, with vanishing gradients, 
(p pe) — 1](d/dx)[(K — 1)Dz] = (K — 1)Dj(d dx)(p p.) 
between, say, the positions O and §S of Fig. 2. Also, 
(1/2)(x — x,)(d/dx)(p/p-) (p pe) 

Hence, approximately, 

Cyl + Uw/Te, -)(K2Dj2/b)[mx/(Meg. 2 — 1)" ]h" = [87K — 1)'7/CM 1)'/*R,'/*K ](b/D;3)'"* 


\t separation, A vV/ 2 by Eq. (25). 
coeflicient at separation. 


that at separation 


Ce,,1 + 38 w/e. -)”’ “(mx Cr 


With m 
ey between the solution 
Ky | 


ol 


0 the wall is flat. The fairly small discrep- 
(3S) for * 


is due to the approximate methods of integra 


with 772 O and 


ion used in the present section. 
Eq. (38) shows that the pressure coefficient at separa- 
40n 1s reduced by convex curvature. For example, con- 


sider two cases with zero heat transfer and values of 


My 


Hence, from Eqs. (18’) and (24), nj 


This gives Dz, by Eqs. (15’), and }, by Eqs. 


compared with the value of 0.01SS for a flat wall. 


is known in terms of Cp,, the pressure 


(24), and it can thus be shown, using Eq. 


O.SS (Mer, = ik 
38) 


0.105(/,; ly )” 4 


AM and R, approximately the same as in the experi 
/ Py I 

ments of reference 4+. For J/; 3 and R 0.0 XK 
10°, Eq. (3S) shows that Cp, is 0.0105 with ma 1/3, as 


For 


Me -= 2, &. ~and mx 1/3, Cp, is 0.0150, 


O4 X& 10! 
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LAMINAR SEPARATION 


4s compared with the flat wall value of 0.0266. Of 


urse the theory is not, strictly speaking, applicable at 


i 
euch low Reynolds Numbers, but more will be said on 
this pomt 1m the following section. 


[he physical process which brings about the reduced 


REA 


essure coeflicients at separation also applies down- 
Hence the result [Eq. (38) ] would 


yl 
stream of separation. 
Jead one to infer that the general level of the pressure co- 
be reduced by 


eficients in the “laminar foot’’* should 


onvex curvature, by the order of 50 per cent probably 


Bt WIA l/s and My, 


> ] 


in the region of 2 or 3. 


) 


COMPARISON OF THE RESULTS WITH EXPERIMENT, 


» 


AND GENERAL DISCUSSION 


As mentioned before, the results obtained in Section 
2) where the inner profile is assumed to be small in 
thickness, should really only be applicable for very 
large Reynolds Numbers. 
heat transfer to the wall, the results seem to fit the ex 
perimental data quite well at low and moderate Reynolds 
Numbers in the range 0.04 X 10° to 0.6 X 10°. 
and 4 show a comparison between two sets of experi 
mental measurements? ¢ 
tions and theoretical curves derived from Table 3. The 
in the x-direction by 


However, for cases with zero 





Figs. 3 
f pressure coefficient distribu 


theoretical curves are positioned 
making the experimental and theoretical pressure co 
eflicients agree at the point 5 where the pressure co 
and where, therefore, 
theory. The 


eflicient is as given by Eq. (51), 
separation cording to the 
theoretical curves terminate downstream at the point 
Tables 2 and 5 ter- 


occurs ac 


where the computed values of 


minate. In Fig. 4, the measured pressures increase 
sharply downstream of the position x = 2.5 in., be 
cause transition to turbulent flow takes place in the 


The Reynolds Numbers in Figs. 


vicinity of this point. 
| 5and 4 differ widely and the Mach Numbers are different 

ilso, but in both cases the agreement between theory 
Fair 
the nondimensional maxi 


ind experiment is good. agreement is also ob 
tained for [x(dC,/dx) ],, 


at the upstream end of the laminar foot. 


num slope 
Fig. § 5 shows a comparison between theoretical curves 
or [x(dC, dx) |nax, according to Eq. (32), and experi 


The experimental points are very scat 


ental values.* 
tered, because they are derived from pressure distribu 
tions in the measurement of which the primary aim ws 
to determine the overall features of the flow, such as the 
accurate 


upstream extent of separation. Hence 


measurement of the pressure gradients was not at 


tempted, which in any case is a difficult thing to do. 


In view of this, Fig. 5 represents a fair agreement be 
tween theory and experiment, though it seems that the 
a little too low. 
Fig. 6 shows the theoretical curves for the pressure co- 

elhcient Cy 


theoretical curves are, on the average, 


at separation, according to Eq. (31). In 
reference 7 the rather imprecise, average experimental 
values C; 0.050, 0.022, and 0.029 are given for Mach 
Numbers of 2, 3, and 4, respectively, 
transfer to the wall, and Reynolds Numbers R, in the 
range 0.2 X 10° to 0.4 X 10°. 


with zero heat 


For the same Reynolds 
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Numbers range under zero heat-transfer conditions at a 
Mach Number of 2, experimental values of Cp, ranging 
from about 0.032 to 0.045, with the larger values occur 
ring at the lower Reynolds Numbers, can be derived 
from references 9 and 10. However, the most complete 
and consistent data for the laminar separation point are 
those presented in references 11 and 12. These results 
were obtained not by means of a surface pitot tube, as in 
references 7, 8, 9, and 10, but from observations of a very 
small quantity of oil smeared on the surface in the re 

gion of separation. <A thin hairline of oil collects at the 
separation point, and the 


thus be determined from the pressure distribution. The 


pressure p, at separation can 


experimental results show that under zero heat-transfer 
conditions at any particular Mach Number, (p pe) p 
(in the notation of the present paper) is proportional to 
the square root of the skin-friction coefficient cy, just 
upstream of the region of interaction. A relation ap 


proximately fitting all the results for Mach Numbers in 


the range 1.1 to 3.5 and Reynolds Numbers R, in the 
range 10? to 5.10° is 
(Pp Pp PeV ¢ = 

0.0206 | VW, (1, Dis ‘ly (¢ 10 
where (¢ N10 isc,;at R 10°. For the case where the 


Prandtl Number o is | and viscosity is proportional to 
absolute temperature, as assumed in the present paper, 
Cy/(Cy)195 = (10°/R 5 independent of Mach Number. 
Hence, according to the experimental results of ref 
erence 12, 

Cr, = 2(ps — pe)/ Me. 7p 


f 


ie f= O°" RS 
under zero heat-transfer conditions, a result close to the 
relation (31). Further support for this relation is given 
by a simple dimensional argument presented in ret 
1 ] 
i/2 R 1 4 


erence 12. This suggests that C; ae x R ,in 
dependent of wall temperature, at any given Mach 


and that for Mach Numbers close to 1, at 
(Mp 2 — 1)7"'*. Likewise, an 
reference 10 and 


Number, 
any rate, Cp, 
equally simple analysis presented 
based on the momentum equation suggests the result 


{ 


Ce, = 2c/7/*(Me 2 — 1)~" 
1.63(Me 2 — 1)7"” R. 


Eq. (33) of Section (2.2) can be compared with 
Lighthill's theoretical estimate'*® for the upstream ex 
ponential rate of increase of pressure with distance. 
If the viscosity is everywhere proportional to the ab 
solute temperature, as assumed in the present paper, 
Lighthill's result, which becomes asymptotically exact 


as R,,—> ©, is 


a = 0.208(Meg .2 — 1)*R,,° "(Ue -/Iw) 


> 


This is identical with Eq. (33), 
small difference in the numerical factor. 

The results obtained in Section (2.4) for the effects 
of surface curvature seem to agree quite well with ex- 


apart from the fairly 
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periment. The analysis, which assumes the inner pro- 
file to be small in thickness, as in Section (2.2), predicts 
that convex curvature reduces the pressure coefficient 
Cp, at separation, and probably reduces the pressure co- 
efficient C,7 at the top of the laminar foot in roughly the 
same proportion. The experimental values of reference 
1 for the proportional reduction in C,7 agree reasonably 
with the calculated values given in Section (2.4) for 
the proportional reduction in Cp,. The experimental 
data of reference 1+ show a similar reduction of C,7 due 
to surface curvature. 

When it comes to predicting the effects of wall tem- 
perature, however, the theory breaks down, at any rate 
at moderate Reynolds Numbers. Fig. 7 shows a com- 
parison between theoretical curves obtained from Sec- 
tion (2.2) and Table 3 and experimental pressure co- 
efficient distributions from reference 6 for cases with 
wall temperatures of 241°K., 272°K. (zero heat trans- 
fer), and 350°K. The theoretical curves have been 
terminated at x = 1.7 in., because it is evident from the 
way the experimental pressures increase that the flow 
becomes turbulent The theoretical 
curves for the cooled and zero heat-transfer cases are 
not very different, because the temperature difference 
is fairly small, and they agree reasonably with the ex- 
perimental points, but the theoretical curve for the 
It may well be, as 


downstream. 


heated wall is definitely in error. 
pointed out in reference 6, that the effect of heating or 
cooling the wall is to reduce or increase the pressure co- 
efficient at the separation point, although the overall 
pressure distribution is approximately unchanged. The 
theory of Section (2.2), as mentioned in that section, 
predicts that cooling or heating makes the boundary 
layer, in effect, more or less difficult to separate. Hence 
theory and experiment may, in a sense, agree qualita- 
tively, but they certainly disagree quantitatively con- 
cerning the pressure gradients. 

In this respect it is interesting to note that the simple 
dimensional argument presented in reference 12 gives 
results that disagree less with experiment. Reference 
12 predicts that a length /;, proportional to the distance 
OS of Fig. 1, is proportional to 6R,! where 6, is the 
displacement thickness just upstream of the region of 
interaction. Since C,, is found to be proportional to 
R,~‘''” independent of wall temperature, [v(dC, ‘dx) Jas 
should, at any given Mach Number, be proportional to 
(x 6) R, * This function is independent of Reyn- 
olds Number and decreases with increasing wall tem- 


perature. According to Eq. (32), 


[x(dC,/dx) |nax © R,/*ly- 2 


which is insensitive to moderate Reynolds Number 
variations but varies considerably more markedly with 
~. Thus, at a 
the 


wall temperature than does (x 6,)R, 
Mach Number with wall temperatures in 
vicinity of the zero heat-transfer value, the variation 


of 3 
of pressure gradient with wall temperature predicted by 
Eq. (32) is about twice as large as that according to the 
method of reference 12. Thus the latter method is less 


inconsistent with the experimental effects of wall tem- 
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perature. However, from the theoretical viewpoint. :) Sectio 
must be borne in mind that although the dimensjop,) expres 
argument of reference 12 is of great interest, becaus; uter 
its attractive simplicity, it is not rigorous. It js ioMbe ac 
instance, assumed in reference 12 that the rate of chang ere 
of friction stress with distance away from the wal] it th 
proportional to the friction stress at the wall divided }) metho 
the boundary-layer displacement thickness; here ayy) the in 
other characteristic thickness of the boundary laye! reason 
could plausibly be chosen, and the choice would aff Eve 


the results for wall-temperature effects. Also the rj) elemet 
sults of Section (2.2) of the present paper presupposf inner | 
that the Reynolds Number is very high, while the 
perimental data are for low and moderate Reynol 
Numbers. If the experiments of reference 6 could | 2/2; 
repeated at very much higher Reynolds Numbers, bet 
ter agreement might well be obtained with the predic.) 
tions of Section (2.2) concerning wall-temperatur mo’ 
effects. However, in practice it would be very difficul nalog 
to obtain in a wind-tunnel laminar flow at the very high| via 
Reynolds Numbers of 10° or more for which Section fo" 
(2.2) is, in theory, valid. ‘ee 
In Section (2.3) results were obtained which shoul a 
not be subject to the Reynolds Number restriction peng 
Section (2.2), since they were not based on the assum] eo 
tion that the inner profile is small in thickness. Hoy ™ 
ever, these results show an even bigger variation of pres- | (pl 
sure gradient with temperature than those of Sectio 
(2.2). Also they indicate unreal solutions at quite | 
values of the pressure coefficient, although experi-| and ab 
mentally laminar flow is found to persist up to much} at the 
higher pressures. Thus the analysis is unsatisfactory} this we 
and the probable reason is that there are not enough} profile 
terms in the power series in relation (16) for it to repre-} sults o 
sent the real shape of the inner profile closely when this ever, t 
forms a large proportion of the total boundary layer} would! 
However, with more terms in the series, extra conditions| in the 
would have to be used to determine the additional co-| (2.2), ‘ 
efficients. At first sight it would seem reasonable to use} bitrary 
for these conditions the relations obtained by successive | agreem 
differentiations of the boundary-layer Eq. (1) at the} differe 
wall. In incompressible flow such differentiation gives | peratu 
dicts t 
[0?(pu)/Oz*],, = 0, [04(pu)/dz*‘),, = stream 
portior 
)(1/m){O( pu) (Oz ][07(pu)/ Oxdz];,, et 
al ver\ 
and analogous, but more complicated, relations can | | robah 
obtained for the compressible case. However, in pra Reynol 
tice the use of such relations in conjunction with mor Stream 
terms in the power series would involve a prohibiti tance 
complication of the analysis. Moreover, it is by cae 
means clear that it would lead to a better representati profile 
of the inner profile. In the first place, the wall co pressib 
ditions do not specify [O(pz)/Oz],,, the most 1mportal 
differential when it comes to determining the shape 
the inner profile. Also, it seems doubtful that matte! t the 
would be improved by satisfying a large number of dil ton 
ferential relations at the wall if at the joint of the inne 
and outer solutions only pu, O(pm) Oz, and O*(pu) 0% hit 
n) Which, 


are made continuous. The outer solution obtained ! 
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1eWpoint. if Section |) could be further differentiated to give 


dimensionsMexpressions for the higher differentials. However, the 


. because yf outer solution 1s only approximate, so that there must 


It is. felhe a considerable loss of accuracy with each successive 


an.@diferentiation, and the use of higher differentials does 
the wall «i pot therefore seem justifiable. Thus it seems that the 


divided byt method is essentially restricted to conditions in which 


here the inner profile is small in thickness, and may be 
dary layel reasonably represented by a simple power series. 
ould aff Even under these conditions, however, there is an 
Iso the pf element of arbitrariness in the choice of the series for the 


presuppos) inner profile. Thus in place of Eq. (16) one might use 


iil  ¢ 
e the ou), = pw. cll, « x 


Reynol 
6 could | 3/23) + b(2?/257) + d(z*/2z*)] + 
nbers, bet (pu), . (16’) 
the predi : - ‘ » = : 7 
[his form, in fact, satisfies the compressible flow 


nl peratur : ‘ : ak 
| oe malog of the incompressible condition [0*(pu) 0z*|w = 
rv difficul 3 ' : 
hion fat VErY high Reynolds Numbers. For the reasons 
‘very hig ‘ : 


h Sect riven in the preceding paragraph, however, this does 
“A Sect 


s 


not mean necessarily that Eq. (16’) would give more 


' 
nt 
| = . 
ch shoal! wecurate results than Eq. (16) does. In any case the 
ich shoul 5 a ; ye 
» OE results given by Eq. (16) would be identical in form 
trictior ie ; ay ; 
1") with those of Section (2.2). As an alternative, one 
€ assum} 


SS. How 
On of pres | 
t Section} 


might replace Eq. (16) by the simpler relation, 
pll = pu UR « x 


wd [a(z/ 2g) + b(s* 237)] + (pu), 2/K 
h_ experi} and abandon the condition that 0°(pu) Oz is continuous 
. to much} at the join with the outer flow. It seems unlikely that 
istactory,| this would give as good a representation of the inner 
't enough } profile shape as Eq. (16), but it would again lead to re 
to repre sults of the same form as those of Section (2.2). How- 
when this ever, the use of these different inner-profile relations 
iry layer. | would lead to changes, of less than 50 per cent, probably, 
onditions }in the numerical coefficients in the results of Section 
tional co-} (2.2), and since the choice of inner profile is clearly ar- 
ble tous} bitrary, within limits, it might possibly be that the best 


uccessive } agreement with experiment would be obtained by using 
|) at the} different forms of inner profile for different wall tem 
ion gives} peratures. If this were the case, Eq. (32), which pre- 
dicts that the maximum pressure gradient at the up- 
stream end of the laminar foot is approximately pro 


ote portional to /y *) would not necessarily be true even 
tials at very high Reynolds Numbers. In fact, however, it 
1s can be) Probably is approximately correct for very large 
in pra Reynolds Numbers, because the result (33) for the up- 


ith more} Stream exponential rate of increase of pressure with dis- 
ohibitive) t€nce agrees in form with Lighthill's asymptotically 


is by nop €Xact result.!’ Further support for the use of the inner 


entation | Profile relation (16) is given by the fact that for incom- 
yall con- | Pressible flow it leads to the formula 

a sie “C, (dC, ‘dx)?], = 0.0108 (39 
shape 

matters) &t the separation point S. This is very close to the rela- 
or of dit ton 

he inne [x?C,(dC,,‘dx)?], = 0.0104 

pu) 0: 


ained in| Which, as shown in reference 15, gives the best overall 
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agreement with all the known exact and quasi-exact 
numerical computations for cases with separation in 
incompressible flow. 

There seems indeed to be a similarity between the 
incompressible and supersonic applications of Strat 
ford’s method,! insofar as in both fields results are ob 
tained which are accurate under conditions far removed 
from those postulated in the theoretical derivation. 
Thus in the incompressible case it is assumed that the 
pressure distribution is of the general shape of Fig. | 
and that the pressure gradients are sufficiently steep for 
the pressure coefficient C,, at separation to be much less 
than 1 8. 
approximately accurate under almost all conditions; 


The resulting formula (39) is, however, 


for instance, in the case where the external velocity de 
creases linearly from the leading edge and C, 0.23. 
Similarly, the results of Section (2.2) should in theory 
be valid only for very high Reynolds Numbers R,, 
greater than 10°, but in practice they agree well with the 
experimental data for the zero heat-transfer condition 
when the Reynolds Numbers are low, of order 10°. The 
agreement is less good when there is heat transfer, but 
this is hardly surprising —what is surprising is that the 
results for zero heat transfer should fit the experimental 
data for low Reynolds Numbers so well. 
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Antonio Ferri* and Martin H. Bloom** : imu 
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Brooklyn, New York eae 
May 20, 1957 tenti 
: An id 
SYMBOLS Biche 
p PIBAI 
thermal conductivity ea 
” molecular weight ; 
VJ Mach Number ul 
P Prandtl Number | 
surface heat flux 
1 absolute temperature ien 
streamwise velocity we 
streamwise coordinate ‘a 
ratio of specific heats 7 
mass density 
7 absolute viscosity | ik 
Subscripts stimat 
local flow outside boundary layer 
1 air 
H helium The 1 
stagnation value = 
quantity evaluated at surface temperature ‘ 
adiabatic surface value Howeve 
. . . F 8 
A WIND-TUNNEL TECHNIQUE for investigating hypersonic fl ipl 
phenomena at high Reynolds Numbers and with relatively] the 1 
large models is described in reference 1. The procedure in-§ 1 de 
volves enveloping the model in a shroud which is contoured would t 
such a way that the proper pressure distribution is impress be 
over the surface, or for base regions, in such a way that a suit yet 
able stream tube in the vicinity of the base is simulated. 1 Thus fi 
ssur 
+ This work was done under Contract AF 33(616)-3978 sponsored by t 
Aeronautical Research Laboratory of the Wright Air Development Center 
* Director 
* Associate Professor of Aeronautical Engineering 
For | 
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rsonic 

| relative 
cedure 

ntoured } 
impress 
hat “ 
fed T 


ored byt 


ent Center J 


o the 
tribu- 
o the 
incon- 
wing 


id, here termed the “shroud method,” is particularly 
suit the study of boundary-layer phenomena such as 
t transfer, mass transfer characteristics, and surface phe 


na. It can also be used to advantage for the study of ma 





1en 
be r and structural characteristics involving thermal 
I S method has been used in connection with the 
le heater air supply of the Polytechnic Institute of Brooklyn 
rsonic facility (references 2 and 3 From this system, 


tion pressures up to 600 psi i, St ignation temperatures up 
2.5 the 


R, and weight flows of 12.5 Ib./sec. under ifore 


ned conditions can be supplied 


DISCUSSION 
For tests involving materials and thermal effects in structures 
ften desiral 
tribution 


a pres¢ ribed 
ot 


model is 


le to obtain high heat fluxes with 
over the model. The magnitude the 


flux achievable with air for a given 


he maximum Reynolds Number that can be applied 


imum heat 


nodel and by the maximum temperature or enthalpy 
ivailable from the heating system. 

\n idea of the maximum flux available with the aforementioned 

PIBAL pebble heater can be obtained by considering the turbu 

t heat rate predicted with a constant pressure relation evalu 

it sonic inviscid conditions. For a settling chamber pres 

of 690 psia, stagnation temperature of 3,000°R and surface 

mperature of 690°R, the following relation can be used con- 

niently for this purpose since specific heat variations will have 


minor influence in this instance 


1 Taw) = 0.029P he /%) (pcttiX/pu)®* (Tf 7 )®* (1 

For a local Mach Number of 1 and x on the order of 0.5 ft., we 
stimate for air 

gq = 1,350 B.t.u. (2 

The value 1,350 B.t.u./ft.2 sec. thus represents the order of the 


heat flux available with air in the present system 


it is of interest that still higher fluxes can be achieved 


maximum 
However, 
ifa gas with a higher thermal conductivity than that of air, for 
xample helium, is used. With helium the mass flux represented 
by the term p,.; in Eq. 1 would be decreased somewhat resulting 
ina decreased Reynolds Number; however, the combined effect 
1 be to increase the heat flux by a factor of about 3. This 
1 and noting that the Prandtl 


woul 
1 be seen by referring to Eq 
ber r 


Thus for equal surface temperatures and the same stagnation 


ind viscosity w of helium and air are about the same 





pressure and temperature 


~ k p Mi ’ 
R Pi Ps) 


| ] )®-57, P87, -0-2 


Wi( 7 


ym 


For local sonic conditions the Mach Number functions p,/f 


I, /T, change negligibly from air (7 = 1.4) to helium (y7 = 
1.66) for present purposes. Likewise, the difference between 
I’,, and 7, is negligible. Hence Eq. 3 gives 

g ~~ kw[(ym 8 { 
With the following approximate ratios 
ki/Ra = 6 Vy/¥q_ = 1.2 muy/ma = 1/7 5) 
btain 
P = 3 a} 
Pywlt = p pil as ().4 7 
i the iximum heat flux available with helium is 
gu = 4,000 B.t.u./ft.? sec 8 
The comparable stagnation point heat transfer rates for p, = 
") psia, 7 3,000°R and 7, = 6OO0°R are given as follows: 
\ose radius (in.) ] yw (B.t.u./ft.* sec 
} 950) 850 
2 350 1,200 


dae 


air is used as a wind tunnel medium pri 


If 


mary emphasis would be placed not on the fluid flow character 


a gas other than 


istics as such but rather on the effects of the combined heat flux 


distribution and simultaneous pressure loading which could be 


impressed on the model. Changes in the heat and pressure 
loading distributions on a given model can be accomplished by 
appropriately varying the shroud contour 

In the PIBAL pebble heater the mass flow is limited by the 


allowable dynamic pressure in the pebble bed If helium is 
heated rather than air and if the same dynamic pressure and 
exhaust temperature 7, and pressure p, are maintained, the 


allowable weight flow with helium is reduced to 4 Ib./sec. as 


compared to 12.5 lb./sec. with air. This is of practical impor 
tance for the shroud method since a sonic throat is usually formed 
between the shroud and model when a sonic point is impressed 


model 
flows 


on the Therefore, taking into account the ratio of sonic 


mass is seen that the sonic area required to 


Eq. 7), it 
produce a sonic point in a shroud with helium is only 0.8 of the 


trea required for air. Thus, despite the fact that it 1s not too 
important to avoid channel flow in the minimum section for 
structural and material tests, it is unlikely that channel flow 


will be encountered 
Finally it is noted that the possible length of run as dictated 





by the heat transfer characteristics of the pebble bed in the 
PIBAL heater would remain the same with helium as with air 
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Criteria for the Generalized Dynamic Stability 
of a Rolling Symmetric Missile 


Charles H. Murphy, Jr.* 

Ballistic Research Laboratories 
Maryland 
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Aberdeen Proving Ground 


A S AN AID to the proper design of dynamically stable rolling 
symmetric this 


the results of the linear theory 


missiles, note will attempt to summarize 


in a rather convenient form and 
to generalize these results by the introduction of an upper bound 
for damping 

The dynamic stability criteria for a missile with both rotational 
and no thrust, are given in 


symmetry and mirror symmetry 


reference 1. For thrusts which are proportional to the dynamic 
pressure these results can be made to include the effect of thrust 
by replacing the drag coefficient wherever it appears by the dif 
ference of the drag coefficient and a thrust coefficient. Since a 
missile with its fins canted to induce roll does not possess mirror 
symmetry, a few comments are required so that the effect of fin 
cant may be included in the stability analysis of missiles flying 
at a constant spin-to-velocity ratio 

In reference 2 
of attack 


demonstration is based on a consideration of the force 


it is shown how a vaw moment proportional to 


angle and to cant angle can exist \lthough this 


normal to 


an airfoil, the inclusion of a drag force or of body-on-fin inter 


ference does not affect the form of this moment For constant 
spin-to-velocity ratio, 
C,,6 = —{¢ + k,? ( pl/ | l 


* Deputy Chief, Free Flight Aerodynamics Branch 
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where &, is the dimensionless axial radius of gyration. Since 


pl/V 
angle of attack, together with its symmetric pitch moment, can 
The 
of 


is proportional to 6, this yaw moment due to cant and 


be included in the Magnus moment term involving Cy, 


total aerodynamic moment then is that given by Eq. (4) 
reference 4: 
M+ iN = (1/2) p V? SI (pl V) Ci pe 


Cu, l\(¢g + in)l/V] — iC ,{ (él V) + 2 (pl/V) é]3 2 


} z ins + 
i Va s 


5 + ta 
With the usual assumption of equal transverse moments of 


where — = 


inertia, small geometric angles, and the neglection of all forces 
except drag and normal force,f the equation for the pitching and 
yawing motion can be written in nonrotating coordinates as 


t” + (H — iP)t'’ — (M+ iPT)t =0 (3) 
wheret 
F sie tel SP dt 
H = (pSl/2m) (Cy, — 2Co — k,~? (Cu, + Cur) 
M = (pSl/2m) [k,~? Cu. + (pSl/2m) X 
. Cy, — Cr) (ke? Cuy — Co?) 
= (pSI/2m) ki? Cu, 

T = (pSl/2m) [Cy — Cp + k,~? CH pa! 
P = (J,/I,) (pl/V) 
ki, = V1, ml? = </ I. mi? = transverse radius of gyration 


and primes denote derivatives with respect to dimensionless are- 
length. 

For constaut spin-to-velocity ratio and constant coefficients 
(linear forces and moments) the motion is described by a linear 
combination of two complex exponentials with exponential co- 
efficients: 


\; + id; = (1/2) [-H + iP + 


V4M + H? — P?+2Pi(2T-—H (4) 

We now seek the restriction imposed on the coefficients in Eq 
(3) by an upper bound A on the damping exponents. (The re- 
sulting generalized dynamic stability criteria reduce to the usual 


dynamic stability criteria when J is taken to be zero.) 


A <A AM<KA (5) 


H+2rx>0 6) 


H+ 20> |R[VW4 + H? — P?+ 2Pi 2T — || (7) 


Inequality (7) can be simplified by routine algebra using the 
t+ For quite conservative estimates, the other aerodynamic forces should 
have no effect on the stability of the motion 
t The second term in the definition of M is important for the stability of 


[(pSl/2m) ™ 1, Cn > 0 


airships which are statically unstable ie 


AL SCIENCES OCTOBER 1957 
identity reg 
R[VW4A+ iB]! = \ | + V/A? + B?| 

Pp? — 4M P? ((2T — H)/(H + 2d)]? — 4A/H \ ; 
For \ = 0 this is precisely Eq. (24) of reference 2 w] the 
forces considered in that equation are the normal force and dro t 
force. Laitone® has stated that this equation in refcrenc 
erroneous due to the use of an incorrect square root approxi 


mation, Since Inequality (9) is exact, this does not seem t 


the case 
missiles the second term on the right 


For heavier than air 


Inequality (9) can be neglected and with the introduction 


generalized dynamic stability factor, sa, the following conve 
form may be obtained: 
$1//P? 2 — sa)s 
where Sa = 2(T + A)/(H + 2d 
The term on the left of Inequality (10) is the reciprocal of : 
gyroscopic stability factor The condition that spinni 


missile, acted on by only the static moment (Cy,_), will “sle 
- a 
on its trajectory is that this quantity be less than unity 


The implications of Inequality (10) are concisely displayed 


Fig. 1. (The choice of the reciprocal of the gyroscopic stabilit Som 
factor as the ordinate was suggested by W. R. Haseltine WI ns 
Inequality (6) is satisfied, the requirements of dynamic stabilit Sa 
may be summarized by the use of this figure: 
(1) A dynamically stable missile must be gyroscopically stat ‘ie 
(2) If sy lies in the interval (0,2), a statically unstable missi sis 
may be stabilized by sufficiently high spin and a statically stabk} . jytii 
missile is always dynamically stable a 
(3) If sa lies outside this interval, a statically unstable miss deities 
cannot be spin-stabilized and a statically stable missile m 
be made dynamically unstable by a sufficiently high spin nlate-1 
It is interesting to note that most of the ways a symmetr oe 
missile can become unstable have been observed.’ | ” 
j 
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Limitations of the Finite Difference Method 
of Structural Analysis 


Bertram Klein 

Head, Structures Theory and Methods Group, Convair Astronauiic 
San Diego, Calif 

May 24, 1957 


T THE May Issue of this Journal, Weingarten has presented “g 
note! showing and advocating the use of the method of finit C 
differences augmented by extrapolation in solving certain t Jourr 
angular plate instability problems. The contention is that this 
method may be superior to the method presented by the autho! 
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READERS’ 


First let it be 





gard 1 mplicity, accuracy, and generality 
tl e author already has used this approach in the 
ferred to by Weingarten There is no doubt that the 
life e technique is a powerful tool in the realm of 
itu vsis. However, as is true for all methods, it has 
in di yacks and shortcomings. For example, the method 
ite ences does not always give the correct answer as 
plific i recent paper by Bijlaard Then again, the ac 
rtifice of extrapolation is not always valid as pointed 
recent paper by Bijlaard.* Certainly, when the 
lems are nonlinear, the finite-difference method may become 
cumb me Another fact to remember is that each prob 


isa different numerical example it is not possible to solve, 
method, a set of geometrically dissimilar plates with one 
Also, different types of meshes are required for treating 
ructures of The 


sensitivity of the matrix defining a plate problem by the 


the same type but geometrically dissimilar 


v be a limitation 


In certain problems of this type, 
of points needed to give the first approximation to 


number 


solution may become excessive It has been found that on 
‘aleulations in 
In fact, 


the worse the condi 


the matrix is ill-conditioned—41.e., the 


ve the taking of small differences of large numbers 
some problems, the more points taken, 
ing Of the matrix 


Some engineers have raised objection to the method of finite 


They claim that little engineet 
On the other 


ifferences on a different basis 
judgment is required in using the method 
nd, in the method of collocation, there is room for good engi 
insight A knowledge of the 

sical behavior of the structure is an important part of the 
Phe 
solution in which relatively few points of collocation 


in solving a problem 


istute choice of deflection functions and regions 


ntion 





required. For example, Weingarten! discusses the solution by 

differences of the simply supported equilateral triangular 
plate under uniform compressive edge thrust acting on all edges 
The author has solved this problem by the method of collocation 
using hybrid functions® and only two points which yields a much 
reference 2. The buckling 


ser answer than that found in 
pe is taken as: 
= |a@ F a F. cos | rxh 2ya 
F, = sin (ry/h) sin (ry/2h F, = sin (2ry/h) sin? (ry/2h 


The two points chosen for collocation are along the y axis in 


region of maximum plate deflection and are 


v = (2/3)h and y = (3/4)h 
rhe buckling coefficient is found to be 39.32 which is less than 1/2 
recent away from the correct solution? on the safe side. Of 
urse, the iccuracy Can be improved, if deemed necessary, by 


ising more points 
One objection raised against the methed of collocation is that 
me does not know whether a calculated buckling load is too high 
» low he author feels that eventually some augmenting 
tistical or similar method of predicting this location as well as 
redicting or increasing the accuracy of the results obtained will 
ve 
In conclusion, it should be realized that, although the method 
finite differences is a fine method, it is subject to certain limita- 
ind disadv intages, and therefore, the engineer should keep 


pen mind as to what other methods are available and can be 











veloped further, or to possible new methods of structural 
Sis 
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Further Comments on the Buckling of 
Triangular Plates 
Hugh L Cox 
Executive Vice-President, Structure Specialties, Inc., Santa Monica 
Cal ¢ 


June 3, 1957 
A RECENT NOTE by Weingarten! calls att 
for 


use of extrapolation formulas in conjunction 


ention toa well-known 
solving certain buckling problems by the 
difference 


technique 
with finite 
methods. Weingarten computes the critical buckling load for a 
simply supported uniform equilateral triangular plate subjected 
the 


boundaries of 


to uniform axial compressive loads along the 


plate, and then makes several statements that are somewhat 
misleading. He states that the finite-difference technique is 
simpler, more accurate, and more versatile than the one pre 
sented by Cox and Klein;? however, his note contains no new 
results or solutions that may help aircraft engineers. If thi 
finite-difference method is so inexpensive and easy to use it 
would seem that some new curves or improvements of those 


given in reference 2 should be presented 


The author questions seriously the possibility of obtaining a1 
1 ] g 


accurate critical buckling load with only a few grid points, as 
used by Weingarten,!' when the plate geometry is such that 
several buckles may occur For example, when //a $) 
reference 2, Fig. la) the true buckled shape of the plate when 


represented by a few grid points will not fit accurately into the 
fourth-degree polynomial which represents a first-order difference 
pattern 

The purpose of the paper presented by Cox and Klein? was to 
provide a reasonably accurate chart for estimating buckling loads 
of isosceles triangular plates of various shapes. Since the solu- 
tions presented involved only computations associated with a 
three-by-three determinant, the authors believed that the collo 
cation solution, for the range of parameters considered, was con 
siderably simpler and faster than solving many different systems 
by finite-difference equations which would lead to considerably 
of 


higher order determinants than three-by-threes for values 


(h/a) > 2.5 

The example! of an equilateral triangular plate involves fewer 
simultaneous equations than another triangular plate, since many 
of the grid points of the equilateral triangular plate have the 
same deflections as a result of symmetry 

The author agrees that the finite-difference technique may be 
considerably more versatile for many other problems than the 
finite-differ 
ence solutions does not always work out as Wein- 
a plate or beam of highly vari 


collocation method; however, extrapolation with 


as accurately 
garten apparently assumes. If 
able stiffness and load distribution is considered, some question 
arises as to how much difference exists between assumed analyti 


cal models for different grid sizes \ plate with highly variable 





stiffness may require that many grid points be taken; thus, 
accurate extrapolation from small-order systems will not be 
possible 
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Nonsteady Surface Temperature Effects on 
Forced Convection Heat Transfer 


E. M. Sparrow and J. L. Gregg 
NACA Lewis Flight Propulsion Laboratory, Cleveland, Ohio 


lune 5, 1957 


INTRODUCTION 


N THE COMPUTATION of forced convection heat transfer from a 
I surface whose temperature varies with time, it is common to 
suppose that quasi-steady conditions exist. That is, it is assumed 
that the actual heat transfer at any instant is identical with the 
steady-state heat transfer corresponding to the instantaneous 
surface temperature 

Further, in the experimental determination of heat-transfer 
coefficients by the transient technique, it is customary to assume 
that the instantaneous measurements may be used to compute 
steady-state values. So, the assumption of quasi-steady heat 
transfer is again invoked 

In reality, there will always be some difference between the 
actual instantaneous heat transfer and the quasi-steady value 
The extent of the deviation is determined by the ability of the 
temperature and velocity distributions to respond to the changes 
in surface temperature. So, the response characteristics of the 
boundary layer and the rapidity of surface temperature varia 
tions both play a role 

The aim of the present study is to compute the first-order 
deviations of the instantaneous heat transfer from the quasi 
steady value. The system chosen for analysis is a semi-infinite 
flat plate (see Fig. 1). 
to be spatially uniform over the surface, but may vary in an 
arbitrary (but The 
stream velocity is steady, and compressibility is included 


The temperature of the plate is taken 


differentiable) manner with time free 
Those 
who are primarily interested in results are invited to pass over 
the section on analysis [see Eqs. (15) and (16)}]. 

Interested readers are referred to the papers of Moore! and 


Ostrach,? who studied departures from quasi-steady conditions 


due to variations in free-stream velocity. 


ANALYSIS 
The equations expressing conservation of mass, momentum, 
and energy for unsteady laminar boundary-layer flow over a 
flat plate are 


(O/Ox) (pu) + (O0/Oy) (pv) + (Op/dt) = VU (1) 


p[(Ou/Ot) + u(Ou/Ox) + v(OUu/Oy)] = (0/Oy¥) [w(Ou/Oy)] (2) 
pc,((OT/dt) + u(OT/dx) + v(OT/dy)] = 
(0/Oy) [R(OT/Oy)] + w(Ou/dy)? (3) 
Time is denoted by ¢ and the static temperature by 7 
It has been pointed out by Moore! and Ostrach? that the con- 
servation-of-mass equation is satisfied by a stream function de- 


fined by 


u=(p p) (Oy Oy) | 


v = —(p,/p) [(Ov/dx) + (0/dt) | (p/p..)dy| \ 
J0 


Introduction of the following new variables 
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Fic. 1. The mathematical model 
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g@=(T—T Te - 


I / / ; 
Y) - Wry —=— wrxr¢ = vy 9 
6, + O[Tu/(Tw — 7 YyOx — yx 
l Py 19, Tt iv k 1 — 7] y¥ } 


X¥, Y, and 7 


represents the time derivatiy 


For brevity, the derivatives with respect to 
noted by subscripts, while 7 
the wall temperature. The fluid properties used in the 


vation of Eqs. (2a) and (8a) are pul constant, pk const 


c, = constant 
The problem considered here is not one where we specify jnj 


seck 


Rather, we focus attention on a particular y 


conditions (at ¢ = ©) and then solutions for subs 


values of time 
of time and inquire about the departure of the actual inst; 


ous conditions from quasi-steady conditions. With thisinn 


we seek solutions of Eqs. (2a) and (8a) having the form 
y¥= VviUX lfoln) + B n) + B n) + | 
6 = Alyn) + BAN) + BxB2(y) 4 
+ [U.2/2cp(Te — T,)|¥(n)) 
where 
sn = (Y/254) 7 (24/1 R, == [7 ,-— 1 V/l 
= |7 | ae | X/U.)? 
The variable 7n is seen to be the Blasius similarity variable, 


the parameters 6), (2, are a measure of the promptness wit 
which the boundary layer responds to impressed changes. WI! 


the Bn are very small, the velocity and temperature profiles 


esseiitially quasi-steady 
Introducing the series expansions (5) into Eqs. (2a) and 
gives the following differential equations for fy, @, 4, 42, and 4 


The functions fi, fo, are found to be identically zero fort 


present problem. 


6 + Pr f,6 = § 
A” + Pri fi)’ — 2f)'0; — 46) = 0 
6.” + Pri fiOe’ — 4f,'0. — 46,) = 0 
®” + Prlfob’ + (1/2) (fo")?] = 0 


The primes indicate differentiation with respect to 7 


boundary conditions that « = 0,7 = 0, T = T, at y = Oandt 
u= U,,T =T, aty = transform to 
f, =f,’ =0,6=0 
n 0 

>= 1,04 = & = = () 

f =2 @= 0 

6, = 6; = & = oer) 
The differential Eqs. (6), (7), and (10) are identical to the g 


erning equations for the steady-state problem. The san 


true of the boundary conditions. So, the solutions of t 


equations are available in the literature, and, in particular, rest 
for Pr = 0.72 are given by Ostrach.? 
were obtained ! 


Numerical solutions for Eqs. (8) and (9 


Pr = 0.72 as part of this investigation 


RESULTS 


The instantaneous local heat-flux rate, dinst., at the plate s 


face may be calculated using Fourier’s Law: g = —[RO7/Oy 
Introducing the series expansion of Eq. (5) and the dimension 


variables of Eqs. (5a) and (4), there is obtained 


where 6)'(0),6,'/(0), .. . are abbreviations for [d@)/dyn|, =0, 


are found from solutions of Eqs. (7) through (10 By firs 


oe 


rents 
stanta 
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the adiabatic wall temperature, 7,,,, from the relation 


] Toi:1 +R — 1)/2|M,,?} 
} — |’ (0) /8, () = yl 2 
the icing that the quasi-steady heat transfer is given by 
R(T. — 7 U op /Va%)/* } [—05'(0)]/2$ (18 
to rewrite Eq. (12) in the form 
/ XY “7 () 7 
P ( 
( lei l 4,0 
] \ 21 @.°(0 | 
beg 
7 — |] [ a,/.0 { 
g numerical results for P?7 0.72: @'(O —().5913 
1.416, A. (0 1.4739; there is obtained 
j . 1 \ 
l 2a) : = 
| T,— 7 l 
z 1 ( \ | - 
SO1L5 + 15 
Tr, —- 7 l ( 
ese results may be given in an alternate form by defining 
1 heat-transfer coefficients as 
Ty. — 1 } = ¢ 7, — 7 


O.S8015 


The form of Eq. (15) is convenient for computing heat transfer 
from surfaces having time-dependent temperatures, while Eq 
16) is in a form useful for interpretation of heat-transfer coeffi 
Whether 


stantaneous heat transfer is greater or less than the quasi-steady 


cients obtained under transient conditions the in- 
value appears to depend upon tle particular values of 7, Tx, 

, and T, — f, As an 
mple, suppose 7, = 10°F T. = —0.5°F . iw = 
Tow = 5O°F., U, = 300 ft x = 1 ft. For these 
onditions, the instantaneous heat transfer deviates negligibly 
from 


which enter a given problem eXx- 


] sec sec 


sec., and 


» 


the quasi-steady value, the difference being less than 0.2 


per cent 
If the analysis were carried out neglecting frictional heating, 
16 


through would still apply with the modification 


that (7, — 7 is replaced by (7, — T 
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Sidewall Influence on Oblique-Shock 
Pressure* 
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INTRODUCTION 


_ SHOCK WAVE EQUATIONS for the treatment of oblique 
shocks in the absence of significant sidewall frictional forces 


well known However, the accuracy of the customary 


* This re University of Illinois in 1954, at 


Mechanical 


earch was carried out at the 
which time Engi 


the author was a member of the staff of its 


ering Department 


RS’ 


ad 


FORUM 


analytical relationships is sometimes questioned when thick side 
wall boundary layers are present in an approaching internal flow 
Nevertheless, it can be anticipated that such boundary layers 
would exert very little influence on the shock relationships if their 


small 


fied here both theoretically and experimentally 


mass and momentum defects are This conclusion is veri 


PHEORETICAL ANALYSIS 


An approximate analysis of the influence of sidewall boundary 
by 


and momentum principles, for a perfect gas, to t 


layers on shocks may be made applying the continuity, 


energy, he control 


domain. Such an analysis has been carried out! for the similar 


‘normal’ shock problem and the following is essentially an ex 
McLafferty’ 


Phe geometry of the flow is shown in Fig. 1 


s theory to the oblique-shock pro 


As the fundamen 


tension of 


tal relationships are icross the domain, informa 
flow det shock 


boundary-layer interaction is not required, and the simultaneous 


ipplic d 


garding the iils in the neighborhood of the 


equations connecting state 1 and state 2 are easily obtained n 








rder to make the problem tractable, a quasi-one-dimensional 
uniform velocity profile is assumed at state 2. This dismisses 
the variables which describe the downstream boundary layer 
On this basis, the two simultaneous equations becom«e 
, , 
Io) 1 1 )/2] Ata? */- 
l 7 May 
: , 1 
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Fic. 2. Boundary-layer influence on pressure ratio 
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Fic. 3. Shadowgraph of airflow (wo = 10°, Mia ~ 1.83, 


6/a =~ 0.51) 


where the subscript a identifies conditions as those in the isen 
tropic core adjacent to the boundary-layer flow. 6, 6*, and 6** 
are, respectively, the boundary-layer thickness, the boundary- 
layer displacement thickness, and the boundary-layer momentum 


thickness. Thus, 


and 
w = fo[M,,, m, (6:/h), 0 (4 


Eqs (3) and (4) result from Eqs. (1) and (2), since 6*/6 and 6**/6 
are dependent only on the Mach Number adjacent to the 
boundary layer and to the shape of the boundary-layer velocity 
profile. Tabulated values of these ratios as a function of 7, and 
n of an “nx” power law-velocity profile are available 

A relation for the oblique-shock pressure ratio which is con 


venient for calculation purposes is 
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Fic. 4. Surface-pressure survey (turbulent boundary layer, 
R, = 5.7 X 10°, h = 0.150 in., 6/h = 0.51) 
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P: Vinatl + [Cr 
Pp Me,{1 4 1)/2] Mz} 1/2 
For practical considerations, Eq. (5) may be placed 


form, stated succinctly as 
P2/P, = f3|.M,, m, (6:/h), w 


Typical calculated results for air are shown in Fig. 2. The pr 
sure increases experienced by flows with turbulent sidey 
boundary lavers are seen to be almost equal to those obt ul 
without boundary-layer influences Indeed, the calculatj 
show that had a fully developed turbulent channel flow existed 
state 1, the pressure ratio across the shock would only be about 
per cent less than that obtained for the ideal case of zero bou 


ary-layer thickness 


EXPERIMENTAL RESULTS 


Experimental data of the pressure history through an oblig 
shock with a thick approaching boundary layer were obtained ir 
small rectangular (0.3800 in. by 1.260 in.) supersonic chann 
The shock was produced by a 20° (total angle) wedge. Tot 
pressure profiles were taken perpendicular to the side walls 
several locations, without the wedge in the stream, in order tod 
termine the boundary-layer thickness throughout the channe 


A shadowgraph of the flow in the neighborhood of the wedge 





shown in Fig. 3. Note that even though the approach boundar 
layer is thick, a fairly ‘‘clean”’ pattern appears 

The pressure history through the shock was obtained from fixe 
multiple sidewall pressure taps. These were located parallel t 
the major axis of the wedge at a distance (“0.10 in.) believedt 
be sufficiently remote from the disturbed field due to the ff 
curvature near the apex of the wedge and the disturbed field du 
to the Mach reflection at the nozzle block walls. The wedge w 
translated along the axis of the channel in order to obtain pres 
sure data at small length increments. The results of the pressur 
survey are presented in Fig. 4 

The difference between the actual pressure rise and that cal 
culated neglecting the boundary layer is less than 2 per cent 
surface pressure tap, positioned on the geometrical center of tl 
wedge face, provided the face data which, as indicated, vari 
slightly while the wedge was translated. The length of the 1 
teraction zone is approximately four boundary-layer thicknesses 
An analogous ‘‘normal” shock in a circular channel would have at 
interaction zone length of from 11 to 16 boundary-layer thick 
nesses.” 

The drop in pressure is due to the expansion fan. It starts 
ahead of its theoretical position, because of the boundary layer 
along the wedge surface. Naturally, it would have been ad 
vantageous to obtain the data free from the effects of any expan 
sion. However, this could not be done with the available test 


tunnel 
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‘‘Vortex Interference Effects. . . 


M. Z. v. Krzywoblocki 

Professor, Department o Aeronautical 
Illinois, Urbana, Ill. 

June 20, 1957 


Engineering, University 


—— VERY INTERESTING PAPER by Sacks! on the vortex inte! 
ference effects on the aerodynamics of slender airplanes aud 
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missiles does not contain any discussion on the possible rigorous 


-olutions of the nonlinear equations pertinent to the problem in 
Although the existing solutions are still mathematically 


quest! mn 


omplicated, nevertheless the high-speed computing devices may 


to use them in the future. Certainly, the application 


nable oO! 
f them will give a more general analytical review of the prob- 
The author calls attention to reference 2 
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Experimental Pitot-Pressure Profiles of 
Supersonic Laminar Boundary Layer in a 
Weak Adverse Pressure Gradient* 


Eugene A. Racicot 
Research Assistant, Institute of Aerophysics, University of Toronto 
June 13, 1957 


HE EFFECT of adverse pressure gradients on laminar boundary 
T layers is of theoretical and practical interest. The aim of 
this investigation was to examine the pitot pressure profiles in a 
laminar boundary layer under a weak adverse pressure distribu- 
tion 

The work was carried out in the 5- by 7-in. supersonic wind 
tunnel at the Institute of Aerophysics. The tunnel has an at- 
mospheric pressure inlet and exhausts into an evacuated sphere 

A pair of nozzle blocks providing a nominal Mach Number of 3 
were mounted toward the upstream end of the test section. The 
side walls and floor of the section are flat, but a special ceiling was 
installed which was a continuation of the upper nozzle block, and 
curved slightly inward toward the tunnel centerline. From this 
inward curving liner a family of compression waves crossed the 
flow to be reflected at a flat plate which spanned the tunnel 

The testing surface for the laminar boundary layer was the 
upper surface of this flat plate. The leading edge thickness was 
0.0017 in. The plate was supported from the floor of the tunnel 
by a central longitudinal support 

The first profiles were taken with rectangular probes. The 


from the 
In the 


profiles so obtained differed from those calculated 
analysis of G. M. Low? by factors of up to fifty per cent 
light of the concurrent work by G. R. Ludwig? on spurious probe 
effects associated with rectangular probes of conventional sizes, 
these measurements were rejected and circular probes were made 
and used. They were 0.0075 in. outside diameter and 0.0064 in. 
inside diameter. These were the smallest probes which could be 


successfully manufactured in the UTIA shop 


*The work summarized herein was supported by the Defence Research 


Board of Canada. A complete account is contained in reference 1 
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The pressure distribution over the model was obtained from 
two sources. One was the direct measurement of the static pres 
sure over the model by means of pressure holes in the surface, and 
the other was a static pressure inferred from the stagnation pres 
sure of the flow and the pitot pressure measured when the probes 
were sufficiently far from the boundary layer to be free from 
boundary-layer effects. The pitot stations were not at the same 
longitudinal positions as the static pressure holes; thus, the num 
ber of stations at which the static pressure was determined was in 
creased. The results are shown in Fig. 1, where pf is the static 
pressure; subscript 7 refers to the reference conditions (taken as 
those of the incident flow) and x is the distance along the plate 
from the leading edge. The values obtained by the two methods 
were in fair agreement. A second order polynomial was found to 
give a good fit to the experimental longitudinal pressure distribu 
I 


tion; this gave the coefficients required for the application « 
Low’s theory 

As the probes used were all of the same diameter, the ratio of 
probe height to boundary-layer thickness increased at those sta 
tions which were nearer to the leading edge. The errors due to the 
probe interference effects were quite small with the circular probes 
when the ratio of probe height to boundary-layer thickness was 
less than 1/4 

The results obtained for one of the test stations are shown in 
Fig. 2. Here h is the probe height, y is the vertical distance from 
the surface of the plate to the physical center of the probe, and 
H is the pitot pressure. Two theoretical curves are shown. The 
solid line is the result obtained by appiving Low’s theory to the 
experimental conditions, with the pressure distribution along the 
The dashed 
line is the result obtained by assuming that the pressure and Mach 
Number outside the boundary layer are constant along the plate 


plate assumed to be given by the solid line of Fig. 1 


and equal to the experimental local values at the measuring sta 
tion, and then applying the Chapman-Rubesin‘ theory for laminar 
boundary layers under no pressure gradient 


Several conclusions can be drawn from this figure. (1) The 
pressure gradient applied was sufficient to produce a noticeable 
difference between the two theoretical results. (2) The predic 


tions of Low’s theory are in good agreement with the experimental 
results. (3) The ‘“‘quasi-local’’ calculation in which the zero-pres 
sure-gradient theory is applied gives an approximation to the 
pitot-pressure profile which, for the applied pressure gradient, 


is in error by about 10 per cent 
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T° REFERENCE 1, experimental results were presented which indi 0 4 8 12 16 20 hos La 
cated that although moderate cooling in cone-cylinder and DISTANCE FROM TiP OF CONE ~ INCHES : ge A 
parabolic-cylinder models resulted in an increase in transition > 7 p tices , 1.3 B Heat Tri 
; : iG: 2 mperature distr on on » 20-in. 10° cor 
Reynolds Number, extreme cooling reversed this trend and re . ” regen pit xi ‘ re saa d a boa as : : Me Ve 
- é iS : bs in steady state with extreme cooling and with zero heat transfer. f 154 
sulted in a decrease in the transition Reynolds Number. The de- VWs = 2.70 : 
4 ane . ° ° “ Jac 
lay in transition was predicted theoretically in references 2 and 8 ies 
and was well-documented experimentally for moderate cooling in temperature ratio, T/T, of 0.60 Fig. 1 shows the variatior NA 
references +-6. In reference 4, experimental reversal in transition with time of the temperature of a monitored thermocouple 6 
was obtained, but only for the case of sufficient roughness, and it in. from the tip of the cone during progressive cooling. The 
I g g 1 
was argued that the effect was probably due to the thinning of the model was cooled internally with a mixture of gaseous an 
laminar boundary layer relative to the surface roughness liquid nitrogen. Included on the figure are schlieren pictures 

Experimental data obtained recently by the present authors (magnified 20 times normal to the flow) taken at each 4 see. in Taste 
on a 20-in. 10° (apex angle) smooth (10 micro-in.) cone tested terval. The pictures show clearly the movement of transitior Pre 
in the 12-in. tunnel of the Jet Propulsion Laboratory at a toward the rear of the model. As cooling progressed, the transi 
supply turbulence level of 0.3 per cent, a local Mach Number tion moved off the model, leaving a completely laminar boundary 

6 - - . o ~ ¢ ss VW 
Mo of 2.70, and a local Reynolds Number per in. of 0.70 XK 10°, laver. It is noted that the change in transition, as observed by R 
° mar ’ “ 4 Nortnrc 
do not show a reversal in transition Reynolds Number when the magnified schlieren, lagged somewhat behind the change i 9 
5 e i ’ 2 He ne 2 
the boundary layer is subjected to extreme cooling such that the monitored temperature at the lower temperature ratios. This si 
the wall-to-total temperature ratio, T/T, lies in the neigh- can be explained by the fact that temperatures downstream of the 
F 25 naeiiees : 5 ee eS ; . ; Tr Is 
borhood of 0.25 (corresponding to a wall-to-adiabatic-wall monitored temperature also lagged behind the monitored tem- | | : 

. ws ‘ cy . Oo7 3 ; : the 
temperature ratio Ty/Tadia of 0.27 and to a wall-to-local perature, owing to turbulent flow downstream where heat Numt 

- ° \umbe 
transfer rates were greater. Indeed, simultaneous measurement ; : 

t This research was supported by the United States Air Force, through the of the temperature at the 6.8 in. station showed a time lag in J ae 
Office of Scientific Research of the Air Research and Development Com reaching steady state of approximately 4 sec. behind the tempera : rhe 
mand 5 7 iret bys 5 : a tee - Reynol 

© Clot Scientist. Technical Sciences ture at the 6.3 in. station. After 44 sec., a steady-state tempera- | ‘inti 

** Senior Engineer, Technical Sciences ture distribution existed along the model This steady-state Sul : 
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Fic. 1. Delay of transition on the 20-in. 10° cone with extreme cooling. 1/5 = 2.70 * Sen 





READI 
f ributic f temperature during extreme cooling is presented 
© cthe lower data in Fig. 2. It is seen from this distribution and 
; rs corresponding magnified-schlieren picture that with extreme 
oling (7 ./Tr = 0.25) the flow was completely laminar, (For 
x mparisol, the steady-state zero-heat-transfer temperature dis 
tribution and its cor responding magnified-schlieren photo are also 
shown in Fig. 2 There, both temperature distribution and 
10000 © shlieren photo indicate transition from laminar to turbulent 





70 x 108 fow.) During the cooling process, it was not until after 3 min. of 
running time that a hazy frost began to appear on the model 

In conclusion, experiments conducted by the present authors do 
not show a reversal in transition Reynolds Number with large 
mounts of cooling, contrary to the results obtained in reference 1 
Rather, it appears that a delay in transition may still be ac 


mplished at very low temperature ratios. 
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it [’ IS THE OBJECT of this note to demonstrate the importance of 
eC em- ° . ° 
heat the usually neglected base pressure variation with Reynolds 
S Ca 


aaa F Number when a turbulent boundary layer exists at the base and 
t Cll . - . 
to present a means of accounting for this effect 
e lag in a he . : , ; , 
f rhe variation of the base pressure of a cylindrical body with 
empera- § er , 7 
noun — Reynolds Number shows a fairly rapid change when the bound 
Yin peta- ' *,¢ . . 
aye iry-layer transition passes forward from the wake to the body 
SLA . °,° . 
Subsequently, as the transition point proceeds forward on the 


body, the boundary-layer thickness at the base increases and the 
base pressure increases correspondingly. The normal decrease in 
local turbulent boundary-laver thickness with increasing Reynolds 
Number eventually begins to offset the thickening effect of the 


forward movement of transition and as shown in Fig. 1, the 


base-pressure curve reaches a local maximum. Any further it 

rease in Reynolds Number causes the base pressure to decrease 

relatively slowly. This latter variation is often considered negli 

gible. Data obtained over small ranges of Reynolds Numbers 

have tended to confirm this hypothesis. However, the seeming] 

negligible slope of the curve of base pressure with Reynolds Num 

ber can become very significant if the Reynolds Number changes 

pF are sufficiently large to cause appreciable variations in the 
boundary-layer thickness 

Consideration must be given to the large Reynolds Number 

Variations which can exist under several important circumstances 

§ For example, the Reynolds Number can decrease by a factor of 

‘pproximately 70 in traveling from sea level to 100,000-ft. alti 


tude. The importance of this effect becomes apparent in the ap 
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Fic. | 3ase pressure vs. Reynolds Number 


plication of flight-test data to other flight conditions. Possibly 
more significant is the error introduced when wind-tunnel data 
are used to define base pressures in flight. Since the Reynolds 
Number is only a measure of the boundary-laver growth, any 
change in the effective fineness ratio between the model and the 
full-scale flight article can cause an important change in the thick 
ness of the boundary layer relative to the base diameter A 
boundary-layer trip can have such effect. Application of the data 
obtained on such a model can lead to serious errors in drag esti 
mation or afterbody design 

Chapman® has correlated base pressures of cylindrical after 
bodies (with turbulent boundary layers) and the function L/D 
(Re)!°, where L/D is the body fineness ratio and Re is the 
Reynolds Number based on body length. This correlation func 
tion is proportional to the ratio of boundary-laver thickness to bast 
diameter. In addition, Chapman has noted that the distribution 
of density and velocity through the boundary laver are also pet 
tinent It therefore seems more appropriate to correlate the 
base pressure with the function 6*/D which is the ratio of the 
boundary-layer displacement thickness to the base diameter 
The displacement thickness, 6*, is evaluated at the base and is, on 
a simple body, directly dependent upon Chapman's correlation 
function as well as on the local boundary-layer profiles. In terms 
of 6*/D, the base pressure varies monotonically throughout the 
entire turbulent regime 

The data correlated by Chapman are at Mach 1.5 and 2.0 
Over the range presented, they show the base pressure increasing 
linearly with his correlation function and the slope of this curve 
increasing rapidly as the Mach Number decrease Phe slopes at 


) 


the two Mach Numbers are presented in Fig. 2 in a modified form 
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Fic. 2. Boundary-layer correction to base pressure for evli 
drical bodies with no exiting jet 
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\. Lambda is the rate of change of the base-pressure coefficient sheets was shown by the writer‘ to be a maximum of per cep fm (and s 
with respect to the normalized boundary-layer displacement for c/t equal to 10, 7 per cent for c/t equal to 20, and negligible | rres] 
thickness c/t of 50 or greater for the problem of a plate with si nply s 
X = dC,,/d(8*/D ported loaded edges and clamped unloaded edges. It can be 
pected that similar discrepancies would hold for the pres er 
The curve drawn through the two data points in Fig. 2 is con problem 
sidered representative of the variation of X with Mach Number It should also be pointed out that the sandwich plate inyey | 
At present it is assumed that A ts independent of the displacement gated by the author is one with a core that has isotropic 
thickness. The base-pressure increment due to the change in properties, that is, the shear modulus is the same for both + B fhe pa 
boundary-layer conditions is then simply the product of \ and yv-s and y-2 planes. Thus, end-grain balsa and cellular celly 
the difference in 6*/D between the actual and simulated flight acetate core sandwich plates can be analyzed by the results of th. : 
ae 2 os : ' i 
conditions. It should be noted that the change in 6*/D is only paper, whereas metal honeycomb and corrugated core sandwic | 
about 20-30 per cent of the corresponding change in boundary plates cannot (or can with decreasing accuracy, if an averag 
laver thickness-to-diameter ratio, 6/D, at Mach Numbers below shear stiffness for the two directions is used and the actual s 
3.0. It is expected that \ can vary appreciably for bodies with stiffnesses differ markedly 
: ; a 3 : é one ; jis i 
various degrees of boat-tatling and for bodies with an exiting jet i : 
= , : : ; : Bot! 
Since application of this boundary-layer correction can produce REFERENCES 
; age s ‘ erime! 
changes in the base-pressure coefficient which are of the same ini ES Be Bixee ad Mee Bines Clintised Saudis - 
order of magnitude as the coefficient itself, it is recommended that Journal of the Aeronautical Sciences, Vol. 24, No. 6, pp. 407-412, June. 1957 ' = 
" ° - . . ‘ » > f Derk Je > . it 1s 
an effort be made to obtain data which would fully delineate these Hoff, N. J., Bending and Buckling of Rectangular Sandwich Plates, NACAB | 
effects rN 2225, 1950 Revno 
Bijlaard, P. P., Stability of Sa ich Plates, Journal of the Aeronaut cept 
| " 3 Sciences, Vol. 16, No. 9, pp. 573-574, September, 1949 t If we 
<EFERENCES ‘Seide, Paul, Compressive Buckling of Flat Re Metalite ; : 
7 } >] th D upphorted led / + ci tiie 
Kurzweg, H. H., /unterrelationship Between Boundary Layer and Base Sandwich Plates with Simply Supported Loaded Edges and Clan 1 Unloa , 
Pressure, Journal of the Aeronautical Sciences, Vol. 18, No. 11, pp. 743 748 Edges (Revised), NACA TN 2637, 1952 the po 
November, 1951 itible 
Crocco, I and Lees, L 1 Mixing Theory for the Interaction Between In re 
Dissipative Flows and Nearly Isentropic Streams, Princeton University Aero ¢ ecthy 
ctl 
nautical Engineering Lab. Report No. 187, January, 1952 
Chapman, D. R 1x Analysis of Base Pressure at Supersonic Velocitte dai 
and Comparison with Experiment, NACA TN 2137 1950 Pai?“ 
; si — a er ee oe On Turbulent Flow and Eddy Heat Transfer 
° “ee . ° ' T 
Diffusivity in a Pipe? 
inal 
+ Raul R. Hunziker* and José S. Florio** In hi 
Mathematical Physicist, Research Division, Reed Research, In Laufer 
. Washington, D.C., and Research Associate, The Johns Hopkir 
sé 1 3 , ' . - near th 
A a - é ye ipa Buckling of University, Baltimore, Md., Respectivel, 
Ont 
amped Sandwic ates July 19, 1957 : 
be defi 
c vabl 
Paul Seide N HIS ANALYSIS of the flow in a circular pipe! or between parall -" 
Guided Missile Research Division, The Ramo-Wooldridge 2 Pai eee bi, sic 
plates? Pai considers simplifying hypotheses which allow hi > 
Corporation, Los Angeles, Calif é i id ; bis ag" < e Re 
to write the evnolds equations in the form first given Db 7 
June 17, 1957 ; ; ees 1 & vall (n 
Kampé de Fériet® for parallel plates i 
. a ° : As in the case of the parallel plates, the equations for the circula tril 
HE RECENT PAPER! by G. A. Thurston contains the useful re . ed . ‘ 3 ributi 
3 , Sfeae : : : : cs pipe reduce to a pair of independent equations which are insu > 
sults of a quite complex investigation which might be modified ime : : . , epres¢ 
é : ficient to determine uniquely the five related dynamical magn oe 
slightly to be more theoretically correct. It should be noted that rae in the | 
: ys ; -_ ; : tudes characterizing the flow baad 
Hoff’s theory? applies to a sandwich plate with a core which hypo : , : mS O 
; : ; dia One of these equations is tite a1 
thetically extends to the middle surface of the cover sheets. For ae 
P . ° P ° I » 
actual sandwich plates in which the core is bonded to the inner pu'w'(r) = tolr/a) + pdi(r)/d bie 
surface of each cover sheet, the change required in Hoff's energy aoe 
expression is merely to change the expression for the strain energy where Phe 
of shear in the core to read Pe ‘adi mditi 
a pipe radius = 
. . . 9 , = ‘ . -ly > 
eas Ge fj u (t/2)(Ow/Ox) ow |° r = radial coordinate solutio: 
> JJ c/2 ox p = dene) . = el " Acco 
. P uu’, wW = radial and axial fluctuation of velocities fnsinets 
. » op , as 1° ° . —— : 
i (t/2 = Oy) Ou | dudy u’w'(r) = correlation of velocities - 
c/2 Ov f pu’w'(r) = Reynolds ‘turbulent stress 
Gi Oe Ee *y ou ow 7? T = shen stress at the wall ieee: 
= 1 + i + m = absolute coefficient of viscosity From | 
2 ( c+ Ox : : : are : ‘rom 
“-s wr) = time average axial velocity distribution 
2v ow |" een - “ | 
- - dxdy Considering the boundary conditions 
crt oy { ‘i 
ee - " ee wea) = 0, u'w(a) = 0 - . 
Pherefore, in the final results, it is merely necessary to replace ere 
the core hear modulus G of Thurston s paper by G(1 t c) to Pai has found the form of a polynomial approximating 17(/) whi 
obtain results for sandwich plates with cores bonded to the inner reduces to the parabolic distribution of Poiseuille when x‘ 
surface of each face plate. It should be noted that this redefinition 
will reduce the value of S calculated for a given plate and will i 
thus lead to slightly higher calculated buckling loads. ee re ee ore 
TI . lo : | F pial P sponsored by the Office of Scientific Research of the ARDC, USAF, unde By i 
a reqeres change in the strain cuere) of shear in the core Contract AF 18(603)-104 
was first pointed out by Prof. P. P. Bijlaard.* The discrepancy * Now with Convair, A Division of General Dynamics Corporation, 4 erall 
between results for the core bonded to the middle surface of th a ca nd 8 
yetween results for le core DONAdCC( © the middie surface of the ** The authors wish to acknowledge valuable discussions with Gordor : 
face sheets and the core bonded to the inner surface of the face C. K. Yeh, Staff Scientist, Reed Research, Inc., Washington, D.¢ =e 
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REA DI 
nd suc it its radial derivative at the wall matches with the 
rresponding shear stress. This polynomial is 

rn mae n ; 
ere 
= "ihe Con*", n= ( 
C, =(s n)/(n 1), Cc. = (1 s)/(n l } 
e parameter s 1s related to the shear stress at the wall by 
s = rna/2wi(O 5 


e distribution of x ‘w"( also has a polynomial form 


, u/ pa iO )fo(n 
n Cyn | 1 n= , C; = 2n(s ] n l 6 
dis identically zero for the case of laminar flow in which s = 1 


Both distributions contain a parameter # which is adjusted ex 


erimentallv, and satisfy jointly Eq. (1) in the whole field (0 < 


< 1) including the so-called laminar sublayver region 

It is noteworthy that in the solution of Pai, the corresponding 
Revnolds shearing stress remains finite in the laminar sublayer 
xcept at the wall itself, where it vanishes 

If we consider the Reynolds equation, the condition of nonslip, 
nd the vanishing of the Reynolds stress at the wall as postulates, 
e polynomial distributions of Pat for any value of ” are com 
patible with them 

In relation to this mathematical consequence which follows 
lirectly from the Reynolds equation and a nonslip condition at 
the wall, Wieghardt* has stated that the velocity distribution of 
Pai? 


the corresponding turbulent shearing stress remains finite in the 


“isan approximate formula and not a rigorous one since, e.g., 


yninar sublavers except at the wall itself.” 

In his recent experimental] investigation in turbulent pipe flow, 
Laufer§ found that high activity of turbulence occurs in a region 
near the so-called laminar sublayer 

On the basis of these empirical facts the laminar sublayer may 


ve defined by the condition y’w’ S yo where is the lowest ob 
servable value of the correlation of ‘Velocities. All these assure 
the physical plausibility of Pai’s analytical description in which 
the Reynolds turbulent shear stress is assumed zero only at the 


, and finite in the laminar sublaver 
of 


tributions which satisfy Eq. (1 


ill(mn = 1 


To each pair values , s, corresponds a pair of Pai’s dis 


and the boundary conditions (2 
representing the observable magnitudes with an approximation 
in the mean which can be optimized. In this sense, the distribu 
tions of Pai constitute an “approximation,” but aiso they consti 
tute a mathematically rigorous solution of Reynolds Eq. (1) and 
Eq. (2) even if it is not the only possible solution, because other 
polynomials, by addition of other terms, may be solutions also 

satisfy these 


of 


The logarithmic velocity distribution does not 
mditions at the 
Reynolds Eq. (1 


solution in the inner lavers 


wall; hence it is not a rigorous solution 


and Eq. (2) even if it may be a possible part of a 


According to the ‘‘momentum transfer theory’’® the eddy dif 
tusivity for heat e,, is defined as 
en(r) = asy(r) = —an'e'(die/d 7 
lere €y, is the eddy diffusivity for momentum 
From Eqs. (3), (4), (6), (7) it follows 
é l (n ( = 
C1 n* (2C, 2nd ~ Ss 
ere vy is the kinematical coefficient of viscosity 
trom Eqs. (1) to (5) it follows 
= 1)/2 9 
ind from the data given in Fig. 4 of references for RK, = 
Za (O)/, 50,000 it results s = DISS 


32 gave a good 


By inspection it was found that the value » = 
erall fitting to the empirical functions represented in Figs. 3, 4, 


na S of reference 5. The values of €y(7)/v obtained for the case 


= 50,000 by application of the relationships used in references 








FORUM 783 
Fic. 1 Different momentum diffusivities for R 50,000 
7, 8, 9, 10, 11 were calculated in reference 11 and are shown in 
Fig. 1 

These values are generally higher than those given by Eq. (8 
for R, = 50,000 and represented by Curve A of Fig. 1 

The eddy diffusivity for momentum, as derived by von Kar 
man,!? is given by 

eu(n)/vy =n | ” ] 10 


and is represented by curve B of Fig. 1 (which almost coincides 


with curve A The application of Eq. (10) to the logarithmic 
velocity profile gives negative values of €),()/v near 7 1, which 
have been excluded and replaced by zero by Marris rhe same 
distribution of €y())/v in the turbulent core was used by Mar 


tinelli.. This function is represented for R, = 50,000 by curve C 


of Fig. 1 


The distributions €,())/v indicated in references 9 and 10 are 


represented as curves D, E of Fig. 1, respectively 


The low values of a = € 4 /€y 1 which reference 13 indicates 


for liquid metals, were considered'* to be in agreement with the 
fact that the theoretical predictions of Martinelli’ and Lyon 
1 give high unrealistic values of the local heat 


based on a@ = 


transmission coefficient and Nusselt Numbers 

Since the eddy diffusivity for heat is ¢, = aey, higher values of 
the local Nusselt Number may be expected, not only if @ is large 
but also if ey is large. Taking a fixed value of a, this fact may 
explain that generally the calculated values of heat transfer rates 
and Nusselt Numbers from the walls of the pipe were higher than 
the experimental ones 

With Eq. (8) and the classic hypothesis a = | a new and better 
fitting to the experimental Nusselt number data has been ob 
tained 


The use of Pai’s polynomial distributions allows a consistent 


determination of the eddy diffusivity for heat e, aéy on the 
basis of the “‘momentum transfer theory” and with a minimum 
of assumptions 
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Note on the Mach Number Effect Upon the 
Skin Friction of Rough Surfaces 


Hans W. Liepmann* and Frank E. Goddard, Jr.** 
California Institute of Technology, Pasadena, Calif 
July 19, 1957 


— “SKIN FRICTION” Of fully rough surfaces is the cumulative 
form drag of the protuberances At a sufficiently high 
Reynolds Number the drag coefficient of bluff bodies becomes 
nearly independent of Reynolds Number and Mach Number. If 
A denotes the frontal area of a roughness element and » the num 
ber of elements per unit area of surface, we can expect the drag d 
per unit area to vary as 

d~ purX(A-n 
where p and uw are density and velocity in the neighborhood of the 
element. If we now compare the skin friction coefficients for 
similar roughness in compressible and incompressible flow, we can 
write 

(py/ pri) [(u/U1).2/(u/ U1),? 
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where p; and l are density and velocity in the free stream ag 
p Pwall 

Consequently the principal Mach Number effect appea 
through the density and hence temperature ratio. This sugge 
instead of Mach Numbe 
Fig. 1 shows recent measurements at the Jet Propulsion Lab 
tory of the California Institute of Technology plotted in thi 
p, or 7\/T,), is eviden 
Thus as a simple rule for the effect of Mach Number on the dra 


plotting the ratio C,,./C,; versus 7)/T. 


manner; the linearity of C,./C;; with p 


of fully rough surface we have 


Wr} = const 


where r is the recovery factor. At extremely high Mach Numbe 
additional effects due to acoustic radiation may become i 
portant 

Thus the density effect, which is largely offset by viscosity & 
fects for laminar skin friction and partly offset by viscosity & 
fects for turbulent skin friction on smooth plates, is com 
pletely dominant for rough surfaces. Hence the decrease in 
friction with increasing Mach Number is most pronounced wh 
the surface is fully rough. Full details of the experimental stud} 
of skin friction of rough surfaces at supersonic speeds will be found 


in a forthcoming paper by Frank E. Goddard, Jr 
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